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ABSTRACT 


This  thesis  Is  part  of  a  more  general  research  program 
1 n vest iga t i ng  the  failure  process  of  a  tunnel  in  a  brittle 
jointed  material.  For  this  purpose  a  model  testing  technique 
on  natural  material  (coal)  is  used.  Two  specimens  have  been 
tested  for  this  particular  thesis.  One  of  the  tests  was 
carried  beyond  the  rupture  of  the  tunnel  wall. 

The  test  results  and  experimental  conditions  are 
reoorted  but  th<*  analysis  Is  concentrated  on  the  study  of 
the  time  —  independent  mechanisms  involved  in  the  rupture. 

Closed— form  solutions  for  the  isotropic  e 1 as t o— p la s t i c 
analysis  of  the  tunnel  closure  measurements  have  been 
adaoted  to  the  model  test  conditions.  Different  hypotheses 
for  the  post— failure  and  volume  change  behaviour  have  been 
considered.  Comparisons  with  experimental  data  have  showed 
that  an  Isotropic  t i me  —  1 n de pen de n t  e lasto— plast 1c  analysis 
cannot  explain  the  failure  process  in  details*  even  if  it  is 
consistent  with  the  elastic  and  strength  properties  of  the 
ma+erlal.  It  is  suggested  that  an  isotropic  elastic  — 
anisotropic  plastic  analysis  be  used*also  taking  into 
account  the  time-dependent  and  time  -independent  anisotropic 
stress  redistribution  processes. 

The  rupture  of  the  tunnel  wall  in  the  second  test  is 
described  and  commented.  Joints  in  the  material  have  been 
proven  to  play  a  role  in  the  rupture  process  by  intercepting 
and  shortening  the  plastic  slip  lines.  The  rupture 
Initiation*  propagation  and  eventual  stabilization  have  been 
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pxpla inpd  by  a  model  taking  into  account  the  joint  density 


and  +he  local  radius  of  curvature  of  the  tunnel  wall*  Local 
instabilities  of  the  tunnel  wall  could  generate  collapse  of 
a  tunnel  and  should  be  analysed  by  an  anisotropic  plastic 
equilibrium  analysis* 
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RESUME 


Cette  these  est  incluse  dans  un  programme  de  recherche  plus 
general  destine  a  etudier  le  processus  de  ruine  d1 un  tunnel 
dans  une  roche  fragile  et  jointi ve.  A  ces  fins,  on  utilise 
des  modeUs  reduits  de  tunnel  dans  un  materiau 
na+areK  cha^bon )•  Deux  blocs  dif ferents  ont  ete  teste  pour 
cette  these  ;  1 1 un  d*eux  presentant  une  ruine  du  tunnel* 

Les  resultats  de  ces  essais  et  les  conditions  dans 
lesquelles  ils  ont  ete  effectues  sont  deer  its  mals  1  'analyse 
crncerre  essentlellement  I'e^ude  des  mecanlsmes  independants 
du  +emt>s  lmpliqu@s  dans  la  ruine  de  l  'ouvrege* 

La  convergence  du  tunnel  a  ete  etudiee  par  des  analyses 
elasto-olastiques  et  isotropes,  adaptees  pour  le 
circonstance  aux  condi  tons  oarticulleres  du  tnodele  t'pdui  t» 
Ces  analyses  classlques  conslderent  dif ferentes  possibilites 
pour  le  eomportement  de  la  roche  apres  se  rupture  ou  en  ce 
qui  concerne  les  variations  de  volume*  Une  etude  comparative 
avec  les  resultats  e xper i ment aux  a  montre  que  ce  type 
d' analyse ,  independante  du  temps,  ne  peut  expliquer 
enmpl  etement  les  deformations  menant  a.  la  ruine  de 
l' oavrage,  meme  si  elle  peut  expliquer  indi vi duel lement 
cheque  convergence  avec  des  parametres  d'elasticite  et  de 
resistance  consi stents  avec  les  essais  de  laboratoire.  Une 
analyse  isotrope  dans  le  domaine  elastlque  et  anisotrope 
dens  le  domaine  plastique  devrait  etre  utilisee  pour  ce 
genre  d' etude  ;  elle  devrait  aussi  considerer  les  processus 
complexes  de  redistribution  des  contralntes  qui  sont 
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anisotropes,  et  qui  dependent  ou  non  du  temps* 

La  rupture  de  la  paroi  du  tunnel  dans  le  second  bloc 
est  dec  rite  et  commentee*  On  a  montre  comment  le 
jointolement  du  milieu  jouait  un  role  important  dans  cette 
rupture  en  interceotant  et  rac c ourc i ssa n t  les  l lgnes  de 
^1  isseme  nt  plastiques*  L'amorce  de  cette  rupture,  sa 


d i f  f us  ion  e  t 

sa  stabilisation  finale  sont  expllquees  par  un 

mode l e  d  on  t 

les  princioaux  parametres  sont  la  densite  du 

jo intpiement 

et  le  rayon  de  courbure  de'  la  paroi*  Des 

i n  stab  ill  tes 

locales  de  la  parol  oeuvent  enrendrer  la  ruine 

cfo  1  1  o'.ivi'ase  et  doivent  etre  etudiees  par  une  analyse 
anisotrope  de  1'equilibre  plastique* 
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CHAPTER  1 


INTRODUCTION 

Gepepp _l 


The 

last  fifty  years 

have 

witnessed  tremendo 

us  pr 

ogre  ss 

t  he 

cor  s  tr uc  t i o n 

of 

stable  tunnels# 

The 

ma  j  o  r 

contributions  are  known  as  the  New  Austrian  Tunnelling 
Method  (NATM)  end  the  Obs er va t i one l  Tunnelling  Method#  New 
technologies  such  as  rock  bolting  and  shotcrete*  new 
construction  procedures  involving  rate  and  sequences  of 
excavation  and  subsequent  supports  with  different 
stiffnesses*  monitoring  of  deformations  and  loads  and 
subsenuent  adjustment  of  the  design  are  the  main  features  of 
these  successful  methods# 

More  recently  the  availability  of  advanced  computer 
techniques  has  permitted  theoretical  analyses  of  some  of  the 
processes  involved  in  tunnel  stability#  Although  an 
imDortant  contribution*  these  analyses  are  now  over  — 
sophisticated  and  need  to  be  supplied  with  observations  on 
the  actual  mechanisms  involved  in  the  time  —  dependent  and 
time  -  independent  failure  of  an  opening#  For  instance 
stress  red i s tr i butl on  due  to  isotropic  or  anisotropic  creep* 
or  to  global  or  local  yielding  around  a  tunnel  are  rarely 
taken  into  account  in  analyses  because  the  processes 
Involved  are  not  well  understood# 

Kalseri  1979  )  initiated  a  research  program  at  the 
University  of  Alberta  in  order  to  investigate  the  time  — 
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dependent  deformations  around  a  tunnel  in  a  brittle 
material*  For  this  purpose  he  chose  to  develop  model  tests 
rather  than  to  analyse  actual  case  histories  (which  have  the 
disadvantage  of  being  rarely  documented  in  sufficient 
details)*  He  used  coal  as  a  model  material  because  this 
natural  material  was  readily  a va 1 le ble , had  suitable 
structural  characteristics  and  had  already  been  studied  by 
Noonan( 1972)*  Even  if  it  has  disadvantages  such  as  lack  of 
control  of  uniformity,  the  use  of  a  natural  material  has  the 
essential  advantage  of  supplying  the  test  with  the  intrinsic 
properties  of  a  true  scale  natural  material  •  Time  - 
dependent  rupture  mechanisms  may  thus  be  controlled, 
monitored  and  recorded* 

The  purpose  of  the  model  testing  undertaken  in  this 
research  is  not  to  simulate  the  behaviour  of  an  ec+ual 
prototype,  but  to  understand  the  failure  process  of  a  tunnel 
in  na+ural  material*  Kai ser( 1 979 )  studied  in  detail  the  time 
—  dependent  deformations  of  the  tunnel  before  failure 
occurred*  This  present  thesis  analyses  the  time  — 
independent  processes  of  rupture*  Future  work  will 
investigate  the  time  —  dependence  of  this  rupture  process* 

i*.2  Behaviour  oj_  the  Eocjc  Mass  near  an  Qnenine  la  a  Brittle 

sLalated  Medium*. 

Excavation  of  a  tunnel  generates  a  tangential  stress 
concentration  In  the  vicinity  of  the  wall*  Because  the  creep 
properties  of  most  rocks  are  also  stress  deoendent,  stress 
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on 

occurs 

near 

the  tunnel  where 

the  confinement 

i  s 

s  t  • 

This  cau 

se  s 

a  time  —  dependent 

redistri but  ion 

o  f 

t  o 

zones  fu 

rt  he 

r  from  the  tunnel 

wa  11* 

a  t 

a  ny  t  i 

me 

the  stress  level 

exceeds  the  time 

t  strength  of  the  rock,  yielding  occurs*  These 
deformations  may  be  global  or  only  local  and  are 
cause  for  stress  redistribution  to  occur*  After 
deformations,  the  tunnel  wall  may  rupture*  This  is 
ed  by  the  structure  of  the  yielding  rock  element  and 
most  likely  local*  Th is  is  a  third  reason  for  stress 
bution  to  occur*  The  r ed 1 s t r i but i on  of  stress  may  be 
c  but  it  is  most  likely  anisotropic* 

pagation  of  the  plastic  zone  end  propagation  of 
terminate  eventually*  Where  the  confinement 

s  (away  from  the  opening)  the  material  becomes  more 
and  the  stress  redistribution  process  becomes 
d*  The  boundary  conditions  of  rock  elements  involved 
ruoture  process  become  less  favourable  for  a  local 
ity*  The  recognition  of  these  stress  redistribution 
s  is  essential  for  a  proper  evaluation  of  the  load 
nel  support* 


!«_3  Scope  of  this  Xheala 
The  objective  of 
independent  mechanisms 
opening  in  a  brittle 
more  general  research 


this  thesis 
involved  in 
rock  mass*  Th 
program  and 


is  to  study  the  time  — 
the  stability  of  an 
is  thesis  is  part  of  a 
for  this  reason,  the 
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testing  procedure  (  mu  It i ole — s t age  creep  tests)  is  similar  to 
the  one  used  by  Kaiser(  1879).  The  time  -  dependent  data  are 
reported  and  described  but  the  analysis  is  concentrated  on 
the  rupture  mechanisms*  Time  —  dependent  aspects  of  rupture 
will  be  analysed  in  a  further  study*  For  the  same  reason r 
testing  nrocedure  and  equioment  are  only  broadly  described 
in  this  thesis*  More  details  have  been  already  given  by 
Ka  i  se  r(  1  97  9  ) . 


Ch  ap  t 
vised  in  t 
by  Noonan 
pressure 

con  t r i but i 
this  mater 
Test! 
Chapter  3 
equ i omen  t 
Two 

mon  th  s  te  s 
Model  test 
The  resul 
B  1  . 

Chap  t 
study  the 
measuremen 
hypotheses 
been  adap 


er  2  summarizes  the  knowledge  about  the  material 
he  model  tests*  This  coal  has  been  already  studied 
(  1972)  t Ka i se r (  1 97 9  )  and  da  Fon t our a(  1 98 0  )  •  High 
triaxial  tests  are  reported  to  constitute  the 
on  of  this  thesis  to  extending  our  knowledge  of 
i  a  l  • 

ng  procedure  and  equipment  are  described  in 
•  More  details  regarding  the  development  of  this 
are  given  by  Kai ser(  1979  ) 

coal  specimens  have  been  tested  during  a  eight 
ting  period*  They  are  called  MC— 3  and  MC  — 4  (  MC  for 
on  £oa l  )•  These  tests  are  described  in  Chapter  4* 
ting  stress  —  strain  curves  are  given  in  Appendix 

er  5  analyses  the  data  near  failure*  In  order  to 
deviation  from  linearity  of  the  tunnel  closure 
ts,  elas  to-»plastlc  analyses  with  different 
of  post— failure  and  volume  change  behaviour  have 
ted  to  the  model  test  conditions*  Appendices  Al  to 
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A4  give  the  comolete  derivations  of  these  methods,  a  list  of 
symbols  and  other  miscellaneous  items*  Theoretical  and 
experimental  data  are  compared  and  the  discrepancies 
analysed*  In  the  same  chapter  an  attempt  is  made  to 
determine  Independently  elastic  parameters  of  coal  from  the 
model  test  data* 

Sample  MC-4  experienced  a  rupture  of  the  tunnel  well* 
The  description  of  this  rupture  and  of  the  mechanisms 
involved  In  Its  initiation,  propagation  and  eventual 
stabilization  Is  the  objective  of  Chapter  6*  Time  — 
dependent  data  of  the  rupture  are  described  in  the  same 
chapter  and  reported  in  Appendix  B2* 

A  summary  of  the  conclusions,  some  implications  and 


recommenda t i on s  for  future  research  constitute  the 


concluding  Chapter  7* 


, 


CHAPTER  2 

COIL  AS  A  MQDEL  MATERIAL 

2x1  Intro due t log 

To  use  a  natural  mat 
material  has  similar  inconve 
any  geotechnical  engineer 
heterogeneous*  anisotropic 
controlled  but  have  to  be  det 
The  purpose  of  this  cha 
origin*  structure  and  mech 
Results  f»*oni  several  triax 
pressures  will  be  summarized 
used  in  the  model  tests  will 
This  chapter  mainly 
Noonan(  1972)  *  Kalser(  1979) 

determination  of  the  prope 
explanations  will  be  found  in 

2x2  Origin  and  &triic.tjire 

2x2x1  Origin  and  geology 

The  coal  samples  were  collected  at  the  Highvale  Mine 
located  on  the  south  shore  of  Wabamum  lake*  75  km.  west  of 
Edmonton.  This  strio  mine  provides  coal  for  the  nearby 
Sundance  power  plant  operated  by  Calvary  Power  Ltd. 

Large  blocks  were  collected  in  the  mine  and  samples 
were  trimmed  from  these  blocks  for  triaxlal  tests  as  well  as 


erial  such  as  coal  as  a  model 
niences  to  those  encountered  by 
in  field  work.  The  material  is 

and  its  properties  cannot  be 

e  rm  i  n  e  d  • 

oter  is  to  describe  briefly  the 
anical  properties  of  the  coal, 
lal  tests  with  high  confining 
and  the  structure  of  two  blocks 
be  de  scr i bed • 

summarizes  the  work  done  by 

and  da  Fontoura(  1980  )  in  the 
rties  of  the  coal.  Detailed 

the  above  references. 
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for  model  tests  (  Ka  i  se  r(  1  97  9  )  )  • 

Wabamum  Lake  is  located  on  the  boundary  between  the 
massive  sandstones  of  the  Tertiary  Paskapoo  Formation  to  the 
west  and  the  Upper  Cretaceous  Edmonton  Formation  to  the 
east.  The  coal  is  in  the  Scollard  member  of  the  upper  Ardley 
coal  zone,  at  the  top  of  the  Edmonton  formation,  which  is 
one  of  the  most  promising;  sources  for  Albertan  coaKHolter 
Q t  al( 1975)  )•  This  coal-bearing  stratum  overlies  different 
tyoes  of  bentonitic  shales.  Pearson(  1959  )  describes  the 
Wabamum  Lake  district  in  detail. 


2±.2l±.2.  Physical  Properties 

The  coal  is  class! fie 
to  the  Canadian  class! ficat 
range?  between  21%  and  24% 


( 1.58)  and  the 

degree  o 

provided  by  Noona 

n(  1972  ). 

Most  coals 

are  or+h 

bedding  planes 

and  two  or 

Pomeroy(  1966  )  , 

Ko  an  d  G 

coal  contains 

horizonta 

d i  sc  on tl nuous  , 

closely  sp 

and  a  poorly  de 

vel oped ,  i 

first  cleat  system  will  be 


d  as  "Sub-bituminous  B"  according 
ion.  Its  natural  moisture  content 
•  Values  for  the  specific  gravity 
f  saturation  (70%— 100%)  have  been 

otropic  in  structure  , containing 
thogonal  cleat  systems  (Evans  and 
erstle(1976)  )•  The  Hlghvale  mine 
l  bedding  planes,  a  vertical, 
acec?(  0. 5-2  cm)  first  cleat  system 
rregular  second  cleat  system.  The 
called  a  joint  system  hereafter. 


2. 2. 3  Alteration  of  the  Physical  Properties 
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? . 2 »  3 « 1  In  troduct Ion 

It  will  be  shown  at  the  end  of  this  chapter  that  the 
two  blocks  used  for  model  test  purposes  have  a  different 
structural  aspect*  One  block  contains  more  mi cro— f ractures 
than  the  other*  However,  they  have  been  collected  and 
handled  together  In  the  same  manner  before  being?  stored*  The 
only  difference  Is  in  the  storage  time  and  in  the  sample 
preparation*  These  observations  create  a  concern  to 
determine  the  reasons  for  the  alteration  of  the  physical 
properties  of  the  coal,  in  order  to  optimize  the  storage  and 
samole  preparation  conditions* 

When  extracted  from  its  natural  environment  a  coal 
sample  may  suffer  two  major  types  of  alteration  : 

1*  moisture  content  change 
2  •  ox i da  t i on 

These  processes  have  been  described  by  Fryer  and 
Szladowl 1 973 )  •  This  section  mainly  summarizes  this  work* 

2  *  2 »  3  *  2  Changes  in  moisture  content 

In  its  natural  state,  coal  is  more  or  less  saturated 
with  water*  When  exposed  to  elr,  the  moisture  content 
reduces  until  a  new  state  of  equilibrium  is  reached*  If  the 
relative  humidity  of  the  environment  increases  ,  the  coal 
will  reabsorb  moisture.  This  cycle  may  be  repeated  many 
t I mes • 

As  the  moisture  transfer  is  gradual  from  the  surface  to 
the  Interior  of  the  coal  sample,  internal  stresses  are 
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generated  by  differential  expansion  and  cause  a  partial 
disintegration  of  the  structure.  As  well  os  causing  this 
alteration,  which  is  called  "decrepitation"  or  "slacking", 
absorption  of  moisture  generates  heat  which  accelerates  the 
oxidation  orocess. 

2jl2j.3  *.3  Oxidation 

While  losing  and  absorbing  moisture  when  exposed  to 
a*r,  coal  will  also  "chemisorb"  oxygen  and  so  become 
oxidized.  This  process  may  lead  to  a  molecular  degradation. 
It  is  a  function  of  the  exposed  surface  and  the  temperature. 

Precautions  to  take  in  storing  coal  samples 

Generally  low  rank  coals  will  slack  much  more  quickly 
and  extensively  than  coals  of  higher  ranks.  To  quantify  this 
phenomenon  a  slacking  index  was  defined  as  the  proportion  of 
—9.265  in  material  formed  when  a  +1.05  in  coal  sample  is 
alternatively  air  dried  and  rewetted  by  immersion  in  water 
(Fryer  and  S la zdow(  1 97 3  )  ) •  This  index  was  statistically 
correlated  to  the  natural  moisture  content  of  the  coal.  The 
coal  used  for  the  tests,  with  an  average  moisture  content  of 
24%  ,  has  a  slacking  index  of  50%. 

It  is  important,  when  storing  samples  to  minimize 
dry i ng— we t t i ng  cycles  and  oxidation.  This  can  be  achieved  by 
storing  the  samples  under  conditions  of  constant  humididity 
or  more  simply  under  water,  keeping  the  temperature  as  low 
as  possible,  but  above  freezing. 
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2«  2.  3.  5  Storage  ancf  Handling  of  the  Model  Test  Blocks. 

The  coal  samples  have  been  stored  in  a  moisture  room 
with  100%  relative  humidity  and  at  a  temperature  of  5°C.  To 
prevent  any  moisture  loss*  the  model  test  blocks  have  been 
protected  after  trimming  by  a  latex  coating,  wrapped  in 
olastic  and  strapped  with  plywood  sheets.  Moisture  content 
determinations  on  samples  stored  in  such  a  way  for  a  year 
have  shown  that  this  method  prevents  any  significant 
al terati on • 

When  in  the  compression  machine,  the  samples  are 
covered  by  a  3  mm  plaster  of  Paris  layer  and  a  sandwich  "wax 
pa oer-plas te r— wax  paper",  and  the  tunnel  wall  is  covered 
with  a  double  coating  of  latex.  In  the  room  where  the 
experiments  are  carried  out,  the  relative  humidity  is 
unfortunately  not  kept  constant.  However,  determination  of 
moisture  content  on  the  core  obtained  when  drilling  the 
t unn  e 1  and  on  the  sample  at  the  end  of  the  test  are 
similar.This  shows  that  the  method  used  to  protect  the 
sample  is  efficient. 

The  samples  will  suffer  much  more  during  the  transient 
stages  of  the  handling  process: the  original  transfer  from 
the  field  to  the  laboratory,  the  trimming  and  the 
pr epa  r a  t ion. 

2 . 2 . 3 . 6  Conclusions 

It  appears  that  repetitive  drying-wetting  cycles  have 
to  be  avoided  in  handling  coal  samples.  The  best  solution  is 
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to  kepp  the  original  moisture  content*  Storage  conditions  in 
the  moisture  room  are  sufficient  to  achieve  this  protection 
against  moisture  loss  for  a  long  period  of  time.  It  is 
essential  however  to  keep  the  sample  protected  at  all  times 
during  the  transient  stages  described  above*  The  most 
efficient  solution  is  to  keep  the  sample  wrapped  in  plastic 
sheets  and  to  uncover  it  tor  minimum  time  periods*  The 
practice  of  ’'watering"  the  sample  during  the  preparation  »  In 
order  +o  reduce  moisture  loss,  has  to  be  avoided  as  it  may 
be  detrimental  for  the  sample  structure*  This  may  have 
caused  the  "decrepltat  ion'1  effect  In  sample  MC  —  4* 

2  •_  3  Meehan  leal  Properties 

Review  of  Previous  Testing 
2*3*  1  *  1  Strength  Properties 

Direct  shear  tests  on  both  precut  and  Intact  samples 
have  been  conducted  by  Noonan(  1972  )( see  also  Morgenstern  and 
Noonani  1974  )  )•  The  main  conclusions  can  be  summarized  as 

follows  : 

1*  For  a  dlscont inuously  open  joint,  the  angle  of  friction 
0  associated  with  the  peak  strength  depends  only  upon 
the  shearing  resistance  of  the  intact  material  and  the 
cohesion  Intercept  increases  with  decreasing  degree  of 
separation  of  the  joint*  The  degree  of  separation  is 
defined  as  the  ratio  between  the  total  area  of  the  plane 
containing  a  joint  and  the  area  of  open  joint  surface* 

2*  The  Youne's  modulus  ,  back-calculated  from  direct  shear 


■ 


12 


test  results*  ranges  between  140  and  550  MPa*  increasing 
with  normal  stress* 

3*  The  ultimate  0— value  of  the  strength  envelope  is  30°* 

Kalser(  1979)  conducted  further  direct  shear  tests  with 
normal  stresses  between  1  and  4  MPa  and  concluded  that  there 
was  a  reasonable  agreement  between  his  results  and  Noonan's 
conclusions*  Assuming  a  linear  Coulomb  failure  envelope  and 
a  friction  angle  of  30  0  *  cohesion  intercepts  were  c  a Icul a  ted 
between  0*87  MPa  and  1*9  MPa  for  the  peak  strength*  and  0*3 
MPa  and  0*85  MPa  for  the  reversal  strength*  The  higher  the 
value  of  the  strength*  the  lower  is  the  degree  of 
seoar a  t i on • 

A  similar  approach  was  used  by  Kaiser(  1979)  to  analyse 
the  results  of  several  triaxial  tests  with  confining 
pressures  uo  to  1  MPa*  The  cohesion  intercept  ranges  between 
0*7  and  2*05  MPa  at  oeak  strains*  and  between  0*14  and  0*84 
MPa  at  strains  of  3  %•  These  values  are  affected  by  the  mode 
of  failure*  the  confining  pressure  level  and  the  degree  of 
seoaratior.  of  the  Joints*  The  coal  was  described  as  a 
elastic  strain-weakening  material  with  a  Young's  modulus 


between  850  and  1300  MPa* 


2  *  3* 1  *  2  T i me— dependent  properties 

Kaiser(  1979  )  and  da  Font our a(  1 980  )  have  investigated 
the  1 1  me  — deoen den t  o"operties  of  this  coal*  Kaiser  conducted 
multiple-stage  repeated  relaxation  tests  and  concluded  that 
coal*  as  a  brittle  rock  with  time-dependent  strength,  can  be 
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described  as  a 

s t ra i n— wea ken ing  material 
From  single-stage  and 
Fontoura  determined  an 
pre-failure  range  of  this 


visco-elast ic »  vi sco-plastlc * 

(Kaiser  and  Morgens te rn(  1 979  )  )• 

multiple-stage  creep  tests*  da 
empirical  creep  law  for  the 
c  oa  1  • 


2 «  3 »  2  High-Pressure  Triaxial  Tests 
2  «  3  «  2  «  1,  Introduction 

Tn  order  to  understand  the  behaviour  of  the  coal  at 
high  confining  pressure*  triaxial  compression  tests  were 
conducted  with  confining  pressures  between  3  MPa  and  10  MPa* 
Sample  oreoaration*  testi ng  and  analysis  ere  described  by 
Kaiser(1979)  and  only  a  brief  desclption  will  be  presented 
here  • 

?  *  3  *  2  *  2  Sample  Preparation 

Cylindrical  samples*  4  cm  in  diameter  and  8  cm  in 
length  were  drilled  from  block®  of  coal  in  the  laboratory* 
Water  was  used  as  the  cutting  fluid*  The  axis  of  the  core 
barrel  was  orientated  parallel  to  the  bedding  planes  and  at 
angles  of  30°*  45°  and  60°  respectively  with  the  major  joint 
set.  The  samples  were  cut  with  an  attempt  made  to  keep  the 
major  visible  joints  from  intercepting  the  ends* 


!■  3*  2*  3  Testing  Procedure 


The  samples  were  tested  in  a  Wykeham  Farrance 


hi gh— pressure 


triaxial  cell  designed  for  rock  specimens  up 


i 


14 


to  50  mm  in  diameter  with  a 
of  10  MPa*  The  compress io 
testing  machine  ( Wykeham  Farr 
All  samples  were  consoli 
rate  with  free  drainage*  Th 
The  test  may  be  considere 
strain-rate(  Kaiser(  1979)  ). 

Test  Results 

In  order  to  determine  th 
a  coal  mass  on  the  direct io 
the  principal  stress  directi 
45°  and  60°  with  respect  t 
for  this  study*  A  30°  sample 
£0°  sample  may  be  expecte 
material,  with  a  strength 
intact  material* 

All  the  samples  exhiblte 
the  majority  of  them  failed 
This  is  unfortunate  as  th 
maximum  strength  exhibited 
results  thus  cannot  be  cons 
s t  rength* 

The  results  of  the  elgh 
2*1*  The  strength  envelope 
diagram  which  is  assumed  to 
stress-strain  curves  are  plot 


maximum  a 
n  machin 
ance  )  • 
dated  and 
e  hi ghes  t 
d  as  f  u 


e  depends 
n  of  the 
ons  t  jo i n 
o  the  s  a m 
should  fa 
d  to  fa 
approach i 

d  a  compo 
partially 
e  results 
by  the 
ide red  as 

t  tests  e 
is  pres 
be  linea 
ted  toget 


llowable  cell  pressure 
e  is  a  relatively  soft 

sheared  at  a  constant 
rate  was  0.56  mm  /hr* 
lly  drained  at  this 


nee  of  the  strength  of 
joints  with  respect  to 
t  orientations  of  30°  » 
pie  axis  were  selected 
il  along  a  joint  and  a 
il  through  the  intact 
ng  the  strength  of  the 

site  failure  plane  and 
along  joint  surfaces* 
will  not  indicate  the 
intact  material*  The 
an  upper  limi t  of  the 

re  summarized  in  Table 
ented  on  o  Mohr  circle 
r(  Figure  2*1)»  and  the 
her  in  the  Figure  2*2* 


' 
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Tefcle  2.1  Properties  of  the  Triaxial  Test  Samples 


Test  ID 

A 

B 

C 

D 

E 

F 

G 

I 

kPa 

9308 

9308 

3792 

7240 

10687 

7240 

3792 

3792 

Back  Pr. 
kPa 

345 

345 

345 

345 

345 

345 

345 

345 

03 

kPa 

8964 

8964 

3447 

6895 

10342 

6895 

3447 

3447 

(O3-O 1) 
kPa 

28108 

25564 

16651 

24409 

28490 

22362 

16179 

19761 

(cr  3-a  1  ) 
e=4% 

23700 

24400 

14940 

19792 

21800 

16490 

7590 

14600 

(e=3%; 

E 

MPa 

1221 

1342 

1096 

1153 

1249 

1518 

1164 

1403 

£  0/ 

c  % 

0.75 

0.904 

0.42 

0.80 

0.73 

0.10 

0 

0.624 

p 

f  % 

3.302 

3.273 

2.359 

3.75 

3.401 

1.693 

1.629 

2.280 

* 

Cf  % 

0.25 

0.464 

0.420 

0.85 

0.39 

0.16 

0.24 

0.25 

£  %/hr 

0.28 

0.75 

0.28 

1.92 

0.75 

0.75 

0.28 

0.28 

G 

as  sum 

1.60 

1.60 

1.60 

1.60 

1.60 

1.60 

1.60 

1.60 

g/cA 

1.43 

1.40 

1.34 

1.38 

1.35 

1.40 

1.40 

g/cfl\ 

1.15 

1.13 

1.08 

1.07 

1.08 

1.09 

1.10 

w 

% 

23.86 

23.55 

24.72 

28.06 

25.30 

27.9 

27.42 

S 

r 

0.9824 

0.9104 

0.8096 

0.916 

0.831 

0.961 

0.968 

e 

0.388 

0.414 

0.488 

0.491 

0.487 

0.464 

0.453 

x  % 

(estim. ) 

50 

70 

40 

80 

50 

60 

40  . 

3.  .  _ 

joint 

20  ° 

30  ° 

45° 

45  0 

25° 

60  0 

30  0 

60  0 

ail . 

40  0 

(The  c 

58° 

ngles  8 

45° 

are  take 

40  0 

n  with  r 

30  0 

espect  t 

28° 

0  the  s£ 

60° 

imple  ax] 

26° 

s) 

~~  'J 
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Two  30°  samples  (F  and  I)  failed  along  Joints  on  a 
s  tep- like  surface*  The  remain! ng  sampl e  s  ,  as  Indicated 
before,  failed  partially  through  the  Intact  material  and 
along  joints*  The  Interpretation  Is  complicated  by  several 
other  factors  : 

1 •  The  load  cap  Is  not  oartlcularly  designed  to  allow 

lateral  deformations  and  the  large  deformations  after 
failure  will  thus  be  influenced  by  the  apparatus* 

2*  Subsequent  to  failure  the  upper  corners  of  the  sample 
were  crushed,  which  makes  the  visual  observation  of  the 
failure  zone  difficult* 

3*  As  described  previously,  the  axis  of  the  sample  is 
parallel  to  the  bedding  plane  direction*  In  some  of  the 
samples  the  composite  failure  surface  included  parts  of 
bedding  planes  in  their  step— like  shape* 

Hence  some  of  the  stress-strain  curves  appear  peculiar 
af+er  a  first  load  peak,  and  before  large  strains  occurred 
(see  Figure  2*2),  the  analysis  will  only  consider  the 
behaviour  of  the  sample  up  to  the  first  load  peak*  The 
corresponding  points  are  plotted  as  circles  ir  Figure  2*1 
( "Peek  values"). 

Assuming  a  linear  Coulomb  failure  envelope  and  a  peak 
friction  angle  of  30°,  all  samples  but  sample  I  were  found 
to  have  a  cohesion  intercept  between  1*9  MPa  end  2*4  MPa* 
Sample  I  has  an  unexpec+ly  high  shear  strength  of  3*6  MPa. 
This  may  be  explained  by  the  fact  that  the  failure  surface 
intercepted  the  upper  end  of  the  sample* 
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Figure  2.1  Mohr  Circle  Diagram  for  Triaxial  Compression 
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Figure  2*2  Summary  of  the  Stress-Strain  Loading  Curves 

the  Triaxial  Compression  Tests 
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The  general  tendency  In  the  volume-change  measurements 
is  to  show  compression  prior  to  failure  (0*3%  to  0,7%  volume 
change)  and  dilation  subsequently*  Since  they  are 
incomplete,  they  will  not  be  presented  hereafter  in  detail. 

2_s.fi  Pescrifilion  of  the  Model  Test  Samples 
2*4*1  Introduction 

The  purpose  of  this  section  is  to  describe  the  coal 
blocks  used  in  model  tests  MC-3  and  MC-4.  A  careful 
description,  in  particular  of  the  heterogeneity  of  the 
samples,  is  essential  for  the  evaluation  of  the  experimental 
re  sul t s • 

2.4*7  Sample  MC— 3 

2  *  4  *  2  *  1  Sample  Without  Tunnel 

A  sketch  of  the  sample  is  shown  in  Figure  2*3*  At  first 
sight  the  sample  seems  to  be  particularly  heterogeneous* 
Vegetable  remains,  apparently  from  a  trunk  or  a  root,  were 
exposed  while  trimming  face  1*  The  tnacerals  present  in  this 
lentical  seems  to  be  harder  than  the  average  and  as  a  result 
of  differential  stiffness  a  radial  fracture  pattern 
developed  centred  on  these  remains*  This  is  particularly 
apparent  on  the  bottom  face*  Moreover,  the  main  joint 
pattern  is  perturbed  by  this  hard  zone  from  its  regular 
direction,  and  seems  to  go  round  the  lentical* 

These  features,  as  Interesting  they  may  be  from  the 


geological  point 


of  view , 


are  obviously  not  desired  in  a 
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Figure  2.3  Sketch  of 


Sample  MC-3  before  Preparation 
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sample  for  mode 
possible*  Howeve 
boundaries  of  t 
homogeneous  with 
fractures  appare 
the  whole  sampl 
MC-3* 


The  joints 

incl inat ion 

the 

"hori zontal 

axis"  (  1-3  ) 

1cm 

)  •  Corner  2  —  3 

exh i b i t  s  an 

makes  it  fragile* 


omogeneous  as 
near  the 
is  fa irly 
that  the 
through 
for  test 

varies  between  40°  and  45°  from 
and  are  uniformly  spaced  (about 
existing  fracture  surface  which 


1  test*  which  should  be  as  h 
r,  these  features  are  located 
he  sample  and  the  central  part 
regular  jointing*  It  appears  also 
nt  on  the  faces  are  not  continuous 
e  and  thus  this  block  was  chosen 


2*4*2*?  Core  Dr i 
The  visual 
structure  of  t 
weaknesses  : 
t  •  A  ma j or  pi 

m i d- he  1 ght • 
2*  The  main  f 

towards  face 
3*  Another  fra 

to  the  fir 
through  the 
D\ie  to  these  maj 
five  separate  pi 
the  median  plan 
core*  measured 


lied  in  the  Centre  of  Sample  and  Tunnel  Wall 
observation  of  the  tunnel  wall  confirmed  the 
he  block*  It  showed  three  major  planes  of 

ane  of  separation  crosses  the  sample  at 

racture  no+iced  on  the  top  face  dips  at  45° 

4. 

cture  noticed  on  the  bottom  face  is  parallel 
st  one*  These  fractures  are  not  continuous 
sample • 

or  discontinuities  the  core  was  recovered  in 
eces*  The  core  was  sheared  in  rotation  along 
e  of  separation*  The  moisture  content  of  the 
after  drilling*  was  found  to  be  23*8%  on 


■ 
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averagej  which  indicates  that  the  procedur 
maintain  the  original  mois tu re  of  the  sample  w 
After  drying,  the  oieces  exhibit  a  well 
pattern  with  a  spacing  of  about  1  cm  and 
bedding  planes  with  an  average  spacing  of  about 


2 «  4 * 3  Sample  MC— 4 

2  *  4  * 3  *  1  Sample  without  Tunnel 
This  block  does  not  show 
Jointing  is  regular(l-2  cm)  and  i 
the  block  is  heterogeneous  : 

1*  Some  inclusions  have  been 

trimming  the  block*  They  h 
surface  than  the  other  parts  o 
does  not  appear  to  be 

inclusions*  The  hard  inclusio 
near  comer  3  —  4  on  the  bottom 
2*  The  material  surrounding 

intensively  Jointed* 

3*  Some  discontinuities  may  be  o 

faces  1  and  4  at  45°*  Th 
material  widening  sometimes  to 
inclusions  of  the  same  type  a 
(1  cm  in  diameter)* 

4*  All  over  the  sample  surfac 
material  are  visible*  It  wa 
material  also  filled  the  disco 


any  major 
s  inclined  at 

exposed  by 
ave  a  smooth 
f  the  sample* 
continuous  t 
ns  are  locate 
face* 

the  Indus! 

bserved,  cros 
ey  ©  re  thin  b 
accommodate 
s  the  ones  de 

e  lenses  of 
s  observed  la 
ntinuities  ap 


e  followed  to 
as  efficient* 
defined  joint 
several  major 
4  cm* 


fracture*  The 
4 5°  •  Howe  ve r 


the  saw 

wh  i  1  e 

e r  and 

harder 

The  jointing 

hrough 

t  hes  e 

d  principally 

ons  is 

mo  re 

sing  the  side 

ands  of 

finer 

sever  a 1 

sma  11 

scr i bed 

above 

a  da  rk 

brown 

ter  that  this 
parent  on  the 


■ 


I  *2pj| 
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Plate  2*1  Sample  MC-4  before  Preparation 
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side  faces#  X-rays  analyses  indicated  that  the  material 
was  mainly  composed  of  kaolinlte  (85%),  mon t mo r i 1 loni t e 
(10%)  and  illite  (5%  or  less),  wi th  no  o  rgn  n 1 c 
ma  terial  )• 

2 » 4 , 3 # 2  Core  Drilled  in  the  Centre  of  the  Sample  and  Tunnel 
Wall 

The  tunnel  was  drilled  by  using  a  core  barrel  which 
transmitted  lateral  vibrations  to  the  material#  As  a  result, 
there  was  no  recovery  from  this  drilling  operation#  The  core 
was  reduced  to  small  pieces,  mainly  horizontal  flakes,  with 
an  average  size  of  3cm#  Some  round— edged  pieces  of  coal  were 
recovered#  The  average  moisture  content  of  the  recovered 
materiel  was  as  low  as  18#1%# 

This  low  moisture  content  indicates  that  the  sample 
suffered  from  drying#  Even  if  the  mass  looks  intact  and 
homogeneous,  the  sample  at  this  stage  contains  some 
m i c r o f re c tur in g •  Actually  the  sample  was  left  uncovered  for 
8  hours  during  the  preparation  and  then  heavily  watered 
which  may  have  caused  the  uslacking"  of  the  structure# 
However  if  highly  confined,  the  samole  should  not  be  very 
much  affected  by  this  alteration,  but  the  structural 
alterations  may  reduce  the  strength  of  the  material  around 


th  e  tunnel • 


CHAPTER 


2. 


DESCRIPTION  OF  THE  H^DEL  TEST 


APPARATUS  AND  PROCEDURE 


3  * 1  In  troduc t Ion 

3  « 1 «  1  Preliminary  Remarks 

This  geomechan  teal  model  study  follows  that  of  Kaiser 
(1079)  who  designed  the  test  apparatus^  based  on  equipment 
developed  by  Heuer  and  Hendron  (1967,  1971)#  The  equipment 

and  the  procedure  used  for  this  thesis  are  identical  that 
described  by  Kaiser  ( 1979)  and  a  complete  description  of  the 
system  will  not  be  repeated  here* 


3*1*2  Required  Sped  f  Icatio 
The  most  Important  asp 
summarized  below* 

A  block  of  coal  (6 
under  plane  strain  condi 
has  to  be  provided  in  the  m 
each  corner  of  the  sample 
measurements  and  visual 
head*  A  surface  pressure 
the  lateral  faces  of  th 
maintained  with  minor  fl 
tests  ( >5  days)*  The  la 
Independent  in  order  to  app 
model  block*  A  large  am 
automatic  date  processing 


ns  of  Te 
ects  of 

0cm/60  c 
tions*  A 
Iddle  of 
,  and  fo 
observe  t 
of  maxi 
e  sampl 
uc  tua  tio 
teral  l 
ly  dl f  f e 
ount  of 
is  essen 


st  Apparatus 

the  testing  equipment  are 

m/20  cm)  has  to  be  tested 
ccess  for  instrumentation 
each  lateral  face  and  at 
r  the  tunr.el  deformations 
ions  through  the  loading 
mum  15  MPa  is  required  on 
e  and  will  have  to  be 
ns  during  long  term  creep 
oading  systems  have  to  be 
rent  stress  ratios  to  the 
data  is  expected  and 
tiel  to  control  the  test* 
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This  chapter 
in  accordance 
to  introduce 
the  remainder 
discussion,  see 


describes  brie 
wi th  the  above 
the  reader  to 
of  this  thesis 
Kaiser  (  1979  ) 


fly  the  testing 
specifications 
the  test  appar 
unde  rst  an  dabl  e 


a 


system  design 
The  purpose 
tus  and  to  m& 
For  a  detail 


ed 
i  s 
ke 
ed 


2jl2  E£S£.Hl£±lon  Of  ttje  Xegting  Eauifimeut 

^>2. 1  Compression  Machine 

A  pi c t ur e  of  the  compress! on  test  machine  is  presented 
on  Plato  3.1.  The  sample  with  the  tunnel  in  the  vertical 
position  is  seated  on  a  steel  base  plate.  Plane  strain  is 
maintained  by  a  rigid  loading  head  which  is  controlled  by 
four  rams  (visible  on  the  picture  in  vertical  position).  Two 
rams  act  on  each  lateral  face.  The  four  rams  in  each 
direction  are  connected  to  each  other  and  pressurize  the 
sample  simultaneously.  An  upper  and  a  lower  reaction  head 
are  connected  by  four  rods  and  constitute  the  reaction  frame 
in  the  direction  of  the  tunnel  (called  hereafter 
"longitudinal  direction").  In  the  lateral  directions  (called 
"horizontal"  and  "vertical"  directions)  the  reaction  frames 
are  cantilevered  beams,  supporting  the  rams,  and  restrained 
from  bending  by  two  upper  movable  rods.  A  window  is  provided 
in  the  loading  head  to  make  the  tunnel  accessible  during  the 
test.  The  load  is  distributed  uniformly  to  the  sample  sides 
by  a  set  of  triangular  prisms  in  a  vertical  position. 
Between  the  ram  piston  and  the  prisms  is  inserted  a  load 
cell  (aluminium  hollow  cylinder  with  strain  gauges  glued  on 
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Piute  3# 1  Model  Test  :  Compression  Machine  and  Data 

Acquisition  System 
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the  Inside  wall)  sitting  on  &  spherical,  lubricated  seating 
to  avoid  the  transmission  of  any  bending  moment# 

The  rams  are  controlled  by  three  Independent  hydraulic 
pumps,  one  for  each  direction  of  loading*  A  maximum  surface 
pressure  of  16  MPa  can  be  aoolled  to  the  sample  side  and 
maintained  with  load  fluctuations  of  less  than  one  percent 
for  days  or  weeks*  This  system  allows  maintenance  of  plane 
strain  conditions  and  loading  with  different  ratios  of  the 
aopl i ed  load : 

N=6^/Sy 

6h  b“ing  the  load  in  the  "horizontal"  direction  and  6y  the 
load  in  the  "vertical”  direction* 

2.*.2jl2  Instrumentation 

The  instrumentation  Is  designed  to  monitor  three 
di f ferent  tyoes  of  displacements  t 
t*  Tunnel  deformations  in  four  directions; 

2*  Internal  strains  around  the  tunnel; 

3*  Overall  external  deformations  of  the  block* 

The  deformations  are  measured  by  electronic  transducers 
(  Linear  Vo 1 1 age— D i sp lace men t  Tranducers  or  LVDTs  )• 

The  device  measuring  the  deformations  of  the  tunnel  is 
presented  on  Plate  3*2  ( upoe r  picture)*  A  vertical  stand  is 
fixed  on  the  steel  plate  below  the  sample*  On  each  diameter 
two  LVDTs  measure  the  displacements  of  the  tunnel  wall 
giving  the  closure  of  the  tunnel  with  an  accuracy  of  ±2* 
The  closure  of  the  tunnel  is  defined  as  the  ratio 


' 


. 
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Plate  3-2 


Inst rumen tat  ion 

and 


:  Tunnel  Closure 
Extensometer 


Measuring 


Device 
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between  the  change  In  diameter  and  the  original  diameter  of 
the  tunnel  and  is  exoressed  in  percent* 

The  average  Internal  strains  are  calculated  from  the 
measurement  of  the  relative  displacements  of  two  points  of 
known  original  separation  within  the  rock  mass*  This  is 
achieved  in  practice  by  use  of  the  extensome ters  shown  on 
the  lower  olcture  in  Plate  3*2*  Two  coaxial  rods  are  grouted 
at  their  ends,  in  a  hole  within  the  sample  mass*  The  inner 
rod  Is  glued  at  the  end  of  the  hole  whereas  the  outer  rod  is 
glued  et  a  known  distance  from  the  end  of  the  hole  (5  cm)* 
The  relative  displacement  of  these  points  is  measured  by  an 
LVDT  fixed  to  the  other  end  of  the  coaxial  rods*  The  third 
txibe  shown  on  the  picture  is  a  protective  rod*  The  strains 
are  thus  calculated  with  an  accuracy  of  ±2*10~3*  Two  rings 
with  eight  ex tensome te r s  are  Installed  radially  in  the 
direction  of  the  principal  stress  exes  and  in  the  diagonals 
of  the  sample*  The  first  ring  is  set  up  as  close  as  possible 
to  the  tunnel  wall  (  tio  at  1  cm  from  tunnel  wall),  and  the 
second  ring  measures  the  strain  further  away  from  the  wall 
( tip  at  5  cm  from  the  wall )• 

The  external  deformations  of  the  sample  are  measured  in 
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direction  elves  the  overall  strain  of  the  sample* 


3*3  Sample  Preparation 

Block  of  about  1  m3  are 
described  by  Kaiser  (  1979  ).  They 
shape  of  a  regular  parallelepiped 
are  oerallel  to  one  of  the  diagon 
the  bedding  planes  are  perpend 
faces  are  carefully  manned  and 
parallel  surfaces  and  to  remove 
depressions  and  broken  corners 
grout  and  sanded  until  a  uniform  f 
each  face* 

The  sample  is  then  placed 
removable  sides*  Faces  and  side 
with  a  3  mm  layer  of  Plaster  of  Pa 
of  wax  naoer  isolate  the  plaster  l 
and  from  the  outside.  This  surfac 
to  ensure  a  uniform  transmission  o 

Holes  are  drilled  and  the  ex 
Accurate  records  of  the  sample 
instrument  locations  are  essent 
interpretation*  The  sample  is 
machine  and  the  loading  heads  ar 
the  sample  by  teflon  sheets  t 
sample  faces*  The  LVDTs  are  mounte 
the  load  cells  are  taken*  The  ram 


collected  from  the  mine  as 
are  trimmed  by  saw  to  the 
on*  The  Joints  of  the  coal 
al  plane  of  the  sample  and 
icular  to  the  tunnel*  The 
then  sanded  to  create 
any  loose  material*  Small 
are  filled  in  with  a  stiff 
lat  surface  is  obtained  on 


in  an  aluminium  mould  wit 
s  of  the  sample  are  coate 
ris*  Two  protective  sheet 
ayer  from  the  coel  surfac 
e  preparation  is  necessar 
f  load  to  the  sample* 
tensometers  are  Installed 
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with  the  loading  heads  and  a  zero  reading  of  the  LVDTs  is 
taken*  The  sample  is  now  ready  for  testing* 


3~2-4  Load  Ing  Hi  story 

The  intact  sample  is  loaded  several  times  up  to  10  MPa, 
about  2/3  of  the  capacity  of  the  machine*  This  is  done  in 
order  to  determine  the  elastic  properties  of  the  sample,  to 
compress  the  plaster  layer  and  to  induce  crack  closure* 
Sustained  load  tests  are  also  carried  out  to  study  the 
t 1 m e— depen de nt  compres s lb i 1 i ty  of  the  material  under 
hydrostatic  compression* 

The  tunnel  is  drilled  after  the  loading  head  has  been 
removed*  The  tunnel  wall  is  immediately  covered  with  two 
coatings  of  latex  to  prevent  it  from  drying*  The  sample  is 
reinstalled  end  loaded  by  2  MPa  increments,  after  an 
original  increment  of  5  MPa  up  to  the  upper  limit  of  the 
loading  range  (16  MPa)*  The  load  is  kept  constant  at  each 
load  stage  for  a  period  of  one  or  five  days* 

If  rupture  of  the  tunnel  wall  does  not  occur  the  sample 
is  unloaded  and  left  for  recovery*  In  Chapter  4  the  actual 


loading  histories  for  the  two  tests  are  given  in  detail* 


3  *.5  Description  of  the  Data  Ero£.essjjig  £y£tem 
3*5.1  Data  Processing 

A  data  logger  scans  automatically  the  measurements  of 
44  instruments  (36  LVDT  and  12  load  cells)  at  regular  time 
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intervals*  The  data  are  recorded  on  a  cassette  on  a  computer 
terminal  and  can  be  transfered  directly  to  a  file  in  the 
main  computer  of  the  University*  Datalogger,  transducers 
conditioner,  and  computer  terminal  are  presented  on  the 
lower  picture  in  Plate  3*1* 

Several  computer  programs  allow  one  to  process  the  data 
and  produce  among  other  things  s tress— s tra i n  and  strain— time 
plots*  These  plots  may  be  visualized  on  a  graphic  terminal 
within  less  than  one  hour  after  the  last  reading  has  been 
recorded*  This  rapid  data  analysis  allows  optimal  control  of 
the  test  orocedure* 

3  *  S«  2  Data  Presentation 

Most  data  are  presented  in  the  following  plots: 

1*  Stress-strain  plot  (or  stress— closure ) ; 

2*  Strain— time  plot  (or  closure-time ); 

3.  Log  strain  rate-log  time  (or  closure  rate-time)* 

The  first  two  types  of  plots  have  linear  axes:  there  are 
four  curves  per  diagram,  each  being  identified  by  a  symbol , 
corresponding  to  the  sketch  in  the  upper  right  corner  of  the 
frame*  This  sketch  represents  the  sample; two  diagonal  lines 
Indicate  the  jointing,  the  solid  circle  the  tunnel  and  the 
dashed  circles  the  two  rings  of  internal  ext ensome te rs •  On 
each  figure  the  N-value  calculated  as  an  average  over  the 
data  presented  on  the  figure  is  given.  On  the  strain-time 
plots,  TZERO  is  the  zero  time  of  the  sustained  load  test 
with  respect  to  the  beginning  of  the  test,  and  SIGMAV  is  the 
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Figure  3*1  Typical  Creep  Curves 


TEST  #MC-4 .2  WITH 


A 

CO 

cl 

x 

\ 

V 

LU 

h- 

CL 

cl 

of 

i — 
CO 


A 

CO 

cl 

x 

x 

vO 

ON 

V 

LU 

h- 

CL 

CL 

CL 
I — 
CO 


RELUME  <HRS> 
MC  4,2  RT  11000  KPR 


Flaure  3*2  Typical  Strain  Rate 


Time  Curves 


36 


average  stress  in  the  "vertical"  direction*  Stress  versus 
strain  curves  are  presented  in  Appendix  Bl •  Typical  creep 
curves  (strain  versus  time)  are  presented  on  Figure  3*1  for 
the  closure  of  the  tunnel  and  the  external  deformations* 
Other  curves  are  presented  in  Chapter  4* 

The  strain  rate  versus  time  plots  ere  presented  on  a 
log-log  scale*  The  strain  rate  is  the  slope  of  the  line 
calculated  by  linear  regression  over  three  points  and 
plotted  at  the  time  corresponding  to  the  mid— point*  Only  the 
absolute  value  is  plotted*  If  the  rate  is  positive  a 
straight,  solid  line  is  drawn  between  the  points  and  if  the 
rate  is  negative,  the  points  are  plotted  as  small  solid 
triangles*  After  30  readings  have  been  taken,  a  "trend"  is 
obtained  by  calculating  the  slope  of  the  line  between  three 
points  separated  by  10  readings*  The  results  are  plotted  as 
plus  for  positive  rate  and  crosses  Inclined  at  45  degrees 
for  negative  rate*  Two  typical  examples  of  strain  rate  time 
plots  are  shown  in  Figure  3*2  and  correspond  to  the  creep 
curves  presented  in  Figure  3*1*  Other  curves  showing 
negative  strain  rates  are  presented  in  Chapter  4* 


CHAPTER  4 


DESCRIPTION  OF  XfiE  TESTS  AND  EXPERIMENTAL  RESULTS 


In  troduc  t  Ion 

Two  different  specimens  were  tested  in 
and  MC  — 4  )  •  Each  sample  was  subjected 
sequential  loadings  (e.g.  Tests  MC— 3.0 
actual  loading  histories  are  described  in 
4.3  and  are  based  on  the  procedure  discusse 
They  are  also  presented  in  graphs  in  Append 
Too  many  data  have  been  acquired  to  b 
here*  All  data  has  been  summarized  in  a 
(1979)  of  the  Geotechnical  Section  of 
Civil  Engineering,  University  of  Al 
st r e ss— st r a i n  curves  are  reported  in  full  i 
are  given  in  Appendix  Bl.  Some  representat 
and  strain  rate  —  time  curves  are  shown  in 
Although  similar,  the  two  tests  d 
aspects  : 
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The  Interprets tl on  of  the  experimental  data  may  be 
subdivided  Into  three  parts;  elastic  properties,  strength 
properties  and  time  —  dependent  ( vlsco  —  elastic  and  vlsco 
-plastic)  properties. 

While  the  loading  of  a  classical  triaxial  sample  can 
provi de  independently  an  elastic  constant  (e.g.  Young* s 
modulus ),  the  loading  of  the  model  test  cannot  permit  the 
distinction  between  two  elastic  parameters  which  influence 
simultaneously  the  measurements  (e.g.  Young*  s  modulus  and 
Poisson’s  ratio).  An  attempt  was  made  to  determine  these 
elas+ic  parameters  independently  (Section  5.2).  It  led, 
however,  to  scattered  results.  Because  of  this  complexity 
the  elastic  deformations  will  not  be  analysed  in  detail.  The 
assumotion  of  a  Poisson’s  ratio  of  0.2  leads  to  values  for 
Young’s  modulus  between  1000  and  1700  MPa.  These  values  are 
si lgh+ly  higher  than  the  results  of  triaxial  tests. 

The  strength  properties  will  be  analysed  in  Chapter  5. 

The  visco  -  elastic  properties  of  the  material  in  the 
pr  e  —  f  ai  lur  <=»  range  of  the  loading  have  been  analysed  in 
detail  by  Kaisp"( 1979)  based  on  the  results  of  previous 
tests.  The  time  -  dependent  behaviour  of  samples  MC-3  and 
MC- 4  are  in  agreement  with  his  conclusions  which  are 
summarized  in  Section  4.4.  The  analysis  of  the  time  - 
dependent  plastic  deformations  is  beyond  the  scope  of  this 
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4  *_2  Samelfe  i  Loading  History 

4 « 2. 1  Introduction 

Sample  MC— 3  was  collected  and  trimmed  during  winter 
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2 

Loading 

H i story 

Intact  sample 

*  Test  MC—  3*0  *  Aug. 30th  *  N=0.75 

—  Loading  :  8500  kPa  t  =  2h. 

—  Unloading  t  1000  kPa  t=4h. 

—  Reloading  t  4500  kPa  t=7h. 

Test  MC-3.01  *  Sept.  1st  *  N=1.01 

—  Loading  :  6000  kPa  t  =  2h. 

—  Unloading  t  0  kPa  t=2.5h. 

*  Test  MC-3.02  *  Sept.  4th  *  N=1.09 

—  Loading  :  9500  kPa  t=3.5h. 

—  Unloading  t  0  kPa  t=4h. 
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*  Test  MC-3.03  *  Sept.  6th  *  N=1.01 

—  Loading  :  8500  kPa  t=2.5h. 

—  Creep  :  8500  kPa  t=20h. 

—  Unloading  :  0  kPa  t=20h. 

—  Recovery  :  t=45h. 

The  tunnel  was  drilled  on  September  8th  1978.  The 
drilling  operation  took  about  10  minutes  •  The  core  was 
recovered  in  five  major  pieces  end  is  described  in  Section 

2. 4.2.2. 

Sample  with  Tunnel 

*  Te« t  MC-3.1  *  Sept.  28th  * 

—  Loading  :  4781  kPa  t=2h. 

—  Creep  :  N=l.06  t=120h. 

—  Leading  :  6917  kPa  t=121h. 

—  Creep  N=1.02  t=238h. 

—  Loading  :  9000  kPa  t=240h. 

—  Creep  N=1.00  t=350h. 

—  Loading  :  11117  kPa  t=350h. 

—  Creep  N=0.99  t=609h.  (double  duration) 

—  Loading  :  13336  kPa  t=610h. 

—  Creep  N=0.98  t=719h. 

—  loading  Z  15000  kPa  tr=715h. 

—  Unloading  t=717h 

—  Recovery  t=1067hr 

*  Test  MC-3.2  *  December  9th 

—  Loading  Z  14902  kPa  t=3.5h. 


Creep  N=0.97  +=6.5h. 
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—  Uncontrolled  unloading  after  the  rupture  of  a  corner  of 
the  samole 
4  *.2  *.3  Comments 

The  complete  stress— strain  curves  are  presented  on 
pages  Bl*l  to  Bl*12  in  Appendix  Bl*  This  test  was  disturbed 
by  several  accidents*  It  happened  several  times  during  the 
tests  that  a  hydraulic  ram  leaked  internally  and  was 

temporarily  unable  to  maintain  a  constant  pressure  on  the 
sample*  This  caused  a  break  in  the  symmetry  of  the  applied 
pres  sures  and  the  sample  rotated  slightly*  In  the  worst 
cases  the  sample  had  to  be  unloaded  and  the  ram  repaired* 
This  hapoened  at  the  end  of  test  MC— 3*1,  preventing  any 
creep  test  during  the  last  load  increment*  This  also 

accounts  for  the  Important  time  lapses  during  the  loading 
program*  New  rams  with  higher  load  capacity  have  now  been 
built  to  prevent  these  delays* 

The  sample  felled  by  shearing  of  corner  2-3  along  a 
plane  parallel  to  the  joints*  After  being  loaded,  the  sample 
rotated  slightly  during  the  creep  test  and  the  rams,  no 

longer  perpendicular  to  the  sample  faces,  created  stress 

conditions  which  sheared  the  sample  along  this  plane  of 


w 


eakness •  The  readings  of  the  strain  gauges  do  not  show 


evidence  of  the  shearing  process  at  the  corner*  This 
suggests  that  the  presence  of  a  major  plane  of  weakness 
parallel  to  the  tunnel  axis  is  not  critical  for  the 
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Sample  M£~4  1  Loading  History 


4«.3  *1 


Sample  MC-4  was  collected  and  trimmed  durln 
1976—77  and  instrumented  in  January  1979*  It  was  tes 
a  2  month  period  between  February  3rd  and  March  23 
The  intact  sample  was  subjected  to  7  loadings  and  th 
with  tunnel  to  two  multiple-stage  creep  tests# 

Because  the  previous  samples  did  not  s 
significant  ruDture  of  the  tunnel  wall ,  it  was  de 
increase  the  tunnel  diameter  by  2*3  cm  •  As  the  beha 
the  sample  in  the  prefailure  range  has  been  investig 
because  most  of  the  creep  occurs  in  the  first  day*  t 
test  duration  was  reduced  to  1  day# 


g  winter 
ted  over 
th  1979# 
e  sample 

how  any 
cided  to 
vi our  o  f 
ated  and 
he  creep 


4#  3.2  Loading  History 
Intact  sample 

*  Test  MC-4.0  *  Feb.  3rd  *  N=  1 . 04 

Loading  :  10000  kPa  t=1.2hr  (test  MC-4.01) 

—  Creep  t=2#2h# 

—  Unloading  :  0  kPa  t=3#5h# 

—  Time  laps  t=5h# 

—  Reloading  I  10000  kPa  t  =  6  •  2h •  (test  MC— 4# 02) 

—  Unloading  :  0  kPa  t=7#4h# 

*  Test  MC-4. 03  *  Feb.  4th 

—  Loading  t  4985  kPa  t=0.9h. 

—  Creep  N=1.08  t=24h. 


Load i n  g 


7004  kPa  t=25h. 


. 
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—  Creep  N=1.07  t=48h. 

—  Loading  :  9294  kPa  t=49h. 

—  Creep  N=1.00  t=71h. 

Loading  :  11356  kPa  t=73h. 

—  Creep  N=1.00  t=95h. 

—  Unloading  :  0  kPa  t=99hr. 

—  Recovery  t=119h. 

*  Test  MC  — 4  •  04  *  Feb.  10th 

—  Loading  :  9168  kPa  t  =  1.8h. 

—  Creep  N=1.00  t=24h. 

—  Unloading  :  0  kPa  t=26h. 

—  Recovery  t=42h. 

*  Te«  t  MC-4.11  *  Feb  22nd 

—  Loading  :  5008  kPa  t=0.8h. 

—  Creep  N+1.07  t=5.5h. 

—  Unloading  :  0  kPa  t=6.5h. 

—  Recovery  t=22h. 

*  Test  MC  4.13  *  Feb  23rd 

—  Loading  :  4995  kPa  t=1.2h. 

—  Creep  N=1.02  t=24h. 

—  Loading  t  7166  kPa  t=25h. 

—  Creep  N=1.02  t=49h. 

—  Loading  :  9419  kPa  t=49.3h. 

—  Creep  N= 1 • 00  t=73h. 

—  Loading  J  1535  kPa  t=74h. 
Creep  N=1.00  t=97.5h. 

—  Unloading  t  0  kPa  t=105h. 
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—  Recovery  t=182h* 

The  tunnel  was  drilled  on  March  5th  In  about  10 
minutes*  The  core  was  recovered  in  small  pieces  and  is 
described  in  Section  2*4  • 

Sample  with  tunnel 

*  Test  MC-4.2  *  March  8th 

—  Loading  :  4499  kPa  t  =  0*7h* 

Creep  N=l*ll  t=24h. 

—  Loading  :  6469  kPa  t=25h* 

—  C^eep  N=t*13  t=49h* 

—  Loading  Z  8525  kPa  t=50h* 

—  Creoo  N=1*10  t=73h* 

Loading  :  11109  kPa  t=74h. 

—  Creep  N=l*02  t=97*3h* 

—  Loading  Z  13354  kPa  t=98h* 

—  Creep  N=0*98  t=123h*  (first  rupture  of  the  tunnel) 

—  Loading  :  15967  kPa  t=128h* 

—  Creep  N=0.95  t=15?h* 

—  Unloading  :  0  kPa  t=153h* 

—  Recovery  t  =  2  3  8h • 

*  Test  MC  —  4*3  *  March  21th 

—  Loading  :  7129  kPa  t=l*lh* 

—  Creep  N  =  0  •  9  9  t=24h* 

—  Loading  Z  9361  kPa  t=25h* 

—  Creep  N=0*98  t=48h* 


Load! ng 


11753  kPa  t=49h • 
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Creep  N=0.95 
Loading  :  141 

Creep  N=0.94 
Loading  :  150 

Creep  N=0.95 
Loading  :  152 
Creep  N=0.95 
Loading  :  166 

Creep  N=  0.  q0 
Unloading  :  0 
Recovery  +=42 


t=73h. 

09  kPa  t=74h. 
t=  9Sh • 

15  kPa  t=99h. 
t=l  15h. 

16  kPa  t  =  1 1 5h 
t= 1 42h. 

28  kPa  t  =  l 43  h 
t=l 69h. 

kP a  t  = 1 71 h • 
lh. 


4  Comments 

The  stress  —  strain  curves  are  presented  in  Appendix  31 
(page*?  Bl.13  to  Bl.30).  This  test  is  characterized  by  the 
ruoture  of  the  tunnel  wall  during  the  5tb  increment  in  test 
MC  — 4.2.  The  rupture  is  described  In  detail  in  Chapter  6. 

During  the  first  sustained  load  test  in  test  MC— 4.01, 
the  Intact  sample  exhibited  high  time  —  dependent 
deformations  (0.12%  in  30  minutes).  At  the  same  time  water 
was  expelled  from  the  sample  as  drops  forming  along  the 
coaxial  strain  gauges.  This  water  had  previously  been 
Introduced  in+o  the  sample  cracks  by  a  heavy  watering  of  the 
sample  during  the  preparation.  Consequently,  the  inner  rods 
of  the  extensomet e rs  rusted.  This  explains  the  peculiar 
results  of  tests  MC-4.03  and  MC-4.04  during  which  an 


. 
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unsuccessful  attempt 
The  gauges  wer 
results  Improve  the 
changed  and  it  was  n 
3-4t  second  ring  o 
installed  to  monit 
applying  the  Dress 
order  to  compare 
displacements  of  t 
Chapter  3)  to  the 
measurement  is  given 
and  the  correspondin 
reference  sketch* 

The  data  for 
Appendix  Bl:  the 

only  and  the  seal 
deformation?;  the 
whole  test  where 
scale • 


was  made  to  clean  th 


e  changed  at  the  end 


change*  Only  one  inst 


ot  monitored  further 


f  ex tensometers ) •  The 


or  the  displacement 


ure  on  the  sample  fac 


the  measurements 


he  sample  given  by  th 


displacements  of 


in  the  figures  prese 


g  symbol  is  in  the  ri 


e  ext ensome t e rs • 
of  test  MC-4.04  and 
rument  could  not  be 
(lower  right  corner 
available  LVDT  was 
s  of  the  triangles 
e*  This  was  done  in 
of  the  external 
e  usliel  method  (  see 
the  triangles*  This 
nted  in  Appendix  Bl 
gh t  position  on  the 


test  MC— 4*2  are  presented  twice  in 
first  series  shows  the  data  for  loading 
e  allows  the  analysis  of  the  small 
second  series  presents  the  data  for  the 
the  large  deformations  require  a  larger 


A  safety  system  protects  the  testing  system  from 
excessive  displacements  of  the  rams*  Micro^swi tches  are 
Installed  close  to  the  loading  heads  of  the  rams  during 
creeo  tests*  These  switches  control  electro  valves  which 
stop  the  air  pressure  supply  if  large  displacements  of  the 
ram  occur*  In  test  MC-4 • 3  this  line  was  accidentally 
switched  off  for  30  minutes  in  the  first  creep  test  at  1500 
kPa •  The  load  in  the  rams  dropped  anisotropically  which 


- 
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exolalns  the  particular  aspect  of  the  stress 
for  this  test*  The  sample  was  reloaded  for 
test  at  an  equivalent  stress  level* 


—  strain  curve 
a  second  creep 


4*4  Tl me- Dependent 

4*4*1,  Reliability  of  the  Time  —  Dependent  Data 

The  testing  equipment  and  procedure  are  satisfactory  to 
study  the  time  —  dependent  deformations  around  the  tunnel* 
Responses  of  an  Instrument  in  a  sustained  load  test  are  very 
accurate  in  particular  In  a  recovery  test*  Typical  curves 
are  given  on  Figures  4*1  and  4*3*  Responses  of  the  closure 
measuring  device  have  already  been  shown  in  Chapter  3*  The 
procedure  to  calculate  the  strain  rates  seems  successful  ! 
most  of  the  da+a  show  a  lineer  relationship  between  log 
strain  rate  and  log  time*  The  slope  of  this  line  is  negative 
and  is  more  or  less  constant  throughout  a  test*  This 
corresoonds  to  the  re  sponse  of  a  rheological  model  composed 
of  a  series  of  Kelvin  models  with  a  particular  retardation 
soectrum*  This  model  has  been  Droposed  to  explain  the  time  — 
deoendent  behaviour  of  coal  (Terry(  1956  )  )  Typical  strain 
rate  —  time  curves  are  shown  In  Figure  4*2  and  correspond  to 
some  of  the  strain  -  time  curves  presented  on  Figure  4*1 

4i4i2  Interpretation  of  the  time  —  dependent  data 

Similar  time  -  dependent  data  in  the  pre  -  failure 
range  have  been  analysed  by  Kaiser(  1979)*  He  explained  the 
time  -  dependent  process  by  a  non  -  linear  bulk  modulus  at 
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Figure  4*1  Creep  Curves  of  Internal  Strain  Gauges 


TEST  #MC-4 . 2  WITH  TUNNEL  1979 
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REL  •  T I  ME  <HRS> 


TEST  #MC-3.1  WITH  TUNNEL  1978 


RECOVERY 


Figure  4.3  Recovery  Curves  end  a  Corresponding  Strain  Rate  - 

Time  Curve 
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low  stress  level,  superposition  of  hydrostatic  and 
deviatorlc  creep  a t  Intermediate  stress  levels  and  by  an 
Isotropic  or  anisotropic  time  -  dependent  stress 
redistribution  at  higher  stress  levels.  He  also  mentioned 
the  Importance  of  the  anisotropic  character  of  the  time  - 
deoendent  process. 

The  time  —  dependent  data  near  failure  are  described  in 
Chapter  6  but  a  detailed  interpretation  is  beyond  the  scope 
of  this  thesis. 


ii.5  Sutamgu^r  and  Conclusions 

Two  model  tests  samples  have  been  tested  during  the 
eight  months  of  the  testing  program  for  this  thesis.  For  the 
first  time  in  this  research  program  a  test  was  carried  out 
beyond  the  initial  rupture  of  the  tunnel  wall.  Because  the 
data  In  +he  pre— failure  range  are  in  agreement  with  the 
conclusions  of  Kaiser(  1979  )  based  on  similar  data,  the  test 
interpretation  has  been  oriented  towards  the  analysis  of  the 
rupture  process.  The  different  ultimate  behaviour  of  the  two 
tests  deserves  a  comparative  study  of  the  failure  process. 


This  is  given  in  Chapters  5  and  6. 
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INTERPRETATION  OE  STRESS-STR AIN  &LXL  NEAR  FAILURE 


E  .  1  Intro  due  t  ion 

5 * 1 * 1  Steps  In  the  Data  Interpretation 

Visual  observations  on  both  tests  indicate  that  plastic 
deformations  occur  before  the  completion  of  the  loading 
program.  Since  this  initiates  on  the  tunnel  wall,  therefore 
only  the  tunnel  closure  measurements  will  be  analysed  in 
detail*  The  Internal  extensometers  measurements  are 
perturbed  by  local  effects*  Therefore  their  interpretation 
is  more  difficult  and  it  will  rot  be  undertaken  here*  The 
behaviour  of  the  sample  near  failure  will  be  studied  by 
comparison  of  the  experimental  data  with  the  response  of  a 


theoretical  elasto— plas t i c  model* 


5.*.JL.*.2  Existing  Methods 

Numerous  methods  employing  e la s t o— pi a s t i c  analyses  are 
available  to  evaluate  the  stability  of  an  opening. 
Dae men ( 1975)  reports  the  most  important  methods  and  gives 
his  own  contribution,  the  use  of  particular  strain-weakening 
models.  Egger(  1973)  was  one  of  the  first  to  investigate  the 
effect  of  the  post-failure  behaviour  of  the  material* 
De  scoeudres(  1 977  )  gives  a  detailed  derivation  of  the  solution 
which  employs  an  elastic  perfectly  plastic  and  elastic 
perfectly  brittle  material  model*  Florence  and  Schweri  197S  ) 
point  out  the  importance  of  the  intermediate  principal 
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stress*  Borsetto  and  R 1 bacchi(  1 979  )  studied  the  Influence  of 
a  strain-weakening  behaviour  after  strength  failure*  A  ring 
of  "weakening"  material  between  the  remaining  elastic  and 
the  "res idual"  material  provides  an  apparent  support 
pressure • 

Most  of  these  authors  analyse  the  unloading  of  a 
prestressed  rock  mass  due  to  tunnel  advance*  They  also  take 
Into  account  an  internal  support  pressure  provided  either  by 
the  nearby  tunnel  face  or  by  an  artificial  support*  The  aim 
of  these  theoretical  methods  is  to  evaluate  criteria  for  the 
design  of  underground  cavities* 

Model  tests  have  the  same  stress  and  strain 
distributions  as  encountered  in  an  actual  tunnel,  but  do  not 
have  the  same  displacement  field*  At  an  Infinite  distance 
from  the  tunnel  the  displacement  in  the  model  test  has  a 
finite  value*  The  test  behaviour  results  from  the 
simultaneous  loading  of  a  finite  intact  plate  and  unloading 
of  an  Infinite  medium  due  to  a  cavity*  As  a  consequence, 
particular  solutions  for  the  model  test  have  to  he  adapted* 

5.1*3  Outline  of  the  Analysis 

When  testing  a  nonhomogeneous  material  with  time  - 
dependent  properties,  it  appears  reasonable  at  the  outset  to 
try  and  explain  the  observations  by  simple  methods, 
preferably  with  closed-form  solutions.  This  allows 
identification  of  the  processes,  the  dominant  parameters  and 
provides  a  basis  for  more  advanced  calculations* 
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Closed-form  solutions  for  the  model  test  are  develooed 
In  the  second  s  ec  t i on  of  this  chapter*  These  solutions 
facilitate  the  analysis  of  the  influence  of  different 
e l as t o— p las 1 1  c  models  on  the  deformations  of  an  opening  near 
failure,  as  a  function  of  material  properties*  They  require 
very  restrictive  assumptions,  such  as  homogeneous  and 
t ! me— I ndependen t  characteristics  of  the  material  and 
axlsymmetric  stress  conditions*  Experimental  results  permit 
the  validity  of  the  solutions  and  the  importance  of  the 
assumotlons  to  be  evaluated*  This  analysis  is  presented  in 
the  third  section*  A  first  section  presents  a  method  for 
evaluation  of  the  elastic  parameters  of  the  material  by  use 
of  the  experimental  data* 


F «  2  Determination  of  Elas tl c  Parame  ter 9 
5*2* 1  Intro duct  ion 


It  will 

be 

shown  .1  n 

the  next  sections  tha 

t 

one 

tatlon  in 

the 

analys i s 

of  the  experimental  data 

i  s 

due 

n  on— l 1  he  ar 1 ty 

of  the  e la 

stic  response  of  the  syst 

em  • 

It 

t  here  fore 

of 

1 n  t ere  s  t 

to  know  the  variation 

of 

the 

elastic  constants  with  stress  level* 

The  slope  of  the  response  of  any  instruments  measuring 
the  deformations  in  the  model  test  is  always  a  function  of 
two  elastic  constants,  and  one  curve  is  not  enough  to  know 
the  variation  of  each  parameter  independently*  It  is 
possible  to  get  this  variation  by  comparing  the  responses  of 
two  different  instruments.  Two  methods  have  been  derived  and 
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are  explained  In  detail  in  Appendix  A4.  Both  are  based  on 
elasticity  and  are  only  valid  before  the  yield  point. 

5.2.2  First  Method 

The  first  method  consists  of  comparing  the  responses  of 
an  internal  strain  gauge  before  and  after  excavation  of  the 
tunnel.  It  assumes  that  the  elastic  constants  are  not 
affected  by  the  loading  history  and  that  the  rate  of  loading 
is  the  same  in  the  two  different  tes+s.  For  each  model  test 
the  method  has  been  applied  to  the  deformations  of  the  16 
internal  gauges.  This  study  is  valid  because  the  rates  of 
loading  are  similar,  and  because  the  effect  of  the  loading 
history  was  not  important  e ny  more,  as  the  sample  had  been 
loaded  several  times  before  the  studied  tests.  The  results 
are  presented  in  Apoendix  A4.  They  are  extremely  scattered. 
The  values  of  Poisson's  ratio  are  between  0  and  0.2  and  of 
Youns*  s  modulus  between  1300  and  2300  MPa. 

It  can  be  seen  from  the  expressions  given  in  Appendix 
A 4  that  the  accuracy  of  the  values  of  Poisson*®  ratio  is  a 
function  of  the  distance  of  the  gauge  from  the  tunnel  wall. 
In  other  words  the  further  the  gauge  is  from  the  tunnel,  the 
less  is  the  effect  of  the  presence  of  an  opening.  Thus  the 
results  are  totally  erratic  for  the  furthest  row  of  gauges 
from  the  tunnel  wall.  The  results  for  the  other  instruments 
are  affected  by  the  Inaccuracy  in  determining  the  exact 
position  of  the  gauge  and  by  the  heterogeneity  of  the 
sa  mnl e . 


. 
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5*2*3  Second  Method 

The  second  method  deals  with  one  test  only,  in 
comparing  the  response  of  a  tunnel  diameter  and  of  +he 
external  deformation  of  the  block  in  the  same  direction*  The 
limitations  encountered  in  the  first  method  by  using  two 
different  tests  with  different  loading  histories  are  thus 
overcome  but  only  two  values  (Instead  of  16)  are  now 
obtained  in  each  model  test  (one  for  each  of  the  vertical 
and  horizontal  directions)*  The  results  are  presented  in 
Appendix  A4*  They  show  a  tendency  for  the  Poisson1 ratio  to 
increase  from  a  value  near  0  to  a  value  near  0*2*  The 
correspond ing  values  for  Young’s  modulus  range  between  1300 
and  2300  MPa* 

The  value  of  Poisson’s  ratio  appears  to  be  lower  than 
+  he  values  published  for  coal  (  Labasse(  1.949  )  ,  Ko  and 
Gers t le(  1 976  )  ),  which  are  around  0*4*  However,  the  high 
density  of  jointing,  the  presence  of  bedding  planes  and  the 
state  of  frecturation  make  that  the  overall  value  of 
Poisson’s  ratio  for  the  rock  mass  to  be  lower  than  this 
published  value.  While  0*4  is  a  value  of  Poisson’s  ratio  for 
coal  as  a  mineral,  0*2  or  less  is  the  value  to  be  input  in 
any  theoretical  analysis  of  the  behaviour  of  the  coal  mass* 

5*2*4  Conclusions 

The  methods  used  to  calculate  the  elastic  constants  do 
not  aive  accurate  results,  but  Indicate  a  trend*  They  could 
be  used  to  characterize  tbe  variation  of  stiffness  in  the 
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samole  on  a  qualitative  bas 
an  orthotropic  structur 
considered  In  further  work 
study  of  an  orthotropic  coa 


is. 

Similar  methods, 

but 

assuming 

e 

for  the  materi 
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must  be 

(  see 
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)  for  the 

1  ). 

.3  Elasto -Plastic  Models 
5 . 3. 1  Assumptions 
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f  i  eld  stress* 

The  material  properties  are  time-independent* 

5*3*?  Elas tic— Perf ec t ly  Plastic  Model 

Only  the  logic  of  the  solution  Is  presented  here*  The 
comolete  derivation  is  given  In  appendix  A?* 

5 . 3 *  ?-  *  1  Elastic  Behaviour 

The  elastic  analysis  of  a  cavity  in  an  infinite  medium 
was  originally  attributed  to  Kirsch*  The  derivation  of  the 
stress  distribution  by  use  of  Airy's  stress  functions  is 
found  in  numerous  textbooks  (  i*e*  Obert  and  DuvalK  1S67  )  )• 

Hooke's  law  will  be  used  to  obtain  the  strain  distribution* 
The  integration  of  the  strain-displacement  relationships, 
with  +  he  appropriate  boundary  conditions,  will  give  the 
displacement  field*  The  elastic  closure  of  the  tunnel, 
defined  as  the  ratio  between  the  change  in  diameter  and  the 
original  diameter,  is  equal  to 

u/a=6g(  1-nu2  )/E 

where  u  Is  the  displacement  at  the  tunnel  wall ,  a  the  radius 
of  the  opening,  5g  the  external  stress,  nu  is  the  Poisson's 
ratio  and  E  the  appropriate  modulus  of  elasticity* 

5 . 3  *  2  «  2  Plastic  Behaviour 

If  the  tangential  stress  at  the  tunnel  wall  exceeds  the 
unconflned  compressive  strength  of  the  material,  failure 
occurs  uniformly  around  the  tunnel*  An  annular  zone  of 
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yielding  material  develops  and  propagates  into  the  rock  mass 
as  the  load  is  increased*  For  a  given  value  of  external 
stress,  the  stress,  strain,  and  displacement  fields  and  the 
radius  R  of  the  plastic  zone  are  determined  in  the  following 
manner* 


The 

combination  of 

the 

equilibrium  equation 

and 

t  he 

fa i lure 

criterion  gives 

the 

stress  distribution 

in 

the 

plastic  zone*  By  assuming  continuity  of  the  radial  stresses, 
the  radius  of  the  plastic  zone  can  be  determined*  If  the 
boundary  between  the  plastic  and  elastic  zones  is  known,  the 
radial  stress  at  this  level  is  given  and  may  be  considered 
as  a  support  oressure  for  the  remaining  elastic  zone* 
Therefore  the  stress,  strain,  and  displacement  fields  can  be 
determined  in  the  elastic  zone* 

Strains  in  the  plastic  zone  are  assumed  to  be  the  sum 
of  elastic  and  plastic  components*  It  is  further  assumed 
that  Hooke  *  s  law  aoolies  for  the  elastic  part  of  the  strain 
in  the  plastic  zone,  with  the  same  elastic  constants  as  in 
the  elastic  zone*  As  the  stresses  in  the  plastic  zone  are 
known,  the  elastic  part  of  the  strain  can  then  be  determined 
by  applying  Hooke's  law* 

In  order  to  obtain  the  plastic  part  of  the  strain,  a 
flow  rule  must  be  postulated.  The  following  relationship 
between  the  radial  and  tangential  plastic  strains  was  chosen 
(  Egger  ,  1973  )• 

6r+a6^=  0 

6 

The  parameter  a  is  a  measure  of  the  dilation  associated  with 
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the  plastic  flow*  If  cr  is  equal  to  one*  the  plastic 
volumetric  strain  is  equal  to  zero  since  the  longitudinal 
strain  is  considered  to  be  zero  during  plane  strain  loading* 
It  can  be  shown  that  the  strain  increments  are  then  normal 
to  the  Von  Mises  criterion  failure  surface  ( SalenjonC  1969  ) 


).  The 

associated  flow  rule  together  with  the  Coulomb 

fa i lure 

criterion  will  lead  to  a  similar  relationship 

between 

principal  plastic  strains,  with  the  boundary 

conditions  stated.  This  relationship  is  equivalent  to  the 
above  expression  where  a  is  replaced  by  m»  coefficient  of 


f a i lure 

earth  pressure  and  also  coefficient  of  the  Coulomb 

f  a  i  lur  e 

criterion.  The  parameter  a  will  then  be  between  1 

and  m • 


The 

combination  of  the  compatibility  equation  and  this 

relationship  will  give  the  plastic  strain  and  the 
displacement  field.  The  closure  of  the  tunnel  under  these 
conditions  is  as  follows  : 


/a  £  \ 

—  -  c  (i+m+^  rfR)m+a-ii) 

a  E  u  m+a  LV  IJj 

The 

normalized  closure  has  been  introduced  as  the  ratio 

be  tween 

the  closure  given  by  the  above  expression  and  the 

closure  of  a  tunnel  in  a  linear  elastic  material.  This  ratio 
is  one  up  to  the  yield  point  and  increases  for  higher 
stresses  as  a  function  of  the  strength  parameters  C  and  0» 
and  of  the  dilation  characteristics  of  the  material.  It  is 
Independent  of  the  elastic  properties  of  the  material.  If 
the  strength  parameters  are  time  independent*  the  normalized 
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closure  is  independent  of  the  rate  of  load! 
however,  only  if  the  loading  is  hydrostatic 
redistribution  occurs  within  the  sample* 


ng. 

1 1 

i  s 

val id, 

and 

i  f 

no 

st  r ess 

5  »  3*2*3  Volume  Change 


The  dilation  parameter  a  has 


sensitivity  of  the  closure  to  this 


The  amount  of  volume  change  induce 
value  of  this  parameter  and  the  var 


for  different  values  of  at  will  be 


study  was  not  found  in  the  llterat 


authors  are  concerned  with  real  tunn 


does  not  present  any  particular  int 


for  the  model  test  where  the  volume  o 


The  volumetric  strain  is  the  s 


principal  strains*  It  is  determined  a 


the  strain  field  is  known*  The  vo 


defined  as  the  change  in  volume  of  a 


the  point*  The  volume  change  of  a  bo 


the  entire  volume  of  the  body  of  this 
can  be  shown  that  if  the  volume 
throughout  a  body,  the  corresoondi 
body  is  equal  to  this  constant  value* 


The  volumetric  strain  is  constan 


and  varies  in  the  plastic  region  by 
In  reality,  this  constant  value  is 


occurs  in  the  plastic  zone*  Fowe 


been  Introduced  and  the 
value  will  be  studied* 
d  by  using  a  particular 
iation  in  volume  change 
given*  Such  a  detailed 
ure  because  most  of  the 
els  where  volume  change 
erest*  This  is  not  true 
f  material  is  finite, 
um  of  the  values  of  the 
t  each  point  as  soon  as 
lumetrlc  strain  is  also 
unit  volume  surrounding 
dy  is  the  Integral  over 
unit  volume  change.  It 
trie  strain  is  constant 
ng  volume  charge  of  the 

t  in  the  elastic  region 
becoming  more  dllatant* 
zero  and  only  dilation 
ver  in  the  model,  this 
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constant  value  Is  compression  due  to  the  displacement  of  the 
external  boundaries  and  this  compression  is  superposed  on 
the  dilation  in  the  plastic  zone* 

The  volumetric  strain  in  the  plastic  zone  is  the  sum  of 
an  elastic  and  a  plastic  component*  If  the  dilation 
parameter  increases,  the  elastic  part  remains  constant  in 
the  olastlc  zone  as  it  is  a  function  only  of  the  stress 
distribution.  The  total  volume  change  is  only  equal  to  the 
elastic  volume  change  if  the  dilation  parameter  is  1*  If  the 
dilation  parameter  Increases,  the  plastic  part  of  the 
volumetric  strain  changes  and  becomes  more  dilatant* 

The  integration  of  the  varying  values  of  the  total 
volumetric  strain,  through  the  plastic  zone,  gives  a 
constant  average  value  which  will  be  called  hereafter  the 
’•average  volume  change”  (  LadanyK  1  974  )  ) 

Pieure  5*1  shows  typical  closure  stress  curves  and  the 
corresoondlng  average  volume  change  versus  closure  curves* 
The  strain  due  to  the  compression  of  the  plate  only  is  given 
as  a  reference.  The  behaviour  of  the  model  test  Is  shown  for 
a  ourely  elastic  material  and  for  an  elastic-perfectly 
plastic  material  with  different  values  of  the  dilation 
coefficient*  The  elastic  and  strength  parameters  are  typical 
for  the  behaviour  of  the  actual  model*  The  friction  angle  is 
30  degrees;  the  corresponding  value  of  the  passive  earth 
pressure  coefficient  m  is  3,  which  is  the  highest  value  for 
the  dilation  parameter  on  the  figure  (associated  flow  rule). 
The  values  1  and  2  correspond  respectively  to  no  volume 
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Figure  5* 1  Theoretical  Tunnel  Closure  vs*  Stress  and  Average 

Volume  Change 


. 
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change  and  to  an  intermediate  stage*  Throughout  the 
remainder  of  the  chapter  these  values  of  the  dilation 
parameters  (1,2,3)  will  be  used* 

The  figure  shows  how  much  the  average  volume  change  is 
affected  by  the  use  of  a  different  value  of  dilation 
pa  rame  te  r • 

5  *  3*2*4  Parametric  Study  of  the  Normalized  Closure 

The  normalized  closure  as  previously  defined  depends 
only  on  three  parameters:  C,  0  (strength  parameters),  and  a 
(dilation  coefficient)*  The  purpose  of  this  section  is  to 
evaluate  the  influence  of  changes  in  these  parameters  on  the 
normalized  closure  values* 

Figures  5*2  to  5*4  present  the  results  of  a  sensitivity 
analysis  in  the  form  of  normalized  closure  versus  external 
stress  olots*  The  following  conclusions  can  be  shown: 

1*  The  normalized  closure  increases  with  decreasing 

cohesion  and  angle  of  friction  and  with  increasing 
dilation  coefficient* 

2*  The  yield  point  where  the  curve  deviates  from  1,  is 
directly  proportional  to  the  cohesion  Intercept  C,  but 
it  is  relatively  insensitive  to  changes  in  the  friction 
angl e • 

3.  In  Figures  5*2  and  5*3  the  curves  have  similar  shapes* 
Therefore  several  groups  of  values  ( C,  0)  could  fit 
given  experimental  data*  However  figure  5*2  shows  that 
the  cohesion  is  clearly  determined  by  the  yield  point 
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ELA-PLAST/ELflST 

Figure  5«2  Normalized  Closure  : Sen si t i vi ty  of  the  Cohesion 

Intercept 
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Figure  5*3  Normalized  Closure  JSensiti vity  of  the  Friction 

Angle 
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Figure  5.4  Normalized  Closure  rSensitivity  of  Dilation 
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where  the  normalized  closure  deviates  from  unity* 

5*3*3  Other  Elasto- Plast ic  Models 
5  *  3  *  3  *  1  Introduction 

Two  other  types  of  elasto-plastlc  models  have  been  used 
to  explain  the  behaviour  of  an  opening:  an  elastic  perfectly 
brittle  model  and  an  elastic-strain-weakening  model*  The 
first  one  leads  to  a  closed  form  solution  of  the  problem  and 
the  latter  requires  numerical  treatment*  Only  the  perfectly 
brittle  model  will  be  used  to  derive  a  solution  for  the 
model  test  case*  Published  results  on  the  strain-weakening 
material  indicate  a  oartlcular  shape  for  the  closure  curve » 
which  is  characteristic  tor  this  model  (see  curve  4  on 
Figure  5*5A*  after  Egger(  1973  )  )• 

^*3*3*2  Elastic  Brittle  Plastic  Material 


Th  is  m 

odel  has  been 

used 

to 

derive  expressions 

for  the 

n  o  rma 1 i ze  d  c 

losur e 

i  n  the 

c  as  e 

o  f 

the  model  test  and 

typical 

results  are 

sh  own 

on  f 1 gur  e  5  • 

5  A  • 

The  derivation*  p 

re  sen  t  ed 

in  Appendix  A2  (Section  4)t  does  not  differ  significantly 
from  the  previous  one*  The  failure  criterion  is  the  same  but 
the  olastlc  criterion  relating  the  stresses  in  the  plastic 
zone  is  : 

6J  =  tn6j  +  s6ff  •••  0<s<l 

where  6£  and  are  the  maximum  and  minimum  principal 

st  resses  T  m  the  coefficient  of  passive  earth  pressure  ^  6  ^ 

the  uncon  f 1 ned  compressive  strength  and  s  a  parameter  to 
evaluate  the  strength  drop  after  brittle  failure*  If  s  is 
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ELA-PL AST /ELflST 


NORMALIZED  CLOSURE  VS  STRESS 


E 


Mater ia 1 

B  : Theore  t ic  a 1  Ground  Reaction  Curves  tor  Different 

Materials  (after  Egger(1973)) 
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e*  qua  l  to 
material# 
f  r i c  t i ona 1 
stability 
pressure  , 
has  been 
Panet(  1 97  6 


It  this 
If  s  is 
in  its 
of  the 
and  the 
used  to 
)  • 


is  equivalent  to  a  perfectly  plastic 
equal  to  0,  the  material  is  purely 
ultimate  state*  In  the  latter  case  the 
opening  is  impossible  without  a  support 
tunnel  collapses  immediately*  This  model 
study  rockbursts  in  deep  tunnels  by 


Figure  5*5A  shows  a  discontinuity  in  the  slope  of  the 
stress  closure  curve  immediately  after  yielding*  This  slope 
after  yielding  decreases  with  decreasing  ultimate  strength 
towards  zero  where  the  ultimate  state  is  purely  frictional* 


5*3*3*3  Elaatic’-S  +  ra 

in 

Weakening 

Materials 

F i gu  r e 

5. SB 

in 

d  1  ca  tes 

t  he  results 

obtained  by 

Fsg er(  1973) 

i  n  us  l 

ng 

di f f er  ent 

e  las  to— pie  s  t  i  c 

mode  Is •  They 

are  presented  as  ground  reaction  curves  resulting  from  the 
orogressive  unloading  of  a  rock  mass*  The  closure  curves 
obtained  in  the  model  test  by  progressively  loading  a  tunnel 
have  a  similar  shape*  The  elastic  response  is  different,  but 
as  far  as  the  normalized  closure  is  valid,  the  results  shown 
in  figure  5*5B  may  be  used  as  a  qualitative  description  of 
the  response  of  a  s t ra i n— we aken i ng  material* 

Four  material  models  are  presented  in  Figure  5.5B, 
characterized  by  the  sketch  in  the  upper  right  corner*  The 
s t r a i n— we a ke n i ng  material  has  a  response  intermediate 
between  the  perfectly  plastic  and  the  perfectly  brittle 
materials*  There  is  no  discontinuity  in  the  slope  of  the 
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response  curve  after  yielding*  but  this  slope  decreases 
quickly  to  become  equal  to  the  slope  of  a  perfectly  brittle 
material  response* 

The  response  of  a  strain-weakening  material*  as 
described  here*  requires  a  numerical  treatment*  However*  the 
c 1 osu re— s tre ss  curves  have  the  characteristic  shape 
described  above  and  this  will  be  recognized  if  evident  on 
experimental  data*  If  necessary*  a  numerical  analysis  could 
then  be  carried  out* 


5*.3*.4  Factors  Influencing  Data  Interpretation 

g  *  3*  1*1  Influence  of  the  Longitudinal  Stress 
Throughout  the  derivation  it  has  been  assumed  that  the 
1 on? i tud Inal  prlncl pa 1  stress  is  the  i n te r medl ate  principal 
stress  and  that  the  Coulomb  failure  criterion  holds  in  terms 
of  the  tangential  and  radial  stresses*  It  is  reasonable  to 
assume*  as  did  Florence  end  Schwe r(  1 97 8  )  *  that  at  a  certain 
distance  from  the  tunnel  the  longitudinal  principal  stress 
becomes  equal  to  the  radial  stress  which  increases*  Further 
the  longitudinal  stress  will  become  the  minimum  principal 
stress*  In  the  latter  case,  the  failure  criterion  is  to  be 
taken  between  the  tangential  and  the  longitudinal  stress*  In 
the  intermediate  zone  where  the  two  smaller  stresses  are 
equal  the  stresses  have  to  satisfy  simultaneously  both 
above-mentioned  failure  criteria*  The  failure  surface 


changes  from  the  original  position,  parallel  to  the 


longitudinal  axis* 
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The  analysis  of  this  problem  is  reported  in  Appendix 
A3*2*  It  shows  that  a  value  of  Poisson's  ratio  higher  than 
0*15  is  required  to  avoid  a  change  in  failure  surface*  at 
the  upoer  limit  of  the  loading  range  of  the  machine*  The 
value  of  Poisson's  ratio  measured  in  the  longitudinal 
direction  and  based  on  the  assumotions  of  isotropic  material 
properties  is  about  0*1  (measured  as  a  ratio  of  the  applied 
loads  )  • 

The  oroblem  will  be  more  significant  if  a  support 
pressure  is  applied  on  the  tunnel  wall*  As  there  is  no 
expert  me  ntal  evidence  of  a  change  in  fai lure  directl on  ,  this 
pf feet  has  been  neglected  in  the  present  study* 

5  *  7, 4 «  2  Influence  of  the  Boundaries 

The  sample  block  used  in  the  orevious  analysis  was 
assumed  to  be  infinitely  large  with  respect  to  the  size  of 


the  opening 

•  Howe ve  r  its 

s  i  ze 

is  finite  and  the 

ratio  of 

the 

width  of 

the  sample  to 

♦h  e 

tunnel  diameter  is 

4*80  and 

4*0 

for  the  tests  MC-3  and  MC-4  ,  respectively.  Theoretically  the 
analysis  should  be  corrected  by  assuming  that  the  external 
stress  is  multiplied  by  a  coefficient! 

69  =  69  x  b 2  /(  b 2  -  a2) 

where  69  is  the  external  stress,  a  the  radius  of  the 
opening,  and  b  the  distance  from  the  center  of  the  tunnel  to 
the  boundary  of  the  sample* 

However,  this  only  means  an  increase  in  the  true 
external  stress  of  6*7%  and  4*5%  for  the  test  MC~4  and  MC— 3, 
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respectively 


This 


erro  r 


is  not  significant 


and 
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? «  3  *  4  *  3  Plane  Strain  Condition 

The  expressions  for  the  normalized  closure  have  been 
established  by  use  of  an  elastic  relationship  valid  only  for 
plane  strain  conditions*  The  experimental  data  would  differ 
from  the  theory  If  this  plane  strain  condition  was  violated* 
As  explained  in  chapter  3,  the  loading  head  is  rigid  and 
controlled  by  four  Independent  rams  which  allow  control  of 
the  displacements  in  the  longitudinal  direction  (  Heuer  and 


Hendron  (  1 969  )  ) 

•  It  can  be 

shown  that  the  minor 

Ion  gi tud 1 

na  1 

di solacements 

occurring 

during  the  test 

cause 

an 

i nsi a n i f i c  an t 

error*  This 

is  partly  due  to  the 

1 ow  va lue 

of 

Poisson's  ratio* 

It  is  Important  that  the  operator  maintain  plane  strain 
conditions  within  reasonable  limits  during  the  whole  test* 
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5*.4xi  Method  of  Analysis 

The  sample  was  loaded  in  increments  to 
was  allowed  to  creep  for  a  constant  period 


stages  where  it 
of  t  ime  (  R  days 


. 
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for  test  MC— 9,  1  day  for  test  MC—  4  )«  In  order  to  get  the 
normalized  closure  from  the  experimental  s t ress— st ra i n 
curves,  some  assumptions  about  the  time  —  dependent 
behaviour  of  the  sample  have  to  be  made*  The  following 
aoproach  has  been  adopted  as  a  result  of  analyses  undertaken 
bv  Kaiser(  1979  )  on  the  ti me— dependen t  behaviour  In  the 
pre-failure  zone* 

5-S.4.X.1  jlI.  Theoretical  Model  Response  at  a  Constant  Loading 
Rate 

The  modulus  of  elasticity  of  rock  is  time-dependent  and 
varies  between  two  extreme  limits:  the  upper  limit 
corresponds  to  Instantaneous  loading  (short-term  stiffness ) , 
whereas  the  lower  limit  corresponds  to  infinitely  slow 
loading  (long-term  stiffness)*  If  a  rock  sample  is  loaded  at 
a  finite  rate  to  a  given  stress  level  and  left  under  a 
constant  load,  it  will  creep  until  the  accumulated  strain 
reaches  the  value  which  would  have  been  obtained  during 
infinitely  slow  loading  of  the  sample*  This  is  achieved 
theoretically  after  an  infinitely  long  period  of  time,  but 
in  practice,  most  of  the  creep  strain  occurs  within  the 
first  day*  The  amount  of  creep  strain  is  a  function  of  the 
rate  of  loading  and  of  the  stress  increase  during  the  latter 
loading  (comparable  to  a  vi sco— elastic  material)* 

Assuming  limits  for  the  modulus  of  elasticity  and  a 
value  of  Poisson's  ratio,  the  corresponding  limiting  curves 
for  the  closure  can  be  determined*  This  has  been  done  in 


' 


* 

r‘- 


75 


F 1  sure  5*6*  From  the  elastic  curves*  which  are  assumed  to  be 
linear  (dashed  lines)*  the  corresponding  e las t o— plasti c 
limits  have  been  obtained  by  muliplylng  the  elastic  closure 
by  the  value  of  the  normalized  closure  at  a  given  stress 
level*  Reasonable  elastic  and  strength  parameters*  as 
indicated  on  Figure  5*6*  have  been  assumed* 


5  *  4 » 1*2  Theoretical  Model  Response  to  Step  Loading 

A  theoretical  closure  curve  for  the  stress  history  of 
the  model  test  may  be  plotted  within  this  framework*  The 
sample  is  loaded  up  to  the  first  Increment  along  a 
constant-rate  loading  curve*  In  order  to  simplify  the 
exolanat 1 ons *  this  rate  of  loading  is  assumed  to  be  the 
Instantaneous  loading  •  Then  the  sample  creeps  along  a 
horizontal  line  until  it  reaches  the  long-term  loading 
curve*  It  is  reloaded  to  the  next  increment  following  a  line 
parallel  to  the  Instantaneous  e last o-pla s ti c  loading  curve 
and  this  cycle  repeats  up  to  the  upper  limit  of  the  loading 
range*  It  can  be  further  assumed  that  the  material  behaves 
as  a  linear  visco-elastic  material  and  then  the  amount  of 
the  creep  deformation  at  any  stage  is  proportional  to  the 
previous  stress  increment* 


5  *  4* l . 3  Application  to  Experimental  Model  Test  Responses 

As  a  consequence  of  the  above  analysis,  it  is  assumed 
that  the  curve  which  is  obtained  by  Joining  the  points 
corresponding  to  the  ends  of  the  creep  tests  is  the 
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TUNNEL  CLOSURE  VS  EXT,STRESS 


Figure  5.6  Theoretical.  Model  Response  for  an  Elasto-Plastic 

Linear  Visco-Elastic  Element* 
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Ion 

calculation  of 

t  he 

later  Figures 

5*7, 

reloading  curve 
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curve,  this  latter  V 


—term  loading  curve* 
initial  portion  of 
g— term  loading  curve 
experimental  normal 
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A  straight  line 
this  curve  is  the 
This  enables  a 
zed  closure  (see 
assuming  that  the 
antaneous  loading 


5*.4*.2  Experimental  Measurements 

The  original  tunnel  closure  measurements  for  the  two 
tes  +  s  have  been  reproduced  in  Figures  F*7  and  5*8*  The 
sketch  in  the  upper  right  corner  indicates  the  orientation 
of  the  joints(  line  portions)  and  the  relative  direction  of 
the  instrument*  From  this  experimental  step— like  curve  the 
e l as t o-p las t  ic  (axis  line)  and  the  constructed  linear 
elastic  ( dashed  line)  long— team  loading  curves  have  been 
determined  in  the  manner  explained  in  the  last  section*  This 
latter  line  intersects  the  abscissa  at  a  strain  which 
corresponds  to  the  strain  necessary  to  close  the  cracks  in 
the  rock  mass  and  to  allow  some  initial  seating  of  the 
instruments*  This  crack  closure  is  deducted  from  the  total 
strain  values  to  calculate  the  normalized  closure,  by 
dividing  the  reduced  value  of  the  elasto-plas t i c  closure  by 
the  reduced  value  of  the  elastic  closure  at  each  stress 
level  (see  section  5*3*4)* 

The  second  solid  line  on  each  plot  corresponds  to  the 
constructed  instantaneous  loading  curve*  Its  shape  confirms 
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Figure  5*7  Test  MC— 3  t  Experimental  Tunnel  Closures 
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Figure  5.3  Test  MC-4  :  Experimental  Tunnel  Closures 
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fhe  validity  of  the  method  of  analysis#  It  has  been 
constructed  in  the  following  manner  (the  references  in  the 
brackets  refer  to  figure  5*6)* 

For  a  linear  vi sc o— e 1 a st 1 c  material  the  amount  of  creep 
strain  Is  proportional  to  the  latter  stress  increment#  The 
amount  of  creep  strain  occurring  in  the  second  creep  test 
(cd)  corresponds  to  a  stress  increment  of  2  MPa#  It  is  used 
to  calculate  what  should  theoretically  be  the  amount  of 
creep  strain  in  the  first  stage  of  loading  corresponding  to 
a  F  MPa  stress  Increase  (  ab)#  The  first  point  (a  )  of  the 
sh  ort  — term  loading  c  urve  is  de  te  rm ined#  It  is  not  a  point  of 
the  experimental  curve  as  some  creep  already  occurred  before 
the  completion  of  the  first  load  increment#  Assuming  that 
the  reloading  curve  between  stages  is  parallel  to  the 
el asto-plast ic  short-term  loading  curve  (in  other  words 
ab=“c  )  the  next  point  (e)  of  this  curve  may  be  determined# 
The  same  process  Is  reoeated  to  get  the  complete  short-term 
loading  curve# 


F  #  4 «  3  Limitations  of  the  Method 

The  actual  shape  of  the  estimated  loading  curves  show 
some  inconsistencies  which  Indicate  that  the  method  of 
analysis  is  limited  by  several  restrictive  assumptions,  such 
as  linear  elasticity,  linear  visco-elasticity,  and 


uniformity  of  loading# 


, 


, 


81 


5  «  4  * 3  « 1  Linear  Elastic  Reference  Line 

A  linear  elastic  reference  line  was  assumed  in  order  to 
calculate  the  normalized  closure*  However,  it  is  evident 
from  the  l oading  tests  on  the  sample  without  the  tunnel  that 
this  assumption  is  not  completely  satisfied*  The  elastic 
properties  of  the  material  tend  to  change  with  increasing 
stress  in  such  a  way  that  the  coal  becomes  less 
c ompress 1 ble •  A  linear  elastic  reference  line  leads  to  an 
o ve re s 1 1 ma 1 1  on  of  the  elastic  closure  and  an  underestimation 
of  the  normalized  closure*  Linearity  is  not  a  criterion  for 
elasticity;  a  straight  line  could  be  Interpreted  as  the 
yielding  line  of  a  non— linear  stiffening  elastic  line* 

Cn  the  other  hand,  the  elastic  line  is  based  on  the 


first 

two  creep  tests 

where  the  material  is 

assumed  to 

behave 

elastically*  If 

this  is  true  for  the  first 

stress 

stage, 

•♦•his  assumption 

may  be  questionable  for 

the 

second 

creeo 

test  at  7  MPa* 

F  *  4  »  3  *  2  Tlme-’Dependent  Strength 

The  method  supposes  that  the  normalized  closure  is 
t *  me- 1 ndependent  and  that  the  strength  parameters  are 
constant  throughout  the  test*  However, coal  has  a 
visco-plastic  behaviour  (Kaiser  and  Morgenstern,  1979);  the 
strength  parameters  and  particularly  the  cohesion  intercept 
C  decrease  with  time,  and  as  explained  in  the  sensitivity 
analysis,  the  normalized  closure  increases  for  a  given 
stress  level*  The  experimental  normalized  closure  has  been 
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calculated  from  Ion?— term  loading  curves  where 
parameters  have  their  lower  but  constant 
analysis  of  the  data  should  give  a  lower 
strength  oarameters* 


the  strength 
values*  The 
bound  of  the 


F  .  4 . 3  «  3  Uniformity  of  Loading 

The  solution  Is  based  on  the  assumption  of  isotropic 
loading  conditions*  However^  if  the  N— value  is  kept  sensibly 
cons+ant  during  each  load  stage,  it  is  difficult  to 
reoroduce  the  same  N— value  at  different  load  stages  and  also 
to  get  the  N— value  exactly  equal  to  unity*  This  analysis  of 
the  data  appears  to  be  extremely  sensitive  to  a  change  in 
N-value  at  different  load— stages*  The  elastic  reference  line 
used  to  calculate  the  normalized  closure  is  no  longer  linear 
as  the  elastic  response  is  a  function  of  N*  The  upper  right 
plot  on  Figure  5*3  is  a  good  example  of  this  feature*  The 
linear  elastic  line  crosses  the  "plastic"  line  which  is 
i mpossi bl e • 

If  N  differs  from  l ,  the  concept  of  normalized  closure, 
indeoendent  of  the  elastic  constants,  is  no  longer  of  value* 
But  as  the  variation  of  N  from  unity  is  small,  values  of  the 
normalized  closure,  corrected  for  a  change  in  N— value  at 
different  load  stages  have  been  calculated. The  method  is 
exolained  in  detail  in  Append! x  A3* 1  •  It  is  assumed  again 
that  the  first  two  points  of  the  experimental  closure  curve 
are  in  the  elastic  range*  Their  secant  moduli  result  from 


the 


multiplication  of  a  function  A  of  the  elastic  constants 
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by  a  function  of  N  depending  on  the  relative  orientation  of 
the  gauge  and  the  maximum  applied  stress#  To  close  the 
cracks  an  original  "closure  strain"  Gel  is  necessary#  The 
secant  moduli  are  calculated  by  assuming  a  translation  of 
the  origin  of  value  Bel  as  in  the  conventional  method#The 
knowledge  of  the  N— value  at  the  first  two  load  stages  will 
permit  the  calculation  of  the  values  of  the  unknown  A  and 
Pci  for  each  curve#  Knowing  these  values  and  knowing  N  at 
each  subsequent  stress  stage,  it  is  possible  to  calculate  a 
reference  elastic  closure  for  each  stress  level  and  then  to 
calculate  the  normalized  closure#  The  application  of  this 
method  to  the  experimental  data  (lower  plot  on  Figure  5#9 
and  5*10)  proved  to  be  successful  in  reducing  the  scatter  of 
the  results# 


5  »  4  #  4  Interpretation  of  Experimental  Results 

Figures  5# 9  and  5#10  present  the  experimental  results 
as  dashed  lines  for  the  tests  MC-3  and  MC-4  respectively# 
The  upper  plots  correspond  +o  data  directly  interpreted  from 
Figures  5#7  and  5.8#  Only  three  measurements  are  given;  in 
each  case  one  of  the  diameters  shows  a  normalized  closure 
smaller  than  i  and  is  not  reoresented  on  the  figure.  The 
lower  plots  show  the  data  corrected  for  a  varying  N-value 
according  to  the  method  noted  in  the  previous  section#  The 
curves  are  smoother  and  the  results  become  more  consistent# 
The  comparison  between  the  experimental  data  and  the 
results  of  the  theoretical  analysis  given  in  section  5#3.2, 
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Figure  5« 9  Test  MC— 3  : Experimental  Normalized  Closure 
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Figure  5.10  Test  MC— 4  2  Experimental  Normalized  Closure 
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has  been  undertaken  In  the  following  manner* 

1*  A  computer  code  permits  the  visualization  on  a  graphic 
terminal  of  theoretical  elasto— plast i c  responses  of  the 
model  with  any  set  of  parameters  C,  0,  a ,  s*  The 
analysis  consists  of  comparing  plots  of  the  experimental 
data  with  different  theoretical  curves  for  changes  in 
parame  ters • 

7m  The  cohesion  intercept  C  will  vary  around  4  MPa,  as  it 
has  been  shown  that  the  response  is  sensitive  to  a 
change  in  C*  Moreover,  as  the  analysis  has  to  take  into 
account  the  relative  orientation  of  the  joints  in  the 
failing  element,  C  is  an  appropriate  parameter  with 
which  to  study  variation  due  to  structural  anisotropy* 

3*  The  angle  of  friction  0  is  kept  equal  to  30°,  which  is 
the  value  given  by  trlaxlal  tests,  as  it  does  not  affect 
substantially  the  e last o— o las t 1 c  response* 

4*  The  dilation  parameter  a  is  taken  equal  to  1  (no  volume 
change)  or  to  3  (associated  flow  rule  if  0  equal  to  30° 
)•  This  is  done  in  order  to  account  for  the  theory  found 
in  the  literature,  by  which  one  may  expect  the  material 
to  dilate  according  to  the  associated  flow  rule  only  to 
a  limited  oortlon  of  the  pos t-f ai lure  strains  and  then 
to  maintain  a  constant  volume*  (Ladanyl  ,1974  ) 

5.  The  strength  loss  parameter  s  is  given  the  values  1,  0.8 

end  0*6*  The  three  curves  then  obtained  will  give  an 
indication  of  what  would  be  the  response  of  a 
strain-weakening  material  with  a  progressive  strength 
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toss# 

This  analysis  leads  to  the  following  conclusions  about 
e  e last o-plas t § c  closure  of  an  opening*  These  conclusions 
e  In  accordance  with  the  generally  accepted  practice  in 
e  Interpretati on  of  field  data*  Right  after  the  onset  of 
elding  the  material  dilates  according  to  the  associated 
ow  rule*  This  effect  is  dominant  but  tends  to  decrease 
en  yielding  proceeds,  until  the  material  no  longer  shows 
lume  change*  This  reduction  in  volume  change  may  be 
plained  by  some  arching  after  dilation  has  taken  place* 
is  would  explain  why  the  curves  B  and  C  on  Figure  5*9  show 
change  in  curvature  and  tend  to  exhibit  less  closure*  At 
e  same  time,  as  the  plastic  strains  increase,  the 
rain-weakening  material  loses  progressively  its  strength 
d  more  deformations  occur*  The  strain-weakening  effect 
come  orogressively  more  dominant* 

This  scheme  explains  individually  every  curve  on 
gures  5*9  and  5*10  and  the  following  conclusions  can  be 


awn  : 

Test  4C-4  results  (where  the  tunnel  ruptured) 
dominated  by  the  strain-weakening  character 
material  than  test  MC— 3  results* 

The  D— curves  presenting  the  closure  in  the  d 
pe r pendi cula r  to  the  joints  are  much  more  aff 
the  strain-weakening  effect  than  the  other  measu 
The  anisotropy  in  structure  and  in  strengt 
argued  to  explain  the  different  behaviour  of 
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gauges  in  the  same  test. 

4.  The  strength  parameters  found  to  fit  the  experimental 
curves  are  In  accordance  with  our  knowledge  of  the  coal. 
The  cohesion  intercept  may  be  considered  as  high  with 
respect  to  the  values  obtained  in  triaxlal  tests.  This 
can  be  explained  as  an  effect  of  the  size  of  the 
rupturing  element  which  is  much  smaller  in  the  tunnel 
wall  than  in  a  triaxial  test. 

5.  This  size  or  shape  effect  which  influences  also  the 

post— failure  behaviour  (see  Chapter  6  Fig.  6.8)  may 
explain  why  the  test  MC— 4  with  a  bigger  opening  and  a 
larger  rupturing  element  may  exhibit  a  larger  strength 
loss  for  the  same  value  of  the  strain  than  the  test 
MC— 3.  (  See  Chapter  6  for  more  details  about  the  shape 

effect  on  the  stability  ) 

However*  this  interpretation  does  not  explain  a  major 
inconsistency.  In  both  tests*  the  diameter  C  parallel  to  the 
joints  shows  much  more  closure  than  the  diameter  D 
perpendicular  to  the  joints.  The  above  interpretation 
indicated  an  effect  of  the  anisotropy.  But  if  this  was  true, 
the  cohesion  intercept  would  have  to  be  much  lower  along  the 
C-dlame+er  (where  the  joints  are  perpendicular  to  the 
maximum  principal  stress)  than  along  the  D— diameter. 
According  to  common  knowledge  about  the  anisotropic  strength 
of  a  jointed  mass,  this  cannot  be  true.  Moreover  visual 
observations  on  the  samples  show  that  in  particular  on 
sample  MC-4,  the  tunnel  wall  was  intact  in  the  C-dlrection 
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after  unloading  and  the  large  closure  cannot  be  attributed 
to  yielding  only.  This  serious  inconsistency  occurs  in  both 
tests,  showing  that  some  assumptions  on  which  the  theory  is 
based  are  not  valid. 

Because  of  the  heterogeneous  and  t  ime-dependent 
characteristics  of  the  material,  it  is  reasonable  to  assume 
that  the  discrepancy  between  the  data  and  e lasto— plast i c 
theory  is  due  to  non— ax isy mmetr ical  yielding  conditions  and 
to  an  anisotropic  time-dependent  stress  redistribution. 
Simplifying  assumptions  were  essential  in  order  to  derive 
closed— form  solutions  to  analyse  this  problem  and  only 
numerical  analyses,  which  are  beyond  the  scope  of  this 
thesis,  will  give  a  quantitative  answer  to  the  problem  with 
thi=  increased  complexity. 

Kaiser  (1979)  analysed  the  results  of  a  similar  test 
and  studied  in  detail  the  time-dependent  pre-failure 
behaviour.  He  was  able  to  explain  the  t i me- independent  and 
lime-dependent  behaviour  by  s t re ss- redi s tr i but i on  due  to 
non— linear  creep  properties,  non— i so t ropi c  creep  properties 
or  local  yielding  and  argued  that  an  "equivalent  opening" 
approach  can  be  used  to  describe  the  observed  non  linearity 
of  the  c 1 osu re— s tre ss  curves.  This  work,  in  connection  with 
the  conclusions  from  this  chapter,  leads  to  the  following 
i nt eroretati on  of  the  failure  process  in  the  model  test: 

1.  elastic  deformations  (more  or  less  isotropic)  , 

2.  anisotropic  stress  redistribution  (creep), 

3.  anisotropic  stress  redistribution  (yielding), 
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4.  anisotropic  rupture. 

Each  of  the  last  three  phenomena  makes  the  tunnel  behave 
like  an  elliptical  opening.  In  the  cases  2  and  3  the 
elliptical  opening  is  only  apparent  while  in  case  4  the 
opening  actually  become  elliptical. 

In  the  light  of  this  interpretation  a  simple  analysis 
of  the  experimental  data  was  carried  out  by  use  of  charts 
given  by  Feder  (1978)  for  the  elastic  deformations  of  an 
elliptical  opening  (Figure  5.11).  The  particular  expressions 
for  the  closure  of  the  model  test  are  indicated  on  the 
Figure  5.11.  These  charts  provide  only  the  closure  of  the 
cavl ty  in  the  direction  of  the  axes  of  the  ellipse  and  the 
experimental  measurements  are  not  necessarily  along  the  axes 
of  the  "apparent"  ellipse.  The  application  of  these  charts 
shows  that  the  curves  C  and  A  on  Flsure  5.9  and  5.10  may  be 
i n t erore ted  as  +he  elastic  closure  of  an  elliptical  cavity 
along  its  small  axis*  while  an  istropic  e las t o-p la s t 1 c 
analysis  interprets  them  as  substantially  influenced  by  a 
plastic  behaviour.  The  curves  B,  and  in  particular  Dt 
correspond  to  the  minor  deformations  along  the  long  axis  of 
a  similar  elliptical  cavity,  under  uniform  loading 
conditions  (N=l).  They  are  influenced  by  yielding  effects. 

It  appears  then  that  the  deviation  from  linearity  of 
the  tunnel  closure  measurement  curves  are  not  essentially 
due  to  the  yielding  of  the  material  but  may  be  due  to  other 
effects  associated  with  the  anisotropic  character  of  the 
material.  As  shown  above,  the  interpretation  of  closure  data 
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Figure  5*11  Elastic  Deformations  of  an  Elliptical  Cavity 

( after  Feder( 1978 ) 


Radial  displacement  of  inside  edge  of  compact  area.  Shown  equation  for  n  refers  to  roof 
displacement.  Calculation  of  wall-displacement  requires  exchanging  the  position  of  F.y  and 

Fir  and  using  the  right  hand  charts 
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by  mean  of  elasto-plastlc  methods  is  apparently  possible  but 
may  lead  to  incorrect  conclusions. 

fLaJE  Summary  and  Conclusions 

In  order  to  explain  the  non-linear  behaviour  in  the 
experimental  closure— stress  curves*  several  elasto-plastic 
models  have  been  adapted  to  derive  the  theoretical  value  of 
the  tunnel  closure  in  the  model  tests.  These  models  are  used 
commonly  in  the  interpretation  of  field  measurements.  They 
aret 

1.  elastic  perfectly  plastic  model  (  c losed— f or n  solution) 

2.  elastic  brittle-plastic  model  (closed-form  solution) 

3.  elastic  strain-weakening  model  (after  Egger  (1973)  and 
Borsetto  (1979)). 

For  isotropic  loading*  the  concept  of  normalized  closure 
(ratio  of  elasto-plastic  over  purely  elastic  closure)  has 
been  introduced  and  the  analysis  of  the  sensitivity  of  the 
strength  oarameters  C  and  0  and  the  dilation  coefficient  at 
on  the  closure*  leads  to  the  following  conclusions  : 

1.  The  closure  Increases  when  decreasing  the  cohesion  C, 
when  decreasing  the  friction  angle  0  and  when  increasing 
the  dilation  parameter  a. 

2.  The  yield  stress  where  the  closure-stress  curve  deviates 
from  linearity  is  proportional  to  the  cohesion  C  but  is 
not  particularly  affected  by  the  other  parameters. 

3.  The  use  of  the  perfectly  brittle  model  causes  the  curve 
to  present  a  discontinuity  in  the  slope  right  after  the 
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onset  of  yielding. 

4.  A  s t rai n— weaken Ing  model  shows  an  intermediate  response 
between  a  perfectly  plastic  model  and  a  perfectly 
brittle  model  and  may  be  recognized  in  this  manner. 

5.  The  influence  of  several  factors  such  as  plane  strain 
conditions*  isotropic  loading*  and  sample  boundaries  has 
been  recognized  and  their  effect  has  been  either 
neglected  or  taken  into  account  in  the  analysis. 

The  interpretation  of  the  experimental  data  after  this 

theoretical  e lasto-pl as t ic  analysis  leads  to  the  following 

conclusions! 

1.  The  closure  of  a  tunnel  seems  to  be  dominated  by  two 

superposed  effects  :  right  after  failure  by  a  volume 

change  following  the  associated  flow  rule  and  whose 
importance  decreases  with  increasing  plastic  strains  and 
by  a  strain-weakening  post-failure  behaviour  which  is 
sensitive  after  some  plastic  strains  occur. 

2.  The  influence  of  the  strainrweakening  is  more  sensible 
in  test  MC-4  where  the  diameter  of  the  opening  is 
larger .This  may  be  exolained  by  the  size  effect  of  the 
rupturing  element  on  the  strength  and  on  the 
oost— failure  behaviour  of  this  element. 

3.  In  the  light  of  this  interpretation  each  individual 
normalized  closure  curve  may  be  fitted  by  a  theoretical 
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intercept )  • 

4*  However  the  group  of  four  measurements  in  the  tunnel 
cannot  be  consistently  interpreted  in  such  a  manner,  and 
the  discrepancy  is  too  Important  to  be  attributed  to 
some  error  in  the  exoerlmental  procedures  or  in  the  data 
analysis*  In  particular  the  measurements  of  the  closure 
along  the  diameters  parallel  to  the  joints  show  a  major 
deviation  from  linearity  as  visual  observations  indicate 
an  intact  material*  Even  if  it  appears  so,  the  closure 
of  the  tunnel  cannot  be  completely  explained  by  an 
isotropic  elasto— plasti c  analysis* 

5*  The  original  assumptions  of  an  isotropic, 

t i me— i ndependen t  material  were  suspected  to  be  not  valid 
and  contribute  to  the  discrepancy  between  theoretical 
and  experimental  closures*  These  assumptions  are, 
however,  essential  for  the  derivation  of  closed-formed 
solutions  and  if  they  are  not  valid  the  solutions 
require  a  numerical  analysis  which  has  to  be  developed 
for  this  purpose* 
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non-linearity  in  the  stress-closure  curves.  A  brief 
analysis  conducted  on  the  experimental  data  in  the  light 
of  this  hypothesis  show  reasonable  agreement  between  the 
data  and  the  closure  of  an  apparent  or  real  elliptical 
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The  most  important  conclusion  of  this  chapter  is  that 
even  if  an  isotropic  elasto-plastlc  analysis  predicts  well 
the  deformations  of  a  tunnel*  it  may  lead  to  incorrect 
conclusions  about  the  failure  process*  An  anticipation  of  a 
global  yielding  around  the  cavity  leads  to  an  Incorrect 
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stress  redistribution  processes  has  to  be  analysed  and 
evaluated  In  future  work* 


CHAPTER  6 


ANALYSIS  OF  THE  RUEXLLEE  PROCESS 

6  «  1  Intrgdugt  jog 

^<1*1  Scope  of  the  Chapter 

Chapter  5  concludes  that  the  recognition  of  the  failure 
mechanisms  is  essential  to  develop  a  design  criterion  for 
the  support  of  a  tunnel*  The  advantage  of  model  tests 
performed  on  real  material  over  a  case  history  study  or  a 
t°st  on  artificial  material  Is  that  it  allows  us  to 
instrument,  control  and  observe  the  failure  mechanism  of  an 
opening,  in  a  real  material,  under  known  stress  conditions* 
However  the  laws  of  similitude  are  particularly  difficult  in 
the  inelastic  range  of  the  material,  and  as  explained  In 
Chapter  3,  they  have  to  be  carefully  investigated  before  the 
conclusions  obtained  on  the  model  could  be  applied  to  real 
ca vi t i es  • 

In  this  chapter,  the  tunnel  failures  observed  during 
the  tests  are  described*  In  particular,  the  actual  slip 
lines  are  described  in  the  context  of  a  plastic  equilibrium 
analysis.  An  hypothesis  on  the  rupture  initiation  and 
pi'ooa  ge  1 1  on  is  given  and  the  role  of  preexisting 
d! scont inu  it ies  Is  emphasized* 

6*1.2  Nomenclature 

In  acccordance  with  B  1  en i a wsk 1 (  1967  )  ,  the  following 

nomenclature  will  be  used*  The  material  Is  linear  elastic  up 
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to  the  yield  point.  As 
irreversible  deformations 
characterized  by  a  loss  of 
Rupture  is  a  failure  which 
the  rock  element  in  sever 
that  the  behaviour  of  a  roc 
and  rupture  may  be  controll 
enough.  In  a  rock  mass 
conditions*  Intact  element 
play  the  role  of  the  1 
rupture.  Rupture  depends  on 
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LjlZ  Description  of  the  Observed  Syfiturgg 

6.2. 1  Sample  MC-3  after  Unloading 

As  noted  in  the  test  history  (Chapter  4),  this  te 
Interrupted  when  a  corner  of  the  sample  sheared  o 
rupture  was  apparent  on  the  tunnel  wa 1 1  dur i ng  the  tes 
block  was  taken  out  of  the  compression  machine.  Pla 
«=h ows  the  crack  pattern  on  the  bottom  face  of  the  s 
The  sheared-off  corner  Indicates  the  direction  of  the 
as  the  sample  failed  along  a  steplike  plane  parallel 
joints. 
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plastic  strains  are  created  in  the  plaster  by  loading*  and 
they  cause  tensile  stresses  during  unloading  as  the  coal 
Imposes  the  displacements# 

The  crack  pattern  is  randomly  distributed  near  the 
boundaries  and  is  radial  around  the  tunnel#  The  assumption 
of  a  uniform  radial  stress  field  is  therefore  supported#  The 
extent  of  the  cracks  Is  a  function  of  the  state  of  stress 
reached  in  underlying  coal  elements •  If  stress 
redistribution  occurs  during  the  test*  the  stress  level 
increases  first  and  decreases  subsequently  on  each  element 
around  the  tunnel#  The  plaster  layer  will  "memorize"  this 
stress  increase  as  an  Irreversible  strain#  A  circular  crack 
pattern  will  indicate  a  uniform  stress  distribution*  and  an 
elliptical  zone  will  show  a  non— uniform  stress 
redistribution*  and  will  give  the  preferential  direction  of 
+hls  redistribution# 

A  careful  observation  of  the  crack  pattern  presented  in 
Plate  6#1  Indicates  a  vertical  preferential  direction  for 
stress  redistribution#  This  corresponds  to  the  behaviour 
recorded  by  the  Instruments  and  described  in  the  last 
section  of  Chapter  5#The  other  major  cracks  correspond  to 
preexisting  faults  described  in  Chapter  2# 

The  plaster  covering  the  sample  and  the  latex  film 
which  prevented  the  tunnel  wall  from  drying  were  removed# 
Some  debris  remalnded  attached  to  this  membrane,  leaving 
cavities  and  cracks  on  the  tunnel  wall#  Some  slip  lines 
opened  up*  indicating  that  even  If  the  tunnel  did  not 
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Plate  6*1  Sample  idC— 3  after  unloading  ICverall 
bottom  face  and  close-up  of  the  z. one  around 


view  of  the 
the  tunnel 
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rupture  9  a  small  plastic  zone  existed*  When  compressed 
locally  by  finger,  the  "Intact”  tunnel  wall  has  a  friable 
appearance,  indicating  a  high  level  of  mlcrof racturing  In 
the  neighbourhood  of  the  tunnel  wall* 

The  second  picture  in  Plate  6*1  is  a  close-up  of  the 
tunnel  wall  and  of  the  jointed  material  near  It*  Some  joints 
are  open,  and  in  the  direction  perpendicular  to  them, 
tensile  cracks  propagate  through  the  sample*  It  is  difficult 
however  to  differentiate  the  preexisting  cleat  system  of  the 
coal  and  cracks  due  to  over stress i ng*  Inside  the  tunnel, 
some  features  presented  in  Figure  6*1  are  recognizable  (near 
axis  4  of  the  figure)* 

Figure  6*1  shows  a  complete  sketch  of  the  tunnel  wall; 
a  shaded  area  indicates  a  cavity;  the  curved  smooth  lines 
are  the  traces  of  preexisting  faults;  in  the  middle  of  the 
sample,  there  Is  the  trace  of  a  major  bedding  plane  which 
was  noted  in  the  core  descrlotlon*  Most  of  the  depressions 
in  the  wall,  resulting  from  rock  fall,  are  located  near  one 
of  the  major  d i sc o n t in u i t 1 es •  This  confirms  the  importance 
of  such  discontinuities  on  the  stability  of  a  tunnel, 
particularly  with  respect  to  rock  fall* 

The  shearing  of  the  corner  is  due  to  a  local  change  in 
stress  conditions  after  a  si ight  rotation  of  the  sample*  The 
measurements  in  the  vicinity  of  the  tunnel  were  not  affected 


by  this  failure* 
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Figure  6.1  Sample  MC-3  alter  Unloading  :  Wap  of  the  Tunnel 

Wall 
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6 « 2 a_ 2  Description  of  the  Rupture  in  Sample  MC— 4 
t>  *  2  *  2  »  1  Rupture  History 

The  sample  was  loaded  in  stages  as  described  in  the 
test  history  (Chapter  4)*  The  rupture  Initiation  occurred 
about  two  hours  after  the  sample  had  been  loaded  to  the 
increment  of  13*3  MPa •  Then*  as  shown  on  the  strain 
rate— time  curves  (Appendix  B2),  the  rupture  proceeded  by 
waves,  until  the  process  stabilized  (see  next  section)* 

Deformations  of  the  tunnel  wall  have  been  observed  at 
mid— height  in  the  tunnel,  on  both  ends  of  a  diameter 
inclined  at  about  40  0  to  the  vertical  axis,  close  to  the 
points  where  the  joints  are  tangent  to  the  circular  tunnel* 
A  square  area  (2cm  by  2cm)  of  the  tunnel  wall  is  broken  with 
a  V-shape,  pointing  towards  the  cavity* 

After  a  creeo  test  of  28  hours,  the  failure  became 
stabilized  and  the  sample  was  reloaded  to  the  next  increment 
(15*9  MPa)*  Failure  proceeded  Immediately  after  loading; it 
was  again  characterized  by  the  wave— like  aspect  of  the 
strain  rate-time  log— log  olots*  The  ruptures  extended  not 
only  in  the  horizontal  but  also  in  the  longitudinal 
direction.  The  debris  contained  by  the  latex  film  was  so 
loose  that  the  points  of  the  LVDT*s  penetrated  into  it,  and 
It  then  occupied  a  large  volume  of  the  cavity*  The 
measurements  of  the  closure  in  the  direction  perpendicular 
to  the  joints  are  thereafter  questionable,  not  only  because 
of  the  crushing  of  the  wall,  but  also  because  of  the  large 
displacement  values,  beyond  the  linear  range  of  the 


< 
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electronic  transducers.  After  a  24  hour  creep  test,  the 

sample  had  to  be  unloaded  because  it  began  to  rotate# 

After  recovery,  the  sample  was  uncovered,  the  debris 
cleaned  out  and  the  measurement  device  was  removed#  The 
sample  at  this  stage  is  shown  on  Plate  6.2#  On  these 
pictures,  several  details  should  be  noted  : 

1#  The  rupture  occurred  in  a  direction  between  the  vertical 
and  a  direction  inclined  at  45  °#  This  deviation  is  most 
likely  due  to  the  error  in  maintaining  isotropic  loading 
at  the  beginning  of  the  test#  The  actual  value  of  this 
angle,  measured  on  the  cleaned  tunnel,  is  35°  from  the 
vertical  direction# 

2#  The  rupture  Is  more  apparent  on  one  side  of  the  tunnel; 

more  displacement  occured  in  this  location# 

3#  The  plaster  covering  the  sample  followed  the 

displacements  of  the  broken  mass  and  cracked  along  lines 
whose  shape  approaches  a  logarithmic  spiral# 

4#  The  close-up  of  the  failed  area  shows  clearly  in  the 

debris  broken  "bands"  of  coal#  The  Joints  are  still 
parallel,  but  the  whole  mass  is  rotated  towards  the 
inside  of  the  cavity# 

5#  Another  remark  concerns  the  crack  in  the  plaster  which 
is,  on  the  picture,  the  closest  to  the  opening#  The 
edges  of  the  crack  are  displaced  and  show  evidence  of  a 
sliding  movement  of  the  underlying  rock#  This  movement 
Is  directed  towards  the  bottom  of  the  failed  cavity# 
Feder(1978)  proposed  a  movement  In  the  other  direction 
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MC— 4  after  the  First  Rupture  : 
Tunnel  Wall  and  Debris 


Plate  6*2  Sample 


fiupt ured 
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(  towards  the  Intact  part  of  the  wall)  in  his  analysis  of 
an  elliptical  failed  cavity* 

6*  Because  of  this  sliding  movement*  the  rupture  is  not 


only 

a  splitting  phenomenon* 

The 

section  of  the  broken  tunnel,  presented  on  Figure 

6*3*  is 

more  or  less  e 11 ipt leal; the  ratio  between  the 

dimensions  of  the  axes  is  1*7*  The  value  of  this  ratio  is  in 
accordance  with  the  estimated  value  presented  in  Chapter  5; 
it  is  smaller  than  the  calculated  ratio  for  the  equivalent 


o  oen ing* 

but  only  the  material  in  its  ultimate  state  was 

r  pmoved * 

leaving  a  zone  of  weakened  material  around  the 

c  a vi ty • 


The 

tunnel  walls  were  protected  by  another  latex 

coat ing 

and  the  sample  was  loaded  by  increments  of  2  MPa, 

after  an  initial  increment  of  7  MPa*  Evidence  of  failure 
apDeared  after  loading  to  11  MPa ;  the  latex  membrane 
containing  the  broken  material  swelled  again  into  the 


c  a  vi  ty  • 

Loading  proceeded  by  steps  up  to  16*6  MPa  where  the 

rupture 

process  seemed  to  terminate  and  the  sample  was 

unloaded  and  left  for  a  recovery  test. 

The  new  cross-section  of  the  opening  is  shown  on  the 
Figure  h.7.  It  seems  that  the  rupture  started  again  from  the 
bottom  of  the  failed  zone  and  created  an  elliptical  cavity 
propaaating  radially  outward  from  this  point* 

Pictures  of  the  sample  after  the  plaster  has  been 
removed  are  shown  on  Plates  6*3  to  6*5*  Plate  6*3  presents 


the  general  aspect  of  the  sample  after  unloading*  before  and 


* 
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Figure  6«2  Sample  MC— 4 


Cross-Sections  of  the  Ruptured 
Tunne 1 


. 
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after  the  debris  was  removed*  The  overall  view  shows  that 
there  are  two  major  sets  of  cracks,  on  each  side  of  the 
ellipse,  along  joints,  proceeding  towards  the  upper  left 
corner*  The  sample  would  have  finally  collapse  by 
” through  —  s tampl ng"  of  the  cavity  roof*  Feder(  1 978  )  conducted 
similar  tests  on  artificial  material  and  described  such  an 
ultimate  collapse  of  the  sample*  The  internal  strain  gauges 
now  visible  inside  the  broken  tunnel  are  in  the  45° 
direction,  and  the  deviation  of  the  major  axis  of  the  new 
elliptical  section  from  this  direction  Is  apparent*  On  both 
oictures  of  Plate  6*3,  the  set  of  silo  lines  and  the  manner 
in  which  it  Is  affected  by  the  joints  is  also  apparent!  see 
further  Fig*  6*6A)*  The  following  explanations  refer  not 
only  to  Plate  6*3  but  also  to  Plate  6*4  and  Figure  6*3*  A 
set  of  tensile  cracks  approximately  oarallel  to  the  free 
tunnel  surface  can  be  observed*  On  one  side  the  debris  is 
highlv  crushed;on  the  other  side,  a  whole  broken  lump 
remains  Intact*  Actually  this  block  failed  by  shearing  along 
the  silo  line,  apparent  in  the  remaining  material,  and  by 
breaking  along  the  other  side*  This  difference  in  the 
cons  1 stency  of  debris  could  be  due  to  a  variation  in 
Structure  and  strength  throughout  the  sample*  It  is  also 
apparent  that  the  rupture  is  more  pronounced  on  one  side 
than  on  the  other  side* 

These  features  are  more  visible  on  the  close-up 
presented  on  Plate  6*4  and  on  the  interpretive  sketch  of 
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face 


Figure  6.3  Sample  MC-4  after  Unloading  interpretive  Sketch 

of  Plate  6.4 
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Plate  6*3  Sample  MC-4  after  unloading  lOverali  views 

Sample 


of  the 
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Plate  6*4  Sample  iiC-4  after  unloading  :Close-up  of  the 
tunnel  with  debris;  characteristic  piece  of  debris  collected 

at  the  back  of  the  excavation 


Ill 


6.5  Sample  MC-4  after  unloading  :Close-up  of  the 
tunnel  after  the  debris  has  been  cleaned  out. 


Plate 
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!•  prooapatlon  of  a  logarithmic  spiral  slip  line  through 
the  joints* 

2*  the  shortening  of  these  slip  lines  due  to  a  joint 
3*  tensile  cracks 
4.  comoact  failed  blocks* 

A  small  olece  of  debris  found  at  the  bottom  of  the  failed 
ooening  shows  clearly  two  failure  facestone  of  them  is 
slickensided  and  results  from  a  shear  whereas  the  other  one 
is  irregular  and  may  result  from  tensile  failure*  The  third 
h i dde n  part  of  this  block  is  a  slip  line  parallel  to  the 
joints* 

Plate  6*5  shows  the  tunnel  after  the  debris  has  been 
removed*  A  strain  gauge  is  now  visible*  The  picture  shows 
the  wall  inside  the  tunnel*  the  granular  aspect  of  the 
broken  section*  and  the  intact  sec+ion  in  the  direction  of 
the  joints*  A  crack  developed  at  the  limit  of  this  intact 
section*  It  may  be  due  to  tension  in  the  "roof"  of  the 
elllotlcal  cavity*  The  material  is  highly  jointed  as  is 
apparent  on  the  second  picture  in  Plate  6*5*  Near  the  end  of 
the  excavation*  the  material  is  more  dislocated  than  near 
the  sides*  Slip  lines  start  from  the  acute  end  in  both 
directions*  sometimes  using  joints  and  cross-cleats,  or 
shearing  the  material*  At  the  bottom  of  the  same  picture*  a 
crack  crosses  the  joints*  The  corresponding  block  will  have 
failed  by  splitting  of  this  crack  and  shear  along  a  joint*  A 
similar  mechanism  caused  the  cavity  visible  on  the  tunnel 


wall  at  the  bottom  of  the  first  picture. 
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The  preceedlng  describes  only  one  test*  The 
Interpretation  of  the  failure  and  of  the  role  played  by 
different  cracks  has  to  be  confirmed  by  other  tests*  carried 
up  to  this  state  of  fracture  prpagation*  The  following 
conclusions  can  be  drawn  from  the  observations; 

!•  A  local  instability  generates  the  rupture  which  then 
propagates  in  giving  to  the  cavity  an  elliptical  shape* 
2*  This  propagation  did  not  mean  collapse  in  this  test,  as 
+he  opening  reached  a  stable  condition  at  each  step 
after  propagation  stopped* 

3*  The  preexisting  discontinuities,  mainly  the  joints,  olay 
an  Lmnortant  role  in  the  rupture  process,  by 
intercepting  and  shortening  the  slip  lines  (see  further 
F  i  sure  6*5)* 

4*  Each  failed  lump  is  separated  from  the  mass  along  two 
crossing  slip  lines;  the  displacement  occurs  by  sliding 
on  the  line  subparallel  to  the  joints  and  splitting  on 
the  other  line  (  see  Figure  6*3  )• 

6.*_2*.2*_2  Description  of  the  Rupture  as  Recorded  by  the 
Inst  rumen  t s 

In  this  section  the  measurements  recorded  during  the 
ruoture  process  will  be  presented*  They  will  only  be 
described  and  commented  on  briefly* 

The  stress-strain  curves  of  the  tests  MC— 4*2  and  MC— 4*3 
up  to  failure  are  given  in  Appendix  Bl  and  described  in 
Chapter  4.  The  strain-time  curves  and  the  strain  rate-time 
curves  are  presented  in  Appendix  B2  for  the  last  two  load 


. 
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steps  of  test  MC— 4*2*  Test  MC— 4*3  has  a  similar  behaviour* 
and  will  not  be  presented  In  detail  here*  As  explained  in 
Chanter  4  all  the  experimental  results  are  summarized  in  an 
Internal  Report  of  the  Department  of  Civil  Engineering* 

Appendix  B2  contains  creep  curves  for  the  extensometers 
and  the  corresponding  strain  rate— time  curves  on  a  log— log 
scale  for  the  first  and  the  second  loading  steps*  The  tunnel 
closures  and  the  external  displacements  as  a  function  of 
time  are  also  presented  for  both  load  steps*  In  each  step 
the  curves  ere  numbered  from  9  to  36*  Fach  page  of  Log 
strain  rate  vs*  Log  time  diagrams  presents  four  different 
curve=:  c  o  r  r  e  spond  i  n  g  to  the  four  curves  grouped  on  a  single 
plot  of  the  creep  curves*  Numbers  and  symbols  correspond* 

In  the  first  load  step  presented  (13*3  MPa)  there  is 
evidence  of  failure  on  the  curves  10,  11,  14,  18  and  19(page 
B2*1  of  the  appendix)*  The  corresponding  instruments  are  all 
located  in  on*»  zone  near  the  tunnel  where  failure  seems  to 
occur  first*  As  a  reference,  this  side  of  the  failure, 
corresponding  to  the  lower  right  corner  on  the  sketch 
representing  the  sample,  is  the  one  where  failure  is  of 
lesser  importance  in  the  upper  picture  on  Plate  6*3  (right 
side  of  the  picture )• 

Large  deformations  occur  for  gauge  14  In  extension  and 
7?  in  compression  (page  B2*l)  •  This  may  be  explained  by 
stress  redistribution  due  to  rupture  of  the  wall.  The 
response  of  the  curves  18  and  19  Is  delayed,  corresponding 
to  a  further  propagation  of  the  ruptured  zone* 
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The  strain  rate— time  curves  give  an  accurate  picture  of 
the  rupture*  Unlike  the  creep  curve,  they  Indicate  that  all 
the  instruments  respond  to  the  rupture*  Curve  14  (page  B2*3) 
Is  typical  for  the  rupture  orocess*  It  shows  successive 
acceleration  of  strain  uo  to  a  peak,  corresponding  to 
different  rupture  stages*  Instruments  12  and  16  should  show 
the  fracture  visually  observed  in  this  region,  but  their 
resoonse  is  not  obvious*  The  first  sign  of  rupture  occurred 
on  curve  10  (page  B2*2)  by  an  acceleration  in  compression 
about  S  minutes  after  the  time  zero  of  the  creep  test,  and 
it  Is  followed  by  the  extension  of  the  instrument  14  (page 
B?*3)»  3  minutes  after*  The  other  instruments  then  responded 
to  a  variable  extent* 

The  creep  curves  for  the  next  load  step  are  different 
( page  B2*7),  They  are  characterized  by  extremely  large 
deformations  *  and  the  c or r esoo nd l ng  scale  does  not  allow 
the  analysis  of  smaller  deformations*  The  strain  rate-time 
plots  are  once  more  useful*  Most  of  the  curves  present  a 
plat*3 au  at  the  beginning,  corresponding  to  a  higher  value  of 
the  strain  rate  than  in  the  previous  load  steo.  The  rapid 
succession  of  ruptures  cause  the  wave-like  aspect  of  the 
plots  to  be  suppressed;  only  peak  values  remain*  The 
instruments  monitoring  the  behaviour  of  the  "upper  left" 
side  of  the  tunnel  respond  now  and  the  measurements 
correspond  to  the  large  deformations  visually  observed*  The 
curves  13  to  16  (page  B2.9)  are  typical;  the  Instruments  14 


an 


d  16  in  the  failure  zone  show  a  constant  high  value  of  the 
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strain  rate,  and  Instruments  13  and  15  show  a  constant 
decrease  In  extension  rate*  If  the  long  axis  of  the  Intact 
elliptical  cavity  Increases  regularly,  the  influence  on  the 
elastic  closure  in  the  direction  of  the  small  axis  decreases 
regularly* 

The  closures  of  the  tunnel  also  reflect  the  rupture 
(page  B2*13)«  They  are,  however,  difficult  to  Interpret  as 
the  fractured  material  disturbs  the  readings  of  the 
Instruments.  The  external  displacements  are  also  affected  by 
the  rupture  of  the  tunnel  wall*  Large  strain  rates  required 
a  change  in  the  strain-axis  scale  which  has  to  be  noted*  The 
discontinuities  in  strain  visible  on  some  curves  (page 
B2*13)  result  from  technical  problems  generated  by  the  large 
de  formations* 


6 . 3  Interpretation  o f  the  Rupture  Prfic ess. 

The  purpose  of  this  section  is  to  explain  the  rviptur 
mechanism  in  the  model  test*  More  tests  are  required  t 
confirm  the  hypotheses  presented  here,  and  to  apply  them  t 
real  openings*  However,  the  recognition  of  some  phenomena  i 
rplevant  as  it  may  help  in  the  interpretation  of  othe 
rupture  cases* 
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6.3.1  Effect  of  Preexisting  Discontinuities  on  the  Slip 
Li  nes 

The  theory  of  plastic  equilibrium  shows  that  the  slip 
lines  around  a  falling  circular  opening  hove  the  shape  of 
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logarithmic  spirals  for  a  Coulomb  material*  In  other  words, 
they  make  everywhere  within  the  plastic  zone*  the  same  angle 
with  the  principal  directions  of  stress*  This  angle  is  the 
friction  angle  of  the  mat er 1 al(  used  in  the  Coulomb  failure 
criterion)  with  respect  to  the  maximum  principal  stress 
direction*  Jaeger(1960)  was  one  of  the  first  to  explain  in 
de+all  the  anisotropic  character  of  the  strength  in  a 
jointed  rock  mass*  Patton(1966)  presented  the  concept  of 
multiple  mode  of  shear  failure  in  rock  :  Depending  on  the 
relative  orientation  of  the  joints  with  respect  to  the 
direction  of  principal  stress,  the  rock  mass  is  sheared 
either  along  a  joint  or  through  the  intact  material* 

'f'he  above  phenomenon  has  been  observed  in  the  sample* 
The  slip  lines  shear  the  material  along  a  joint  if  the  angle 
of  Incidence  is  large  enough*  Otherwise  they  follow  the 
encountered  joint*  A  joint  constitutes  a  "shortcut"  for  a 
slip  1 i ne • 

A  qualitative  interpretation  of  the  influence  of  joints 
on  +he  slip  lines  around  a  tunnel  is  presented  in  Figures 
6.4  and  6*5*  This  method  was  already  suggested  by 
Rechsteiner  and  Lombar d i(  1 97 4 )  .  On  a  Mohr  circle  diagram 
three  different  circles  are  represented,  corresponding  to 
the  state  of  stress  at  three  different  points  along  a  slip 
line*  As  the  slip  line  is  within  the  plastic  zone,  the 
circles  are  tangent  to  the  Intact  rock  failure  envelope, 
corresponding  to  the  failure  criterion*  Because  of  the 
assumed  axl syme t ri cal  conditions  all  the  points  located  in 


* 
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th**  reck  mass  on  a  circle  concentric  with  the  tunnel  have 
the  same  stress  level*  Hence  each  Mohr  circle  on  Figure  6*3 
represents  the  state  of  stress  for  every  point  on  one  of 
these  concentric  circles*  The  set  of  three  Mohr  circles  are 
then  available  for  the  study  of  any  slip  line  around  the 
tunnel • 

Circle  1  has  a  minimum  principal  stress  equal  to  zero* 
1+  represents  the  state  of  stress  on  the  tunnel  wall*  Mohr 
circles  2  and  3  correspond  to  particular  circles  concentric 
with  the  tunnel*  For  the  points  on  a  particular  logarithmic 
spiral  slip  line,  the  principal  directions  of  stress  will 
rotate  with  Increasing  the  radial  distance  of  the  point* 
This  particular  slip  line  is  intercepted  by  circles  1,2  and 
3  at  points  Al,A2  and  A3*  Circle  2  will  be  characterized  by 
the  rotation  in  the  principal  directions  between  Al  and  A2 , 
circle  3  between  Al  and  A3*  These  angles  are  20°  and  40°  for 
the  circles  2  and  3  respectively*  The  relevance  of  such  an 
unusual  characterization  is  that  it  makes  the  determination 
of  the  traces  of  a  discontinuity  on  the  Mohr  circle  easier. 

Four  couples  of  slip  lines  are  chosen  arbitrarily  for 
this  study;  points  A,  B,  C  and  D  are  the  intersections  of 
these  slip  lines  with  the  tunnel  wall  (Figure  6*5).  The  Mohr 
diagram  Is  symmetrical  with  respect  to  the  abscissa  axis* 
From  each  of  the  four  points  a  pair  of  slip  lines  are  then 
Initiated  within  the  rock  mass,  symmetrical  with  respect  to 
the  radial  direction*  Along  one  slip  line  the  shear  stress 


is  positive  and  the  failure  corresponds 


to  the  envelope  P  in 
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Figure  6.4  Influence  of  Preexisting  Discontinuities  on 

Plastic  Slip  Lines  (1) 


' 


120 


Figure  6.5  Influence  of  Preexisting  Discontinuities  on 

Plastic  Slip  Lines  (2) 
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Figure  6  •  4  •  Along  the  other  slip  line  the  shear  stress  is 
negative  and  the  failure  corresponds  to  the  envelope  N  in 
Figure  6.4* 

Up  to  this  point,  the  anisotropic  character  of  the 
material  has  not  been  considered*  A  double  linear  envelope 
fo~  the  iilure  along  a  plane  of  weakness  is  given  in  Figure 
6*4#  ■  >  e  trace  of  the  joints  will  be  plotted  in  the  Mohr 
circles  as  indicated  on  the  sketch  in  the  left  corner  of 
F i sure  6*4*  If  the  trace  of  t he  joints  is  above  the  joint 
envelooe  the  material  will  be  sheared  along  a  joint*  The 
joints  make  respectively  an  angle  of  0,  30,  60  and  90 
degrees  with  the  direction  of  the  maximum  principal  stress 
(  tangential  stress)*  If  on  a  elven  Mohr  circle  (i*  e*  circle 

1  )  the  traces  of  the  joints  are  plotted,  they  will 
correspond  to  the  points  a,  b,  c  and  d  as  Indicated  on  the 
circle  1  in  Figure  6*4*  Actually,  the  principal  directions 
a~e  rotating  giving  four  superposable  Mohr  circles  which  are 
represented  by  a  unique  circle  on  the  figure* 

These  conventions  being  established,  the  traces  of  the 
joints  on  the  circles  2  and  3  will  be  plotted  in  the 
following  manner*  To  move  along  a  slip  line  from  the  tunnel 
wall  to  the  point  where  the  state  of  stress  is  given  by 
circle  2,  corresponds  to  a  rotation  of  ±20°  of  the  principal 
directions  of  stress*  The  traces  of  the  joints  on  the  circle 

2  are  given  by  a  rotation  of  ±20°  on  this  Mohr  circle  with 
respect  to  the  former  radial  position  on  circle  1*  The  same 
procedure  is  used  on  circle  3* 


* 
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Initiated  at  each  point  a*  b,  c  and  d  a  double  path 
corresponding  to  a  particular  couple  of  slip  lines  may  be 
established  in  such  a  manner*  The  dashed  lines  correspond  to 
a  positive  rotation  along  the  Mohr  circle  and  to  a  failure 
with  respect  to  the  envelope  P,  the  dotted  lines  to  a 
negative  rotation  and  to  a  failure  with  respect  to  the 
envelope  N*  If  a  point  of  the  trajectory  Is  over  the  joint 
envelope*  the  corresponding  slip  line  will  shear  the  rock 
along  the  joint*  Otherwise*  the  line  will  shear  across  the 
intact  material*  Corresponding  slip  lines  initiated  on  the 
tunnel  wall  at  A *  B*  C  and  D  are  shown  in  Figure  6*3* 

The  anisotropic  redistribution  of  stress  resulting  from 
the  failure  along  a  joint  has  not  been  taken  into  account*  A 
more  sophisticated  anisotropic  plastic  epui librium  analysis 
will  be  necessary  to  explain  the  complete  phenomenon* 

The  actual  slio  line  net  is  affected  by  the 
n on  — ax  1 symme tr ic  characteristics  of  the  material*  Only  lines 
near  A  or  3  will  develop*  This  simple  analysis  neglects  the 
effect  of  the  other  discontinuities  like  the  cross-cleat 
joints  and  existing  faults  which  have  been  proven  to  affect 
the  mechanism* 

The  recognition  of  the  Influence  of  the  joints  is, 
however*  Important  near  point  A*  as  i+  is  one  of  the  reasons 
for  a  rupture  mechanism  to  be  initiated* 


6*3*2  Hypotheses  for  the  Rupture  Development 
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6  * 3 «  2 * 1  Introduction 

Test  MC— 3  did  not  show  any  rupture  but  test  MC— 4  was 
conducted  beyond  the  rupture  of  the  tunnel  wall*  The  nature 
and  the  properties  of  the  material  are  similar  for  both 
tests  and  the  stress  history  Is  comparable*  It  Is  assumed 
that  the  difference  in  structure  between  the  samples 


descr!  bed 

in  the 

section  2*4  will  not 

in  terfere 

wi  th 

the 

rupture 

process* 

The  only  difference 

is  a  chang 

e  in 

the 

tunnel 

di  a  meter* 

It  has  been  shown 

in  sect 

ion 

5.3 

(  " Influence  of  the  Boundaries")  that  this  diameter  change 
does  not  Induce  a  significant  change  in  the  stress  level 
around  the  tunnel*  A  different  shape  of  the  viscoelastic 
stress  r ed 1 s t r ibut i on  could  also  Influence  the  stresses  but 
it  Is  assumed  that  it  is  not  significant  enough  to  cause 
this  difference  in  ruoture  behaviour*  Actually  sample  MC— 3 
was  finally  loaded  to  the  maximum  capacity  of  the 
compression  machine  with  a  high  loading  rate(  this  should 
partly  overcome  the  eventual  problem  of  the  influence  of  the 
stress  red  is tr ibut i on  ) »  but  the  tunnel  wall  did  not  rupture* 
I ♦  seems  that  the  rupture  mechanism  is  influenced  by  the 
local  radius  of  curvature* 

6  *  3  *  2  *  2  Rupture  initiation 

Three  types  of  rupture  initiation  may  be  recognized* 
They  are  shown  on  the  Figure  6*6  B  and  are  as  follows  J 
1*  Splitting  :  A  column  of  material  is  isolated  by  two  slip 
lines  and  a  joint  parallel  to  the  tangent  to  the  tunnel 


. 


■ 
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Possible  Shape  For  a  Slip  Lines  Net  Observed  Fractures 


Buckling  Sliding 


Figure  6.6  A  :  Possible  Slip  Lines  around  an  Elliptical 
Cavity  and  Observed  Fractures* 

B  :  Different  Types  of  Initiation  for  the  Rupture  of  the 

Tunn  e l  Well 


125 


surface*  It  Is  subjected  to  compression  by  the  maxi 


mum 


principal  stress  and  fails  by  buckling.  The  width  of  the 
column  and  the  slenderness  ratio  are  Influenced  for 
geometrical  reasons  by  the  local  radius  of  curvature* 

2*  sliding  :  A  pair  of  slip  lines  are  initiated  at  the 
point  A  of  Figure  6*5*  They  isolate  two  symetrical 
blocks  which  are  subjected  to  the  same  compression  as 
the  block  in  the  previous  description*  The  blocks  slide 
on  one  side  and  are  split  along  their  bases  on  the 
second  set  of  slip  lines. 


3*  Kirschkern  C "Cherrystone") 


This  phenomenon  is 


described  by  Feder  (  1978)*  It  is  the  expulsion  of  a 
block  out  of  the  wall  by  sliding  on  both  sides* 

These  effects  result  from  local  instabilities*  They  are 
three-dimensional  and  depend  on  the  degrees  of  freedom  in 
the  longitudinal  direction*  It  is  difficult  to  make  a 
distinction  between  the  three  mechanisms  by  visual 
observations*  The  arrow— like  aspect  of  the  observed  rupture 
may  indicate  either  buckling  or  sliding*  The  important  fact 
is  that  the  three  mechanisms  are  dependent  on  the  radius  of 
curvature  of  the  opening*  A  smaller  radius  will  decrease  the 
area  involved  in  any  of  the  rupture  process  and  cause  the 


"eooarent  strength"  of  the  tunnel  wall  to  be  higher* 


6*3*2. 3  Rupture  Development 

A  rupture  initiation  gives  a  preferential  direction  for 
the  progression  of  this  phenomenon*  The  opening  takes  an 
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ellip+ical  shape.  The  stress  concentration  increases  at  the 
ends  of  the  growing  long  axis  of  the  ellipsPf  and  rupture 
proceeds  by  one  of  the  described  mechanisme.  The  blocks  tend 
to  slide  and  rotate  towards  the  long  axis  of  the  ellipse. 

Figure  6.6A  shows  on  an  elliptical  cavity  a  possible 
slip  line  net  without  the  influence  of  discontinuities  and* 
on  the  other  half,  the  main  features  observed  on  the  actual 
opening.  As  rupture  proceeds,  new  slip  lines  corresponding 
to  the  elliptical  opening  are  generated.  The  shearing 
process  on  +he  sides  of  the  opening  is  complex  as  it 
'nvolvps  these  new  lines  as  well  as  the  lines  generated  when 
the  opening  was  circular.  This  particularity  is  illustrated 
on  F i gur e  6.3. 

As  the  ellipse  propagates,  and  if  the  stress  stays 
constant,  the  tunnel  would  collapse  immediately  and 
completely  if  the  strength  was  constant.  However,  the 
rupture  process  was  slower  after  some  hours  of  creep  in  test 
MC-4.2  and  stable  in  test  MC-4.3  (reloading  of  the  broken 
tunnel).  In  order  to  achieve  that  ptabil izetion,  the 
aoparent  strength  must  increase.  Actually  the  radius  of 
curvature  tends  to  decrease  continuously  if  the  long  axis  of 
the  elliptical  cav? +y  grows.  At  the  same  time,  and  because 
of  this  decrease  in  the  radius  of  curvature,  the  rupture  of 
the  wall  becomes  more  difficult.  This  could  be  generalized 
in  saying  that  the  apparent  strength  of  the  wall  increases. 

The  rupture  mechanism  has  been  recognized  to  be 
curve  ture~deppndent  and  this  relates  to  the  stabilization  of 
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the  phenomenon* 
association  with 
jointed  material* 


This  concept  will 
the  size  effect  on 


be  visualized  in 
the  strength  of  a 


6*.4  Shape  Effect  and  the  St abi  11  ty  of  &  Tunnel  Wall 

6*4*1  Introduction 

The  rupture  initiation  and  propagation  have  been 
described  in  the  previous  section*  The  actual  physical 
orocesses  were  described*  It  is  of  interest  to  Introduce  a 
macro  -  analysis  of  this  rupture  phenomenon*  For  this 
purpose  the  "apparent  strength”  is  introduced  as  the  stress 
required  to  Initiate  a  rupture  on  the  intact  or  already 
damaged  tunnel  wall*  In  this  section  the  change  in  apparent 
strength  with  propagating  the  rupturing  opening,  will  be 
discussed  in  order  to  study  the  stabilization  of  the  rupture 


pr oc  ess* 


£*,4.2.2  Shape  Effect  on  the  Strength  of  the  Tunnel  Walt 

The  curvature  of  the  opening,  or  in  other  words  the 
shape  of  the  element  involved  in  a  failure  process,  has  been 
shown  to  affect  the  rupture  mechanism*  The  element  of  the 
tunnel  wall  involved  in  a  rupture  will  be  assumed  to  be 
limited  by  an  elementary  solid  angle  da,  its  width  being 
constant  as  a  function  of  the  joint  spacing.  This  is 
Illustrated  on  the  upper  sketch  in  Figure  6.7*  If  the  radius 
of  curvature  is  larger,  the  element  will  be  longer*  Hudson 
at  aJL (  1972)  have  shown  the  Influence  of  the  shape  of  the 
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Size  of  a  Failing  Element 

hi=ad0  :  circular  opening 
h 2  =  a2d0/b  elliptical  " 
normalized  size  : 

0<  h2/hi  =  a/b  <  1 


strength  vs.  Size  of 
the  failing  element 

Intrinsic  curve  for 

the  material 


depend  on  : 


-  Size  of  the  original 
tunnel 

-  Size-dependent  Strength 


Curves  2  : 

V2  \  a„  =  2a„ib 

depend  on  : 

-  Opening  Shape 

a/b 

-  Stress  Conditions 

a0 


Figure  6.7  Illustration  for  the  Model  of  Stabilization  of 

the  Rupture 


. 
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tpcted  element  on  the  complete  stress  strain  curve  of  a 
rock*  The  results  obtained  for  marble  samples  in  uniaxial 
compression  are  reproduced  in  Figure  6*8  and  are  given  for 
Illustration  only*  It  can  be  seen  on  these  plots  that  the 
shape  of  the  rock  element  affects  not  only  the  peak  strength 
hut  also  the  post— peak  behaviour*  A  longer  element  will  have 
lower  strength  and  a  steeoer  post— peak  strength  curve* 

A  similar  behaviour  may  be  expected  for  a  rock  element 
on  a  tunnel  wall  even  if  the  boundary  conditions  are 
dl f f eren t  «  This  will  be  used  to  qualitatively  analyse  the 
rupture  process* 

6*4 *3  Ti me— Dependent  Rupture  Process 

After  the  sample  has  been  loaded  near  rupture,  the 
elements  near  the  tunnel  wall  experience  a  strength  failure 
either  immediately  or  after  some  creep*  Rupture  subsequently 
occurs,  creating  a  new  elliptical  opening*  Elements 
initially  further  from  the  wall  are  suddenly  at  the  surface 
of  the  new  opening*  As  the  stress  concentration  increases, 
they  support  more  load.  The  radius  of  curvature  decreases, 
the  apparent  strength  increases,  end  the  new  cavity  may 
supoort  the  load*  If  the  stress  increases  more  rapidly  than 
the  strength,  the  new  opening  will  collapse*  This  process 
goes  on  until  the  curvature  of  the  opening  provides  a  higher 
strength  than  the  stress  level  and  the  rupture  is 
s  t abl 1 i zed • 

This  repetitive  process  explains  the  wave-like  aspect 


■ 
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SIZE  EFFECT 


SHAPE  EFFECT 


FI gure 
Strain 


6*8  Size 
Curve  of 


and  Shape  Effects  on  the  Complete  Stress 
Marble  Samples  in  Uniaxial  Loading  (after 
Hudson  sJL  ILU  1972  ) 
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of  the  Instrument  response  described  In  the  second  section 
of  this  chapter# 

6.i4#.4  Model  for  the  Stabilization  of  the  Rupture  Process 

The  stabilization  of  the  rupture  process  depends  on  the 
relative  Importance  of  the  stress  Increase  due  a  change  in 
opening  shape  and  the  strength  increase  due  to  a  change  in 
element  shape#  This  has  been  recognized  to  a  certain  extent 
by  Falrhurst  and  Cook(  1966) 

The  model  description  Is  Illustrated  on  Figure  6#7#  The 
height  h  of  a  rupturing  element  is  equal  to  the  local  radius 
of  curvature  mullplled  by  the  elementary  solid  angle  da#  The 
radius  of  curvature  at  the  end  of  the  long  axis  of  an 
ellipse  Is: 

r=  a  2  /  b 

where  a  and  b  are  resoectlvely  half  the  lengths  of  the  small 
and  the  long  axis#  The  normalized  size  of  the  element, 
defined  as  the  ratio  between  the  length  of  the  actual 
rupturing  element  to  the  initial  rupturing  element  length 
(circular  cavity)  is  equal  to  the  ratio  of  the  ellipse  axes: 

h/ada=a2da/b  a  da  =a/b 

If  the  behaviour  of  the  materiel  is  time-independent,  the 
tangential  stress  at  the  wall  in  the  direction  of  the  long 
axis  is  : 

6=262#b/ a 

if  69  is  the  external  stress# 

Then,  for  the  opening  to  be  stable,  the  unconfined 


< 
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compressive  strength  has  to  be  at  least  equal  to  this  value 
of  the  tangential  stress  which  will  be  called  here  11  required 
strength”*  A  relation  is  found  between  the  required  strength 
for  an  opening  to  be  stable  and  the  normalized  size  of  the 
rupturing  element.  This  curve  is  called  curve  2  in  Figure 
6*7.  Figure  6*9  shows  the  corresponding  family  of  hyperbolas 
for  different  values  of  the  external  stress  (solid  lines* 
c  urves 2  )  • 


These  curves  a 
rock  mass  and  show  t 
the  elliptical  cavlt 
The  Influence 
its  strength  has 
Fienlawski( 1968)  or 
this  variation*  Th 
strength  of  the  s 
assumed  thereafter  t 
size  h  of  a  rupt 
strength  of  the  mate 
for  the  material* 
function  of  the  s 
the  opening  and  as 
curve  expressing  the 
of  the  normalized 
strength  curve  by  an 
axis  and  whose  rat 
on  Figure  6*7*  A  p 


re  a  function  of  the  stress  level  in  the 
he  evolution  of  the  required  strength  as 
y  propagates* 

of  the  shape  of  tested  coal  samples  on 
been  shown  by  several  authors* 
Hus tru l i d(  1975  )  gave  relationships  for 
ey  found  a  unique  curve  relating  the 
ample  to  the  size  of  the  element*  It  is 
hat  a  similar  curve  exists  to  relate  the 
uring  element  on  the  tunnel  wall  to  the 
rial*  This  curve  is  unique  and  intrinsic 
Given  this  curve  of  the  strength  as  a 
lze  of  the  sample*  given  the  radius  a  of 
turning  a  value  for  the  solid  angle  da*  a 
variation  of  the  strength  as  a  function 
size  may  be  obtained  from  the  intrinsic 
axial  affinity  parallel  to  the  abscissa 
lo  is  1/adar?  This  process  is  Illustrated 
oint  of  Curve  1  is  obtained  from  a  point 
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Figure  6*9  Model  for  the  Stabilization  of  the  Rupture 

Process  in  the  Model  Test 
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of  the  Intrinsic  curve  by  keeping  the  ordinate  constant  and 
multiplying  the  abscissa  by  the  ratio  1/ador#  The  curves  1 
depend  on  the  material  properties  and  on  the  size  of  the 
opening#  Different  opening  sizes  define  another  family  of 
curves  in  the  strength  versus  normalized  size  axes#  These 
curves  show  the  variation  of  the  provided  strength  when 
rupture  proceeds# 

The  stress  ,  material  and  geometrical  conditions  will 
give  a  curve  of  each  of  the  two  families.  The  comparison  of 
■♦•he  provided  strength  and  the  required  strength  will  give  an 
idea  of  the  stability  of  the  tunnel  wall# 

6#  4#  ?  Experimental  Data  on  the  Rupture 

From  the  load  and  rupture  history  of  sample  MC-4,  and 
from  the  measurements  of  the  tunnel  section  after  rupture, 
it  is  possible  to  get  experimental  points  In  the  previously 
described  axes#  These  points  are  on  the  experimental  curve  1 
corresponding  to  a  15#2  cm#  tunnel  in  Sundance  coal#  This 
has  been  done  in  Figure  6#9  (  dashed  line)#  Failure  initiates 
for  an  external  stress  of  13  MPa; the  sample  was  then  loaded 
up  to  16  MPa,  unloaded  and  reloeded#  The  co r re soond i ng 
values  of  the  axes  ratio  give  the  abscissa  (the  normalized 
size  has  been  shown  to  be  equal  to  the  axes  ratio  of  the 
ellipse)#  When  the  sample  is  reloaded  after  the  debris  has 
been  removed,  the  new  rupture  initiates  for  values  of  the 
external  stress  between  11  MPa  and  13  MPa#  This  shows  that 


the 


debris  provides  a  small  support  pressure  which  is  not 
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taken  into  account  here,  but  seems  to  have  an  Important 
e  f  fee  t « 

From  this  experimental  curve  an  intrinsic  curve  for  the 
strength  of  the  material-  as  a  function  of  size  may  be 
obtained#  Furthermore  a  back  calculated  curve  1  may  be 
obtained  for  sample  MC— 3*  An  horizontal  affinity  with  ratio 
a^dar  will  transform  the  experimental  curve  for  test  MC—  4  to 
the  intrinsic  strength  vs*  size  curve  for  Sundance  coal t  if 
a*  is  the  radius  of  the  MC— 4  tunnel*  The  curve  1  for  test 
MC— 3  results  from  a  further  transformation  of  the  intrinsic 
curve  by  an  horizontal  affinity  with  ratio  l/a3da,  where  a3 
is  the  radius  of  the  MC— 3  tunnel*  In  other  words  the  curve  1 
for  sample  MC— 3  may  be  calculated  from  the  curve  1  for 

■=ample  MC— 4  by  keeolng  at  each  point  the  ordinate  constant 
and  by  multiplying  the  abscissa  by  the  ratio  a*/a3*  The 
results  of  this  transformation  are  presented  on  Figure  6*9 
as  a  dotted  line*  According  to  this  curve  it  is  normal  that 
the  sample  MC-3  did  not  rupture,  but  it  is  also  reasonable 
to  assume  that  plastic  yielding  already  accurred,  in 

agreement  with  the  plastic  lines  observed  when  the  sample 
was  unloaded  from  l7MPa.  The  intrinsic  curve  relating  for 
Sundance  coal  strength  and  size  of  the  rupturing  element  is 
not  presented  here  but  it  is  in  agreement  with  an 
extrapolation  of  published  results  to  small  size,  even  if 
these  published  curves 
sa  mpl e • 


relate  strength  and  size  of  a  free 
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6 > 4 ♦ 6  Discussion  on  this  Model 

This  model  explains  the  behaviour  of  the  model  test 
samples  with  respect  to  the  rupture*  It  is  based  on 
slmolifled  hypotheses  and  needs  to  be  extended  to  Include 
for  instance  time-dependency  or  s t ra i n— s o f ten i ng  effects  in 
order  to  be  applicable  to  the  reality*  However,  some  basic 
conclusions  may  be  drawn* 

1*  For  the  same  material,  a  smaller  tunnel  will  have  a 
strength  versus  normalized  size  curve  higher  than  that 
for  a  larger  tunnel*  Thus  for  a  given  density  of  joints 
a  small  tunnel  is  more  stable  than  a  larger  one*  This 
corresponds  to  real  observations  described,  for 

Instance,  by  Talobrei  1957  )  (  p*  318)* 

2*  The  stability  of  an  opening  will  depend  on  the  relative 
slopes  of  the  provided  strength  (Curvel  )  and  required 
strength  (Curve  2)  curves*  The  lower  grach  on  Figure  6*7 
illustrates  that  if  curve  1  tends  to  be  above  curve  2 
the  tunnel  is  stable,  otherwise  the  rupture  will 


propaga 


te*  A  material  for  which  strength  does  not  depend 


on  element  size  will  exhibit  collapse  and  a  complete 
rupture  as  the  provided  strength  curve  is  flat  and  the 
required  strength  curve  remains  unchanged* 

3.  The  same  conclusion  is  valid  for  a  large  opening  in  a 
dense  material  as  the  provided  strength  curve  is  also 
relatively  flat* 

4*  The  rupture  of  a  tunnel  is  seal e— dependen t  and  a  study 
of  the  laws  of  similitude  applicable  to  this  case  will 


* 
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be  necessary  to  show  under  which  conditions  the  same 
type  of  controlled  failure  may  occur  in  reality* 

6jL.fi  Sumifl£2  and  Conclusions 

The  following  conclusions  may  be  drawn  from  the 
observations  and  analysis  of  the  rupture  process  ? 

1*  The  observation  of  the  nlaster  covering  the  sample 
Indicates  that  the  stress  distribution  in  the  model  test 
is  radial  around  the  opening*  An  elliptical  crack 
pattern  observed  on  one  of  the  samples  indicates  an 
anisotropic  stress  redistribution  • 

2*  Hven  if  the  tunnel  wall  does  not  rupture,  a  small 
olastic  zone  develops  around  the  tunnel,  allowing  some 
blocks  to  be  separated  from  the  rock  mass  ; this  rock 
fall  phenomenon  is  Influenced  by  preexisting  faults  on 
+he  tunnel  wall* 

3*  Rupture  of  the  tunnel  wall  is  triggered  by  a  mechanism 
involving  the  joints  of  the  rock  mass*  It  initiates  in 
the  direction  perpendicular  to  the  joints,  and  is  caused 
either  by  a  sliding  or  a  buckling  of  a  “band"  of  coal 
(between  two  joints)  towards  the  cavity* 

4.  The  joints  Influence  the  rupture  mechanism  by 
intercepting  plastic  slip  lines  end  in  this  way 
shortening  the  length  of  their  resistant  section* 

5*  The  rupture  mechanism  depends  on  the  radius  of  curvature 
of  the  surface  to  be  broken*  k  surface  with  a  smaller 
radius  of  curvature  has  a  higher  resistance  ©gainst  such 


a  rupture  mechanism* 
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6*  The  rupture  propagates  In  the  direction  perpendicular  to 
the  joints*  creating  a  new  elliptical  opening* 

7*  The  rupture  is  better  recorded  by  strain  gauges  inside 
the  rock  mass  than  by  closure  measurements*  The 
t l me— dependency  of  the  rupture  is  well  recorded  on  a  log 
strain  rate  versus  log  time  representation  of  the 
internal  strain  gauges  measurements*  These  curves  have  a 


characteristic  wave— like  aspect* 

8*  The  rupture  process  stabilizes  in  the  model  test*  A 
sample  model  taking  into  account  the  shape-dependent 
strength  of  a  rupturing  element  has  been  developed*  It 
compares  the  strength  provided  +o  an  element  on  the 
tunnel  wall  by  its  shape  with  the  strength  required  to 
keep  an  opening  stable*  if  this  opening*  Involved  in  a 
rupture  process*  changes  its  shape  from  circular  to 
elllp+ical* 

<3  •  This  model  explains  the  stabilization  of  the  rupture  in 
t he  test  MC— 4  *  a  nd  why  the  sample  MC— 3  stayed  intact*  It 
shows  also  that  for  a  given  jointed  material,  a  small 
tunnel  is  more  stable  than  a  larger  tunnel* 

It  appears  from  these  observations  that  an  elastic  — 
anisotropic  plastic  analysis  should  be  undertaken  to  analyse 
the  stability  of  a  tunnel  in  a  jointed  brittle  rock  mass. 
The  importance  of  the  shape  of  the  excavation  with  respect 
to  +he  orientation  and  the  density  of  the  joints  has  to  be 
recognized  in  order  to  control  a  local  rupture.  A  local 
instability  could  initiate  a  complete  collapse  of  the 
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CHAPTER  7 


CONCLUSIONS*  IMPLICATIONS  AND  RECOMMENDATIONS 

7  . \  Conclusion? 

An  extensive  testing  program  of  eight  months  duration 
on  two  different  model  test  samples  was  conducted.  Rupture 
of  the  tunnel  was  obtained  only  in  the  second  sample  and  the 
different  ultimate  behaviour  oermlts  a  comparison  between 
the  two  samples. 

Material ,  testing  equipment  and  procedure  are 
sufficiently  developed  to  study  the  time  -  dependent  failure 
process  of  a  model  test  of  a  tunnel  in  a  natural  material. 
The  testing  technology  is  now  improved  and  will  allow  more 
rapid  collection  of  relevant  additional  data  in  future 
tests. 

In  the  pre-failure  range  the  elastic  and  visco-elastic 
deformations  of  the  sample  are  In  agreement  with  the  results 
of  previous  experiments  on  model  test  and  triaxial  samples. 
Because  of  the  ultimate  behaviour  of  the  second  sample  the 
Interpretation  of  the  test  data  has  been  orientated  towards 
+he  analysis  of  the  near  -  failure  behaviour  of  the  tunnel. 
An  isotropic  e las t o— p la s t i c  analysis  of  the  model 
deformations  has  been  conducted  with  different  hypotheses  of 
post  -  failure  behaviour.  The  experimental  data  are 
consistent  with  the  results  of  this  theoretical  analysis; 
the  corresponding  values  of  the  strength  parameters  are 
reasonable  for  coal  and  the  volume  change  and  post  failure 
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agree  with  the  generally  accepted  practice  in  the 
a 1 1  on  of  field  data* 

vet*  *  the  conclusions  obtained  from  this  analysis 
isual  observations  of  the  failure  mechanism  are  not 
e*  Direct  observations  indicate  that  this  isotropic 
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11. 

This 

occurr  ed 

in  zones  where  the  joints  are 

to  the  tunnel  wall  and  it  was  influenced  by  the 
of  jointing  in  the  material  and  by  the  local  radius 
ture  of  the  cavi ty*  Initial  rupture  may  lead  to  a 
collapse  of  the  tunnel  or  to  a  stabilizing  process* 
model  explains  that  this  ultimate  behaviour  is 
d  by  the  structure  —  dependent  character  of  the 
of  a  brittle  rock*  For  the  same  jointed  material  a 
nel  will  be  stable,  a  larger  tunnel  will  experience 
zing  local  failure  and  a  larger  tunnel  a  collapsing 
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Even  if  the  laws  of  similitude  in  the  inelastic  range 
do  not  allow  at  the  present  time  a  direct  application  of  the 
model  test  results  to  reality,  this  thesis  has  some 
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implications  In  the  practice  of  tunnel  design* 

For  a  tunnel  in  brittle  jointed  material  the 
interpretation  of  field  data  should  no  longer  be  undertaken 
by  use  of  isotropic  elasto-plastic  analysis  only*  It  is 
advised  that  analytical  models  be  developed  taking  into 
account  the  anisotropic  character  of  the  strength,  the 
t 1  me— depen dent  and  time  —  Independent  stress  redistribution 
processes  and  the  implications  of  local  instability* 

The  evaluation  of  rock  characteristics  as  well  as  the 
stress  conditions  Inferred  from  a  pilot  tunnel  could  be 
misleading  if  stress  red  is tr ibut ion  and  scale  effects  are 
not  taken  into  account* 

The  continuing  research  program  should  take  advantage 
of  the  experience  gained  in  the  previous  tests  In  order  to 
conduct  further  experiments  with  maximum  efficiency. 

The  detrimental  effect  of  excessive  drying  -  wetting 
cycles  on  the  structure  and  the  strength  of  a  sample  has 
been  recognized.  This  practice  is  to  be  avoided  by  always 
protecting  the  sample.  It  was  assumed  that  the  rupture 
process  itself  was  not  controlled  by  the  change  in  structure 
of  sample  MC-4  resulting  from  these  effects.  A  complete 
exoeri mental  Investigation  of  the  effect  of  successive 
changes  in  moisture  content  could  not  be  undertaken  within 
the  time  frame  of  this  thesis.  It  should  be  completed  If 
further  experiments  on  model  tests  prove  the  Importance  of 
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this  effec t. 

Similarity  between  the  model  test  and  the  reality  would 
be  improved  if  the  tunnel  could  be  drilled  through  a  sample 
under  load.  However,  under  these  conditions  the  tunnel 
closure  will  only  be  measured  after  the  time  delay  necessary 
to  Install  the  measuring  device. 

Time  —  deoendent  stress  redistribution  has  been  proven 
to  be  essential  in  the  analysis  of  tunnel  deformations.  It 
is  difficult  (  if  not  impossible)  to  monitor  this  local 
effect  i ndependen t ly  In  the  model  test.  A  separate 
experiment  could  be  undertaken  to  study  the  rate  and  the 
amplitude  of  the  stress  redistribution  process.  This  might 
be  achieved  by  pushing  an  element  of  coal,  cut  by  two 
longitudinal  sections  from  a  hollow  cylinder  of  coal, 
against  a  rigid  wedge— shaped  frame  instrumented  with  load 
cells.  The  major  oroblem  fo r  the  feasibility  of  this 
experiment  is  the  development  of  a  non  —  frictional  contact 
between  the  sample  and  the  rigid  frame.  Other  experiments 
could  be  designed  to  study  the  time  —  dependent  response  of 
a  rock  sample  subjected  to  non  uniform  loading. 

The  responses  of  the  instruments  in  the  model  test  are 
sensitive  to  a  change  in  loading  rate  and  in  N-value  (ratio 
between  the  horizontal  and  the  vertical  stress).  Uniformity 
in  the  loading  conditions  depends  on  the  ability  of  the 
operator. The  next  step  in  the  Improvement  of  the  compression 
machine  should  be  the  design  of  an  automatic  and  reliable 
loading  system  where  the  N-value  could  be  constantly 
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monitored*  A  servo— con t ro 1  for  the  plane— strain  condition 
could  also  be  Included*  A  system  or  computer  code  allowing 
rapid  monitoring  of  the  N— value  would  be  useful  for 
experiments  in  the  near  future* 

In  the  stress  stages  preceeding  the  rupture  there  was 
no  clear  indication  of  the  future  local  i ns  tab  1 1 1 1 y • 
Therefore  the  time  —  dependent  aspect  of  the  rupture  was  not 
correctly  controlled  and  analysed  in  test  MC— 4 •  This  may  be 
the  most  Important  point  to  be  studied  in  further 
experiments*  It  is  advised  for  further  experiments  to  keep 
the  same  procedure  up  to  11  MPa  field  stress*  An  amount  of 
creep  for  the  tunnel  closure  in  the  order  of  0*2%  during  the 
first  day  at  this  stage  seems  to  preceed  rupture  but  does 
not  ensure  this  type  of  ultimate  behaviour*  The  next  stress 
Increments  should  only  be  1  .MPa  and  the  corresponding  creep 
tests  should  be  carried  on  for  a  two  day  period*  even  after 
the  eventual  initiation  of  rupture*  The  data  should  be 
analysed  and  termination  of  the  creep  or  rupture  process 
should  be  confirmed  by  v lew ing  the  strain  rate  -  time  plots* 
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APPENDIX  A1 


CONTINUUM  MECHANICS  RELATIONSHIPS  FOR  PLANE  STRAIN  CONDITIONS 


The  equations  of  solid  mechanics  under  plane  strain  conditions  are 
summarized  as  follows. 

For  plane  strain  conditions  the  longitudinal  strain  is  equal  to  zero 
and  the  value  of  the  stress  in  this  direction  is  given  by  : 


o33=  v(  oll+o22  ) 

->Equi  1  i  bri  urn  equation  :  (2  eq.) 

(1) 

3a.  . 

3j  +Fi  =  °. 

-^Strain-displacement  (3  eq.) 

,  3u.  3u. 

P  =  —  (  — 1  h — 3  \ 

eij  2  3j  3i  1 

(2) 

(3) 

->Hooke '  s  law  (3  eq . ) 

o .  .-X  6..  e,.  +  2G  e .  . 
ij  ij  kk  ij 

(4) 

A,G  Lame  constants  ,  6^  Kronecker  delta 


These  8  equations  are  sufficient  to  calculate  the  values 
of  the  8  unknown  : 


o 


i  j 


(3) 


(3) 


(2) 


For  convenience  the  solution  is 
displacements  or  stresses.  This  will 
di ti ons . 


usually  expressed  either  in  terms  of 
take  into  account  the  boundary  con- 


a)  Solution  in  terms  of  stresses 

The  compatibility  equation  is  given  by  . 


a2 


3i  3  j 


+ 


(5^ 


.  • 
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From  equations  2,4  and  5  the  following  relationship  between  the  stresses 
and  the  body  forces  yields  : 


3F. 

l 

3i 


For  a  two-dimensional  problem  the  Airy's  stress  function  (p  is  introduced 
as  follows  : 


3  cb 


The  body  forces  are  said  to  derive  from  a  potential  if  a  function 
V  may  be  found  such  that  the  following  equation  be  verified  : 


The  substitution  of  functions  cp  and  V  into  equation  6  leads  to  a 
governing  equation  for  the  problem,  called  biharmonic  equation  : 


1  -2v 

TTT 


(7) 


b)  Solution  in  terms  of  displacements 

From  equations  3  and  4  the  following  relationship  between  stress  and 
displacement  is  obtained  : 

3 u ,  ,  3u •  3u  . 

_Ui  +  G( — J.  4—1) 

3  k  3  j  3 i 

Substitution  of  cm.  into  the  equilibrium  equations  will  give  Navier's 
equation  :  J 


o.  ,=X  6.  . 
l  J  l  J 


3u .  .  3u  . 


GV2ui  +  (X+G)-gT  (gi-i  +  gj3)  +  F.  =  0 


(8) 


c)  Solution  in  terms  of  polar  coordinates 

The  stresses  are  transformed  from  polar  to  cartesian  coordinates  as 
follows  : 


a  =o  cos26+aQsin2e-T  ftsin20 
x  r  0  r0 

a  =o  sin20+ancos20+i  Asin20 
y  r  0  r0 

t  =(a  -aQ)sin0cos0+i  .(cos20-sin20) 
xy  r  0  r0 


i  ' 


,,  ..  ■  ■'  «.o  »*i 
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The  equilibrium  equations  become  : 


!!i  + 1  3Tre  +  !r^e  +  F 

3r  r  39  r  r 

1 _ 6  +  _ Hi  +  2  -i  +  F  = 

r  B6  3r  r  9 


=  0 


0 


(9) 


and  Hooke's  law  may  be  expressed  as  follows 
-£ -  (  CJ„  -  (~ )°a) 


l-v;  6' 


Ee=  1^1  (  °e  '  (T^r)or) 


1 


Yr0  Gt r6 


(ii; 


The  strain-displacement  relationships  become 
3u 

r 

er  3V~ 


y 


u  i  3uQ 
r  +  1  ;9 

Or  r  36 

.  liv  +^e 

r6  r  36  3r 


u. 


(12) 


and  Navier's  equations  reduce  to  : 


(A+2G)~  -  ^  j£  +  F..  =  0 


36 

<A+2G>7|  +  2GI?+  Fe  =  0 


where 


r=£  +c  -  — 

L  ~r  9  r 


1  3(rur) 


3r 


+  — 


1  8U0 

r  36 


,  3( ru  )  '3u 

_  1_  (  r  + _ t  ) 

u)  -  2r  (  3r  "Q 


36 


(10) 


Therefore  the  bi harmonic  equation  becomes  : 


^(g.ti|F^SrHS*+Flr  +  pS*J*-0 


(13) 
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Stresses  in  polar  coordinates  may  be  derived  from  the  Airy's  stress  func- 
-tion  by  the  following  expressions  : 

rr  =  1  ii  +  1  ^ 
r  r  9r  r2  90 z 


_  9  /  1_  3({)  \ 

Tr0  ’  3rl  r  96  j 


The  general  solution  of  the  biharmonic  equation  in  polar  coordinates 
is  given  by  the  following  expression  : 

(after  Obert  and  Duvall  (1967) 

<j>(r,0)  =  A°r2+B°r2log  r  +  C°  log  r  +  D°  .. 

+(E°r2  +  F°r2log  r  +  G°  log  r  +  H°)0 

+  (Axr  +  -1  +  Cir3  +Dir  log  r)sin6 

+  (Ejr  +  — *  +  Gir3  +-  h^r  log  r)cos0 

+  ( J i r  +  log  r)0sin0  +  ( L ! r  +  Mir  log  r)0cos0 

+  !  (A  rn  +  B  r'n  +  C  r2’n  +  D  r2+n)sin  n0 
2  n  n  n  n 

+  ?  (E  rn  +  F  r‘n  +  G  r2_n  +  H  r2+n)cos  n0 
2  n  n  n  n 


. 
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APPENDIX  A2 


ELASTO-PLASTIC  ANALYSIS  OF  THE  MODEL  TEST  RESULTS 


2.1  ASSUMPTIONS 

The  stress  conditions  are  isotropi c(N=l )  creating  an  axi symmetric 
stress  field  .  ( u Q=0 ) 

List  of  symbols 


o°  external  pressure 

o  radial  stress 

o0  tangential  stress 

oc  unconfined  compressive  strength 

E  Young's  modulus 

v  Poisson's  ratio 


r 

Lame  constants 


( 1 +v) ( 1 -2v ) 


a  radius  of  the  tunnel 

r  current  radius 

R  radius  of  the  plastic  zone* 


u  radial  displacement 
u^  tangential  displ acement (=0) 


e  radial  strain 
e*  tangential  strain 


angle  of  internal  friction 
C  cohesion  intercept 

m  coefficient  of  passive  earth  pressure  (  =  kp) 
a  dilation  coefficient 

List  of  Superscripts 


in  the  plastic  zone  :  plastic  part  of  ... 
in  the  elastic  zone  :  elastic  part  of  ... 
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2.2  Elastic  Relationships  for  the  Model  Test 


a)  Stress  distribution 


The  original  derivation  of  the  stress  distribution  around  a  circular 
hole  is  attributed  to  Kirsch.The  stress  field  obtained  by  using  Airy's 
stress  functions  is  found  in  numerous  books  (Obert  and  Duvall  1967). 


a 


v°0(1 


) 


5 


T 


re 


=  0 


o  +a  =2a° 

r  6 


(14) 


b)  Strain  distribution 


The  strains  are  obtained  from  the  stress  distribution  using 
Hooke's  law  (11)  : 


c  _  o°( 1+v)  / 

er  E  ^ 
_  =  g°(l+v)  ( 
V  E  ^ 


1-2V--2  ) 
r 

1  - 2v+”2  ) 
r 


(15) 


An  intact  plate  would  have  the  following  uniform  strain  distribution  if 
stressed  under  the  same  conditions  than  the  model  with  tunnel  : 

c  =  (1  -2v) 

c)  Displacement  field 

The  displacements  are  obtained  by  integrating  the  strain-displa 
-cement  relationships  (12)  : 

V°  ;  ur=a-(iiil  (d.2v)r  +  i2  )  (16) 

d)  Elastic  closure  of  the  tunnel_ 


The  closure  of  the  tunnel  is  defined  as  the  ratio  between  the 
displacement  at  the  tunnel  wall  and  the  original  radius  of  the-  tunnel . 
It  has  the  dimension  of  strain. The  closure  of  the  tunnel  in  the  test 
is  given  by  : 

,'V  2a°(1-v2)  (17) 

'r-  r=a  E 

t 

The  closure  of  a  real  tunnel  excavated  in  a  prestressed  uody  is  : 


155 


/ ur n.  _  2 q°  (1  +v) 

V  }  E 

Assuming  a  Poisson's  ratio  v=0.25,which  leads  to  the  so-called 
Poisson's  relation: 

A=G 

the  relationship  between  the  closure  of  the  tunnel  in  the  model  test 
and  the  closure  of  the  tunnel  in  the  field  is  given  by  : 

(remodel  =  |  (V)  ,  . . 

v  2  V  'reality 


2.3  El asti c-Perfectly  Plastic  Model 


a)  Coulomb  failure  criterion 


The  Coulomb  failure  criterion  relates  the  minimum  and  maximum  princi¬ 
pal  stresses  as  follows  ; 

a1=ma3  +ac 


where  :  o1,o3  are  the  maximum  and  minimum  principal  stress 

a  is  the  unconfined  compressive  strength 

c  1+sincJ)  .  cos$ 

m  =  T  vyva  and  a  =  2Cy  ■  ■ 

1  -  s  l  n  cf  c  1  -  s  i  n  <p 


In  the  case  of  the  model  test,  this  criterion  is  between  the 
radial  and  the  tangential  stress. 


VmV°c 

or  f=re-mar-a(,=0  (18) 

b)  Stress  distribution  in  the  Plastic  Zone 

The  combination  of  the  equilibrium  equation  (9)  and  the  failure 
criterion(l 8)  gives  the  following  differential  equation  : 


r 


ar(m-l )  +  ac 


(19) 


whose  solution  is  : 

VT-m  +  Cte) 

r  1  -m 


156 


The  boundary  conditions  give  the  value  of  the  constant  : 

1-m  a- 


1-m 
a  a 


Cte  =  a'-m(p.- JE_  )  =___ _ 

1-m  '  1-m 

i f  at  r=a  a  =p.=0  . 

r 

Then  the  stress  distribution  in  the  plastic  zone  is  given  by 

a 


_  _ c  /  / a \ 1 -m 

Vm  (  (7>  -]) 


O 


cr, 


c  /  #a\l-m 


fniff  ' '  -1) 


(20) 


c)  Radius  of  the  Plastic  Zone  R 

In  order  to  estimate  the  radius  R  of  the  boundary  between  the  plastic 
zone  and  the  elastic  zone,  the  radial  stress  is  assumed  to  be  continuous 
across  this  boundary.  The  value  of  the  stress  from  the  plastic  side  of 
the  boundary  is  given  by  : 

vVr=R  m-1  1  lRj  ) 

The  value  of  the  stress  from  the  elastic  side  of  the  boundary  is  : 


(or)r=R  =m4r(  2a°-°c) 


The  latter  equation  is  obtained  after  the  following  relationships 

KW  m(or)r=R  +  °c 
^el.  +  «°r)el.  ■2o° 


The  radius  R  is  determined  from  the  elimination  of  (ar)r=R  between  the 
above  equations. 


1 

r  ,  2  o°  (m-1  )+ar 
a  ( Qn+O  ac 


(21) 


2  •  *  I 


■ 

- 
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d)  Stresses  in  the  elastic  region  surrounding  the  plastic  zone 


The  stress  distribution  in  the  elastic  region  is  identical  to  the 
one  encountered  in  a  body  with  a  R-diameter  tunnel, which  is  supported 
by  a  uniform  pressure  a  under  the  same  field  stress  conditions. 


At  the  elastic-plastic  boundary,  the  stresses  are  given  by  : 

ac  ,  /R.m-l  n  x 
arR  m-1  ^  ” 

(22) 


c  /  / R \ m~ 1  i  \ 

°0R  =nPT  (  m(T>  -  1  > 


°0R+arR=^a' 


(23) 


and  in  the  elastic  region  the  stresses  are  given  by  : 

on  R%.  R2 
O  =0U(l-'H)+0  n  — 2 
r  r  rR  r 

ae=o°(Hp)-orR-p  (24) 


Hooke's  law  (11)  and  the  integration  of  the  strain-displacement 
relationships(12)  give  the  strain  and  displacement  fields. 


o°( 1+v) 

V  E 

_o°( 1+v) 
£e  E 


a°(l+v) 

U  = - r - 

r  E 


(l-2v  +  (JR-  l)j£  ) 

(1  -2v  -  (  -Jr  -  l)p  ) 

((l-2v)r  -  (-^  .  ■,)£') 


(25) 


The  substitution  of  equations  (22)  into  (23)  and  (25)  will  give  the 
exact  expressions  for  the  stress  .strain , di spl acement  fields  in  the 

elastic  zone. 


•  ' 


. 


158 


e)  Strains  in  the  plastic  zone 

The  plastic  strain  increments  are  obtained  by  application  of  the 
associated  flow  rule  definition: 


•p 

r  = 


3f 


da 


=  -Am 


ep  = 

e0 


*P 

£z  = 


v 3f  -  \ 

H 

-  n 

A9a  "  0 

z 


where  f  is  the  failure  criterion  given  by  equation  (18). 

It  follows  that  : 

ep  +  meP  =0 
eP  =  0 

The  integration  of  these  equations  yields  : 

er  +  m£e  ■  g^r) 
eP  =  g2(r) 

Initially  the  plastic  strains  are  equal  to  zero 
g-j  (r )  =  g2(r)  =  0 

Hence  the  associated  flow  rule  which  causes  the  plastic  increments  of 
strain  to  be  normal  to  the  Coulomb  failure  surface  may  be  summarized 
as  : 

ep  +  mep  =  0 
r  0 

ep  =  0 

In  order  to  control  the  plastic  flow,  the  following  relationship  will 
be  used  between  the  plastic  strain  values; 


ep  +  aeP  =0 


eP  =  0 


(26) 


The  parameter  a  isa  measure  of  the  dilation  occuring  during  the 
plastic  flow. If  there  is  no  plastic  volume  change ,a=l  ;i f  a=m,the 
flow  rule  is  associated  to  the  Coulomb  failure  criterion. 
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The  total  strain  in  the  plastic  zone  is  the  sum  of  the  plastic  and 
elastic  strains.  Plane  strain  conditions  have  to  be  satisfied  in  the 
plastic  region.  As  tbecplastic  longitudinal  strain  is  equal  to  0  (see 
26)  ,  the  elastic  longitudinal  strain  has  toffee-equal  tolO. Hence  the 
stress  az  in  the  longitudinal  direction  is  : 


a 


z 


VG 


c 


m-1 


(27) 


The  elastic  strains  in  the  plastic  zone  are  determined  from  the 
combination  of  Hooke's  law  (11)  and  plane  strain  condition  : 


Er=  7gT  CO-On+l  )v)ar  - 

f(m-(m+l  )v)ar  +  (l-v)ac] 


(28) 


The  compatibility  equation  for  the  total  strains  is  : 


de^ 

dr 


+  e9  -£r  =  0 


(29) 


From  this  relationship  the  following  differential  equation  is  obtained 
for  the  plastic  strains  : 


de 

d7 


§■  +  £^l+a)=-^(l-v)(m+l)(f)m'1 


'0 


(30) 


The  solution  of  this  equation  is  (with  e^(R)-O)  : 

p  1-v  m+1  a  /rxm-1  r /R^m+ot  Tn 

£e  in+oT-  c  [(7}  -1] 

(31) 

P  P 

er  ae0 

The  total  strain  distribution  in  the  plastic  zone  is  determined  from  the 
combination  of  equations  31  and  28  : 

°c  r,r*m-l  ,m(l-v)-v.M  ..\m+b  fR\m+a.  ■,  u  +  2vJ_  n 

(  m-l  (1'VWV  m  m-1  J 


_°c  r,r>m-l  ,lwnv^i  ,  \"ltL((B.)m+0l-i))  +  ] 

¥  (T7T^  01(1  VWV  n  m-1 


(32) 


■ 


'  t 
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It  follows  from  the  integration  of  the  strain-displacement  relation¬ 
ships  (12)  : 


ur=ree 


The  closure  at  r=a  is  given  by  : 


Ua_  Qc(1  ~V^n  m+1 , ,  RNm+g  ,  n 
?  2G  Ll  m+aUa;  "UJ 


For  a=l  this  expression  yields  : 

ua  .  °c(1  "v) rR-iHi+1 
a  2G  LaJ 


(33) 

(34) 


f)  Comparaison  between  the  el asto-pl asti c  and  the  elastic  closure 


The  closure  may  be  normalized  by  dividing  the  "el asto-pl asti c" 
closure  by  the  elastic  closure  :  (34/17) 


For  cel 


For  a/1 


^u/a^e.pl.  .  °c  ,R,m+l 

(u/a)el .  "  2o°  a 


(u/a)e.pl. 

(u/a)el 


[1  +!5ll(  (-)m+ct.-| )] 

L  m+avva;  ;J 


(35) 

(36) 


This  ratio  is  independent  of  the  elastic  constants. 


2.4  Elastic-brittle  plastic  material 


a)  Failure  criterion 

The  same  Coulomb  failure  criterion  still  applies,  but  now  the  stresses 
in  the  plastic  zone  are  related  by  the  following  "plasticity"  criterion: 

ae=mVsac 

0<s  <1  (37) 

The  parameter  s  measure  the  intensity  of  the  strength  loss  occuring 
immediately  after  the  strength  failure. The  strength  after  the  failure 
is  reduced  to  s-times  the  value  it  was  before  failure. 


1 
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The  solution  is  similar  to  the  one  described  in  the  previous  section 
for  an  elastic  perfectly  plastic  material .Only  some  equations  are  af- 
-fected  by  the  "plasticity"  criterion  and  are  given  thereafter. 


b)  Stresses  in  the  Plastic  Regi 


on 


The  stresses  in  the  plastic  region  are  now  given 
expressions  : 

Vsrf  ft?)1 'm  -  i] 

Vl=f  K  7)1_m  - 1] 


by  the  following 


(38) 


The  radial  stress  is  continuous  but  the  tangential  stress  experiences 
a  stress  drop  of  (l-s)cr  at  the  elastic-plastic  boundary. 

c)  Radius  of  the  Plastic  Region 


The  radius  of  the  boundary  of  the  plastic  region  is  given  by  : 


£  =  r.  2 
a  L(m+1 )s 


+  m(s-1 )+l+s 
2 


1 

) ]m-l  (39) 


If  the  material  looses  all  its  cohesion  after  failure  (s=0),  this 
radius  increases  to  infinity  and  the  tunnel  col  lapses .Actually ,  a  tunnel 
in  a  purely  frictional  material  would  collapse  without  a  support 
pressure . 


d)  Strains  in  the  Plastic  Region 


The  derivation  leads  to  a  differential  equation  which  is  identical  to 
equation  30 . However , because  of  the  boundary  conditions,  the  solution  is 
di fferent . 

The  total  strains  have  to  be  continuous  across  the  elastic-plastic 
boundary .Because  of  the  brittle  character  of  the  material,  the  stresses 
in  the. pi  as  tic  regi on, given  by  equations  38  are  smaller  than  in  the  case 
of  a  perfectly  plastic  material .Hence  the  elastic  strains  which  are  rela- 
-ted  to  the  stresses  by  Hooke's  law  are  also  smaller. Then  for  the  conti- 
-nuity  of  the  total  strains  the  plastic  tangential  strain  has  the  follo- 
-wing  finite  value  at  the  elastic-plastic  boundary  : 


e 


P 

0R 


2§  0-v)(l-s) 


(40) 


. 


* 
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The  strain  field  in  the  plastic  zone  is  : 


o 


£d=  2GL  'a 


m+a '  r* 


m-1 


...(l-v)(l-s) (~)a+1  ] 


_  crwrNm-l  ,1-mv-v 


v  #s(7r-(- 


m-1 


.a(1-v)!5+l((R)"W«  _U)  +  illy- !,) 

m+a v  v  r  ‘ ™  1 


m-1 


..  -a(l-v)(l-s)(ji)“+1] 


(41) 


The  displacement  field  in  the  plastic  region  is  determined  from  the 
integration  of  the  following  strain-displacement  relationships  : 


In  this  case  the  closure  of  the  tunnel  is  : 


~  =  ^fO-v)[s(l+Eli((f)m+a-l))  +  (l-s)(f)a+1]  (42) 

And  the  corresponding  normalized  closure  is  : 

~ ^  ■  J(^(l^((f)m+a-D)  +  0-s)(f)a+1]  (43) 

el . 


2.5  Elastic  -  strain  weakening  material 


This  model  assumes  that  the  strength  loss  after  peak  is 
strain-dependent .The  "pi astic"cri terion  is  then  strain-depen¬ 
dent  ;it  is  : 

ae=mar+oc(i— ' ’)  ^e"£r^<Y<  00  <44) 

The  solution  is  more  complicated  but  the  system  is  isostatic. 

The  crefficient  y  is  a  measure  of  the  after  peak  slope  of  the 
stress-strain  curve  of  the  strain- weakening  model 


The  constitutive  relationships  are  : 


equilibrium  equation 


do 


oQ-o 
0  r 


dr  r 

plastic  criterion 

a  =mo  +o  (1 — ■ — -) 
0  r  cv  y 

Hooke's  law 

e  1  -v2 


V  E 


(O  Q  ) 

0  1 -v  r  ' 
e_l -v2 v  x 
er~  E  ^ar~l -v  V 


flow  rule 
e^+ae^O 
total  strain 


e 

£e  e  £e 


P 

£  =£r+£ 

r  r 


compatibility  equation 


de 


0 


r  -7-~  +  e.-e  =  0 
dr  0  r 


e 

r 


This  system  is  isostatic:8  unknown  and  8  equations  .  However , the 
solution  is  difficult  as  the  derivation  ends  up  to  a  differential 
equation  of  the  second  order  with  variable  coeffi cients ,and  needs 
numerical  treatment. 


■ 


. 


'  % 
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APPENDIX  A3 


A3.1  Correction  of  Data  for  a  Non-Isotropic  Loading 


The  method  which  is  described  in  the  fol lowing, shows  how  to 
find  an  elastic  reference  value  to  calculate  the  normalized 
closure  if  the  N-value  (N=a?/a°)  varies  slightly  around  1 .This 
method  does  not  overcome  the  problem  of  the  non-validity  of  the 
normalized  closure  concept  if  N  is  different  from  1, but  permits 
to  get  more  consistent  results  for  the  experimental  normalized 
closure . 

The  regular  method, as  described  in  Chapter  5, takes  as  a  elas- 
-tic  reference  line  a  straight  line  through  the  two  first  points 
of  the  long  term  loading  curve. The  N-value  is  assumed  to  be  cons- 
-tant  during  the  test .Actual ly  at  each  experimental  point  cor- 
-responds  a  particular  value  of  the  parameter  N,calculated  from 
the  load-cells  measurements. 

At  each  point  the  secant  modulus  of  elasticity  is  theoretical¬ 
ly  given  by  : 


where  X  is  a  function  of  the  ratio  N  only  and  depends  on  the 

directionaof  the  measurement  with  respect  to  the  external  stress 

axes. If  a  is  the  angle  between  the  direction  of  the  closure 

measuring  gauge  and  the  direction  of  the  maximum  external  stress, 

X  is  given  by  : 
a  3  J 


(2) 


Considering  the  first  two  points  (1&2)  of  the  experimental 
curve, the  knowledge  of  the  corresponding  N-values  will  allow 
to  calculate  their  appropriate  X  -values  ( X  i  2 ) .  If  SL1 
and  SL2  are  the  secant  moduli  of  the  curve  at  points  1  and  2 
respecti vely .they  are  : 


bL  l  _  a 

SL2  '  X  2 

a 


SL1 


If  e1  and  c2  are  the  values  of  the  total  abscissa  of  the  points 
1  and  2, and  if  o1  and  a2are  the  values  of  the  stress  at  points  1 
and  2,  SL1  and  SL2  may  be  determined  as  a  function  of  the  strain 


- 


' 


- 
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corresponding  to  the  closure  of  the  cracks  e 


cl 


e  -e 


SL1  =- 


cl 


E2-£ 


C 


SL2  = 


cl 


(4) 


o 


From  the  equations  3  and  4  it  is  possible  to  get  the  expression 
for  the  crack  closure  strain  : 


£cl" 


E 1 X  2<32-E2X  lO1 

a  a 

"  X  2o2-X  i o1 

a  a 


(5) 


The  value  of  the  elastic  parameter  can  be  also  evaluated  : 


AJ.(l-v2)  - 

H  E  X  2o"-X  j o 1 


e2  -e  1 

7T 


(6) 


CL 


a' 


For  the  next  points  the  value  of  N  will  give  the  X  -value, which 
will  allow  to  calculate  the  elastic  secant  slope  by  multi- 
-plying  by  the  A-value  independent  of  the  N-parameter.  : 


SL  =  A  X 


a 


The  elastic  reference  strain  is  then  known. The  substraction  from 
the  total  experimental  el asto-pl asti c  strain  of  the  crack  closure 
strain  e  will  give  the  el  asto-pl  asti  c  closure  to  be  used  to 
calculate  the  normalized  closure  : 

£t'ecl 


A  X  o 

a 


norm.clos . 


(7) 


’  rt  0  ,  #  l 
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A  3.2  Influence  of  the  Intermediate  Principal  Stress 


In  the  model  test, the  longitudinal  stress  o  is  assumed  to  be  a 
function  of  the  other  stresses, as  the  model  sampfe  is  loaded  under 
plane  strain  conditions. 

o  =v(o  +o  )  (1 ) 

z  r  6 

This  relation  is  also  valid  in  the  plastic  zone. But  in  this  zone  the 
radial  and  tangential  stress  are  related  by  a  failure  criterion; 


o  =mo  +o  (2) 

6  r  c 

The  combination  of  these  equations  will  give  the  following  expres¬ 
sion  for  the  longitudinal  stress: 


a  =v(o  (l+m)+o  ) 
z  r  c 


(3) 


When  loading  proceeds, the  elastic-plastic  radius  R  attains  a  value  R 
such  that  : 


or 


°zR_0rR 


VG 


G  = 


1 -v(l +m) 


(4) 


In  the  elastic  zone  the  relation  (1)  is  equivalent  to  : 

g  =2vo  (5) 

Z  0 

where  Gn  is  the  external  stress. In  the  elastic  zone, the  longitudi¬ 
nal  stress  is  constant. To  the  critical  value  o  corresponds  by  the 
relation  (5)  a  value  of  the  external  stress  oQ  T 

ac  (6) 

°o  2  ( 1  -  ( 1  +m )  v ) 

At  the  upper  limit  of  the  loading  range  of  the  compression 
machine, the  external  stress  is  equal  approximately  to  17  MPa. In 
order  to  have  the  longitudinal  stress  as  the  intermediate  principal 
stress, in  loading  the  sample  ,the  Poisson's  ratio  has  to  be  higher 
than  a  critical  value  given  by  the  relation  (6). 

Numerical  application  :  m=3  ( 4>= 30 0  )  ,C=4MPa,oc=13.87  MPa  : 

d0=17  MPa  +  v=0 .1 48 

Poisson's  ratio  may  be  expected  to  be  bigger  than  0.15, and 


2  '  • 


■j 
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the  failure  criterion  (2)  will  be  assumed  to  be  valid  throughout  the 
test . 


This  section  is  based  on  a  paper  by  Florence  and  Schwer(l 978) . 
The  authors  explained  the  formation  of  three  different  regions  around 
an  opening;ri ght  around  the  cavity, a  first  annular  region  obeys  the 
criterion  (2);away  from  the  tunnel  the  longitudinal  stress  is  the 
minimum  principal  stress  and  the  failure  criterion  is  then  : 

o0=moz+oc  (7) 

The  third  zone  is  located  between  the  two  other  zones, as  the  longi¬ 
tudinal  and  the  radial  stress  are  equal  in  this  third  zone. The 
stressfijobey  both  failure  criteria  (2)  and  (7). These  three  plastic 
regions  develop  one  after  the  other  .  ,and  in  the  model  test, it  is 
assumed  that  only  the  first  zone  develops  within  the  loading  range 
of  the  compression  machine. 


* 

. 
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A  3.3  Volume  Change 


The  total  volume  affected  by  plastic  deformations  around  the 
tunnel  is  : 

V=tt(  (-^)2  -  1 )  (1 ) 

At  each  point  of  this  body  under  plane  strain  condition,  the 
volumetric  strain  e(r)  is  given  by  : 

e(r)  *  £r+£0  (ez=0)  (2) 

The  change  in  volume  is  equal  to  the  integration  of  the  ' 
volumetric  strain  over  the  whole  yielding  volume. 

AV  =  /2lT  /R  e(r)  r  de  dr  (3) 

6=0  r=a 

In  the  plastic  zone  the  volumetric  strain  is  given  by  : 

e(r)  =  £^+£0+  £^+£q  =  £r+£0  (4) 

elastic  pi  as t i c 


Volumetric  strain 


a)  Elastic  part 


Equation  28  in  Appendix  A2  gives  the  elastic  strain  distribu¬ 
tion  in  the  plastic  zone.  From  these  expressions  it  follows  : 


e(r) 


el . 


oc(l-2v) 
2G(m-l ) 


[;I)m-1(m+l)  .  2  ] 


(5) 


This  expression  may  also  be  written  as  follow 

m+1 


eWe1--  (rV)ext  + 


l-2v  °c 


m-1 


(6) 


where  ~  .  is  the  constant  volumetric  strain  in  the  elastic 

V  ext 

zone  and  also  the  external  volume  change  due  to  the  displacement 
of  the  external  boundaries. 


b)  Plastic  part 

From  the  flow  rule  in  the  plastic  zone  (equation  26  -  A2),  the 
volumetric  strain  in  the  plastic  zone  is  given  by  : 

e(r)p=eP+eP  =  (l-a)cP 

for  the  plastic  part  of  the  strains. 


(7) 


. 
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The  substitution  in  equation  7  of  the  value  of  ej?  given  by  the 
equation  A2-31  gives  the  following  expression  : 

e(r)P=(1.a)^-^E±i(^-l[(f)m+“-l]  (8) 

Integrati on_of_the  volumetric  strain 


The  integration  of  equation  6  and  7  according  to  the  expression 
given  by  equation  3  will  give  the  change  in  volume  occuring  with 
plastic  deformati ons . Thi s  change  in  volume  may  be  divided  by  the 
volume  of  the  plastic  zone  (  equation  1  )  to  obtain  an  "average 
volume  change". 

In  the  fol 1 owi ng, expressions  for  the  average  volume  change 
are  given  (for  the  elastic  and  plastic  component  of  the  strain). 

The  combination  of  these  expressions  will  give  the  total  volume 
change.  This  is  illustrated  in  Figure  5.1  in  Chapter  5. 


-elastic  : 


/  AV  x 

ee  V  ext 


°rR  R2 


(9) 


where  o  D  is  the  radial  stress  at  the  boundary  el asto/pl asti c 
r  K  . .  i 

(equation  A2  -  22)  and  r, —  the  external  volume  change  given 


'  ^  '  y  Gxt 

by  the  following  expression: 

/ AV x  2  g°(l+v) (l-2v) 

V  ext  E 


(10) 


-plastic  : 


-i  o 

1  -V  c 


VS*  25  ^^)(f)m+1-(m+1)(I)^+1-a)  (11) 

r  —2-  I 


=  e  +  e 
e  p 


e 


(12) 
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APPENDIX  A4 


DETERMINATION  OF  ELASTIC  PARAMETERS  IN  THE  MODEL  TEST 


Introduction 


The  purpose  of  this  appendix  is  to  present  two  methods  which 
have  been  developed  to  determine  the  elastic  parameters  of  the  coal 
specimen  used  in  the  model  test.As  the  loading  is  not  uniaxial,  the 
elastic  response  of  the  instruments  is  a  function  of  two  elastic 
parameters.The  methods  presented  hereafter  compare  the  responses  of 
two  different  instruments  in  order  to  eliminate  one  of  these 
unknown  elastic  parameters. 


Instrument  responses 
-  Internal  Strain  Gauge  : 


o° ( 1  -2v)  a0  a2  r_  a2 

e  =  2G -  '  2G  d^  C“^Fr 


(1) 


where  -  c  is  the  engineering  strain  resulting  from  the  analysis 
of  the  data  by  the  computer  codes  and  appearing  on  the 
pi  ots 

-  d  is  the  distance  from  the  centre  of  the  tunnel  to  the 

middle  of  the  strain  gauge 

-  e  is  half  the  length  of  the  strain  gauge 

-  a°is  the  applied  external  stress  assumed  hydrostatic. 

In  the  remainder  of  this  appendix  the  internal  strain  gauge  will 
be  characteri zed  by  a  constant  C  which  is  determined  with  the 
location  of  the  gauge. 


-  Tunnel  closure  : 

2  o°(l-v2) 
e  E 


(2) 


-  External  Displacement 


£  = 


o°(1+v) 


( 1  - 2v ) 


a 


(1+v)  a2 


(3) 


where  b  is  the  distance  from  the  centre  of  the  tunnel  to  the 
boundary  where  the  displacement  is  measured. 


t 
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Fi rst_method 

This  method  comparesthe  responses  of  the  same  internal  strain 
gauge  before  and  after  the  tunnel  is  drilled. 

-Assumpti ons 

*  e1  is  the  response  of  the  gauge  before  the  tunnel  is  drilled 

*  e2  is  the  response  after  the  tunnel  has  been  drilled. 

*  the  loading  rates  are  similar  for  the  two  loadings. 

*  no  yielding  occurs 

*  the  elastic  properties  are  not  significantly  affected  by  the 
loading  history  after  the  specimen  has  experienced  several 
previous  loadings. 

*  the  specimen  is  homogeneous 

-  Derivation 

If  the  function  A (o)  is  supposed  to  be  : 


(4) 


the  responses  c1  and  e2  may  be  written  as  follows  : 

e1  =f  A(o° )  (l-2v) 
e2  =  A (a0)  ( 1 -2y-C) 


(5) 


The  elimination  of  the  function  A(a°)  between  the  two  equations  5 
leads  to  an  expression  of  the  Poisson's  ratio  v  : 


with  -l<v<0.5 

This  method  is  applied  to  each  of  the  16  internal  strain  gauges. 
Second  method 

This  method  compares  the  external  deformation  to  the  tunnel 
closure  in  the  same  direction  during  a  loading  of  the  specimen. 

-  Assumptions 

The  stresses  in  the  model  test  are  equal  to  the  stresses  in 
an  infinite  medium  around  a  tunnel  .mulitlied  by  a  coefficient 
function  of  the  geometry  of  the  specimen  used  in  the  model  test. 


where  b  is  half  the  width  of  tie  specimen. 
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But  as  shown  in  Chapter  5,  the  change  in  stress  is  not  significant 
and  is  neglected  in  this  analysis. 


The  restricting  assumptions  resulting  in  the  first  method 
from  the  use  of  responses  during  two  different  loadings  are  over¬ 
come  .But  instead  of  16  measurements  the  actual  method  gives  only 
2  measurements , one  in  the  vertical  and  one  in  the  horizontal 
di recti  on . 


-  Derivation 


The  overall  deformation  of  the  sample  (  e1  )  is  given  by 
equation  3  of  this  Appendix. The  tunnel  closure  (  e2  )  is  given 
by  equation  4. The  elimination  of  the  same  function  A(a°)  between 
these  two  expressions  leads  toan  expression  for  Poisson's  ratio 


v 


£2  jO+p) 


(?) 


-  c 


1 


Results 


The  results  are  too  erratic  to  be  presented  here  in  detai 1 .However , 
some  remarks  may  be  made  about  them  : 

-  Poisson's  ratio  is  increasing  during  loading  up  to  a  value 
around  0.2. 

-  The  knowledge  of  Poisson's  ratio  and  of  the  instrument 
responses  make  possible  the  determination  of  any  other  elastic 
parameter . 

-  Young's  and  Bulk  moduli  are  also  increasing  during  loading 
up  to  values  around  2200  MPa  and  1000  MPa  respectively. 


Con cl  us i on 


The  only  valuable  conclusion  of  this  analysis  is  that  it  con- 
-firmed  a  value  of  Poisson's  ratio  lower  or  equal  to  0.2  in  the 
planes  parallel  to  the  bedding  planes. 
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STRESS  -  STRAIN  CURVES 
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